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ABSTRACT 


ENGINEL is a computer program which was developed to",generate the design 
and off-design performance of a s.ingle rotor turbojet engine' with- or without 
afterburning using a cycle match procedure. It is capable 'of producing 
engine performance over a wide range of altitudes and Mach numbers. The 
flexibility, of operating with a variable geometry turbine, for improved 
off -design fuel consumption or with a fixed geometry turbine as in conven- 
tional turbojets, has been incorporated. In addition, the option of 
generating engine performance with JP4, liquid hydrogen or methane as a fuel 
is provided. 



FOREWARD 


This document presents in users format, with a brief description of 
computational procedures and examples, a Single Rotor Turbojet Engine Cycle 
Match Performance Program designated ENQINEL, developed by Mr. Arvid L. Keith, 
Jr., of the NASA Langley Research Center. 

The Vought Corporation Hampton Technical Center provides technical 
support to the SCAR Program Officer of the Aeronautical Systems Division under 
Contract No. NASI -13500. The work was performed under the LaRC Project 
Manager, Mr. Cornelius Driver, Head of the Vehicle Integration Branch, Aero- 
nautical Systems Division. This report was prepared by Mr. W. A. Lovell, 
under the direction of Mr. C. W. Pearce, the Vought Hampton Advanced Aircraft 
Technology Projects Supervisor. The report was edited by Mr. Emanuel Boxer, 
of the Vehicle Integration Branch, LaRC, and Mr. B. F. Downie of Vought 
Hampton. 
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SUMMARY 


The computer program ENGINEL was developed to provide engine performance 
for a single spool turbojet engine with or without afterburning, using a 
cycle match procedure. ENGINEL will generate, as specified by the user, 
engine design point performance at sea level static standard day conditions 
and will both develop and retain the parameters required to generate off 
design performance. It is then capable of producing engine performance over 
a wide range of altitudes and Mach numbers as requested by the user. 

In the cycle match computational , procedure ENGINEL utilizes a non- 
dimensional compressor map, constant efficiencies for the combustor, turbine, 
afterburner and a constant thrust coefficient for the convergent-divergent 
nozzle. A turbine map is not used, instead the cycle match is .accompl ished 
by flow matching using a constant turbine flow function, determined in the 
design point calculation. 

ENGINEL has the option, exercised by user request, of generating 
performance with a variable geometry turbine, for improved off-design fuel 
consumption, or a fixed geometry turbine as in conventional turbojet engines. 
In addition, the option of generating engine performance with JP4, liquid 
hydrogen or methane as a fuel is provided. 

Included as Appendix A to this document is a complete FORTRAN IV 
listing of ENGINEL. Sample problems, with input explanation, showing the 
various printed output options along with the normal engine control options 
are provided in Appendix B. 


INTRODUCTION 


The computer program described in this document provides engine 
performance for a single spool turbojet engine with and without afterburning. 
This program, ENGINEL, was developed by Mr. Arvid L. Keith, Jr., of the NASA 
Langley Research Center (LaRC), to generate the engine performance presented 
in the Reference 1 document. It was originally programmed by Mrs. Francis 
T. Meissner, also of NASA (LaRC), for the LaRC-CDC 6400-6600 computer 
complex using the ICOPS system. 

Since its inception, this program, has proved to be an exceptionally 
useful tool. It has been used extensively in , the study of Advanced Super- 
sonic Transport type aircraft. During the course of its use, it has been 
modified by Mr. W. A. Lovell of the Vought Corporation Hampton Technical 
Center to be operational on the present LaRC-CDC Cyber computer complex 
using- the NOS system. It has also been modified several times by Mr. Lovell 
to provide the additional engine performance required for engine/aircraft 
acousitc studies. A complete listing of the program is provided in 
Appendix A. 
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The program is capable of four modes of operation based on user 
selection as follows: 

° Fixed geometry turbine with or without afterburning. 

° Variable geometry turbine with or without afterburning. 

With any of these modes of operation the user may select full power operation 
only or full and part power operation. This is explained in more detail 
elsewhere in this document and is shown by example in Appendix B. 
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.sea level static conditions 
compressor surge 
free stream station 
compressor inlet 
compressor outlet 
exhaust nozzle throat 

COMPUTATIONAL PROCEDURE 


The engine performance computational procedure described below has been 
extracted from Reference 1, rewritten, and expanded so as to provide a basic 
explanation of the methods used in the ENGINEL program. 

This computer program utilizes a cycle match procedure to generate the 
performance of a single rotor turbojet engine. The program computes the 
design point cycle performance by conducting progressive, steprby-step 
thermodynamic calculations from the undisturbed free-stream, to the air 
intake, through each of the engine components and finally through the 
exhaust nozzle. In the cycle match iteration, a turbine map is not used, 
instead it is assumed that the v/ork output is available at a constant 
efficiency (input). The corrected flow parameter, determined from the 
design point calculation at the turbine inlet, is, therefore, retained and 
used to define, by flow matching, all possible operating points of a 
given design cycle. The nondissociated thermodynamic properties of gases 
from Reference 2 are used in the cycle calculation. 

The nondimensional ized compressor map characteristics, shown in Figure 1 
are typical of a single rotor supersonic turbojet engine compressor. These 
are input to the program in its present form, however, may oe replaced 
by another compressor should it be desired. Nondimensional ization is 
accomplished as defined in the input parameter identification list. 

The nondimensional ized compressor map is constructed so that the surge 
line pressure ratio (TABPT2 versus TABNST - Figure la) is defined from the 
design pressure ratio, PT2T1D (input) and the design surge margin, Y (input). 
The X value (input) sets the lower limit of pressure ratio for all rotational 
speed lines (TABNST) relative to the surge pressure ratio of each speed 
line. With the PT2T1D, Y and X parameters, a design point compressor 
pressure ratio function (Z) is formed. This Z function is then used to 
establish the compressor design point shown by the target symbol (on 
Figures lb and Ic). 

For the usual case, Y = 1.1, the design surge pressure ratio margin is 
ten percent greater than the operating pressure ratio and X = 0.1, the lower 
pressure ratio limit of the map is one tenth of the local surge pressure 
ratio. The purpose of the X parameter is to allow the intervals of the local 
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values of pressure ratio (TABZ, 0 to 1,0) to be squeezed together, for 
maximum accuracy in matching or spread apart, to insure that the matching 
values (as determined by the compressor map) of pressure ratio and corrected 
flow can be obtained. 

Off-design compressor operating points are established by iteration 
along constant corrected speed lines until the corrected flow at the turbine 
inlet is matched with the value determined in the design point calculation. 
Off-design corrected speed is defined by Mach number, altitude and the engine 
power level. 

Definition of off-design operating points by constant corrected turbine 
inlet flow was selected as a reasonable simplification to the computational 
procedure. The alternate method would be to include a map of the turbine 
work and flow characteristics. Variables which compromise the .choice of 
constant corrected turbine inlet flow as a matching parameter are: 

1. Turbine inlet Mach number less than unity. 

2. Unknown effective turbine inlet flow area. 

3. Variations in the ratio of specific heats and the specific gas 
constant, due to changes in turbine-inlet temperature and gas 
constituency, from the design value. 

Of the three effects, only item (3) was considered to be significant. 
Examination of the effects of the specific heat ratio and the specific gas 
constant when varied from their design values indicated a maximum variation 
in the corrected compressor inlet flow of less than two percent of the design 
value, when the turbine inlet temperature is varied from high values at 
design to very low values at off-design. The variation in compressor corrected 
inlet flow with Mach number and associated altitude resulting from the 
selected flow match procedure is shown on Figure 2 for a wide range of 
conditions. 

Control for the fixed geometry turbine mode of operation is accomplished 
by giving values to two and solving for the third of the following parameters: 

1. Compressor physical rotational speed. 

2. Turbine inlet temperature. 

3. Exhaust nozzle throat area. 

Maximum power at all off-design conditions is determined by setting the 
compressor rotational speed at 100 percent of design, turbine inlet 
temperature at the design value and computing the required nozzle throat 
area. Partial power is defined by a schedule of the ratio of nozzle throat 
area (Figure 3) as a function of compressor rotational speed. This manner 
of engine control requires the computer program to iterate to satisfy flow 
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matching, at both the turbine i.nlet and nozzle throat by varying turbine 
inlet temperature. • • 

Control for the variable geometry turbine mode of operation is initially 
the same as for the fixed geometry turbine in that the maximum power point 
is computed in the same manner. On completion of this solution, the compu- 
tational procedure differs in that the compressor map operating point is 
retained in the computer memory. Part power operating points are produced 
by reducing the turbine inlet temperature to provide lower than maximum 
power, while maintaining the compressor map operating point constant. The 
turbine inlet area is then reduced to match the design corrected turbine 
inlet flow and the exhaust nozzle throat area required, to pass the internal' 
flow, is calculated. The program recognizes a maximum variation in turbine 
area by the value assigned to the variable geometry turbine control parameter 
(A4A4DN). To solve for this limit, iteration of turbine inlet temperature 
and area is conducted until the limiting A4A4DN input value is obtained. 

The inlet total pressure recovery schedule as a function of Mach number, 
presently input to the program is shown on Figure 4. This schedule is 
considered typical of fully variable internal -external compression inlets, 
however, it may be replaced with any other suitable schedule if so desired. 

The exhaust nozzle, incorporated in the program, is a fully variable 
convergent-divergent nozzle, that operates in the fully expanded mode for all 
operating conditions. The nozzle schedule for controlling part-power engine 
operation (Figure 3) is typical of turbojet engines designed for supersonic 
speeds. 

. Current FAA regulations require jet engines and aircraft to meet specific 
noise requirements. Advanced turbojet engines, while operating at maximum 
power do not meet these regulations without the aid of jet noise suppression. 
The performance characteristics of an advanced jet suppressor have been 
extracted from Reference 3 and incorporated in the program. These charac- 
teristics are shown graphically on, Figure 5, with jet suppressor effective- 
ness (dB) as a function of exhaust jet velocity (ft/sec). 

A complete listing of the ENGINEL computer program is provided in 
Appendix A. Provided in Appendix B are several sample problems with the 
required input data and examples of the three printed output formats. 


DISCUSSION 


The ENGINEL computer program provides the performance of single rotor 
turbojet engines, as defined by the user selection of sea level static design 
component performance. It requires a storage capacity of approximately 
60,000 octal words and computation times as follows: 
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Time to compile « 5.26, seconds 

Run time for a design point = .13 seconds 

Run time per data point = .03 seconds 

With this program engine performance can be generated for wide ranges of Mach 
number, altitude and engine power setting with and without afterburning,. 

Incorporated, in the program are the nondissociated thermodynamic proper- 
ties of gases from Reference 2, which will provide the combustion character- 
istics with air of Methane (CH^), Hydrogen (H 2 ), and simulate those of JP fuel 
(CH 2 ). Computation of engine performance with any of these three fuels can 
be accomplished by input of the appropriate control parameter (IFUEL). The 
thermodynamic properties of air are also taken from Reference 2, .therefore, 
the temperature tables for air are also applicable to the fuels used, requir- 
ing the input of only one temperature table. 

All of the input to the program is by the namelist method and all of the 
parameters shown in the following list must be input for every design or off 
design point computed. The namelist method of input, however, provides the 
flexibility of only changing those parameters which are required to generate 
the desired engine performance. These parameters are normally Mach number, 
altitude, and the program controls which determine the number of afterburning 
and/or dry part power operating points. This is explained in the examples 
shown in Appendix B. 

As shown in Appendix B the user may select any combination of internal 
component design values, such as overall compressor pressure ratio, compressor 
efficiency, turbine inlet temperature, turbine polytropic efficiency, airflow 
rate, etc. The program will then dimensional ize the nondimensional compressor 
nap, determine the turbine flow function to be used as a matching parameter 
when computing off-design performance. It will also establish a nondimensional 
exhaust nozzle throat area to be used as a reference when computing all of the 
afterburning and part power data points requested. 

Output from the program may be obtained in several forms as listed below 
and explained later under the output options control coding. 

1. Long form printed output contains state point data at each engine 
location defined in Figure 6 as well as the overall engine performance. 

2. Long form printed output and card deck for use in LaRC airplane 
mission analysis computer program. 

3. Short form printed output - contains overall engine performance and 
key data for acoustic analysis only. 

4. Short form printed output and card deck for use in LaRC airplane 
,, missions analysis computer program. 
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•: 5. Printed output at a specified thrust level, providing performance 
data specifically for use in acoustic analysis. Data^ for acoustic 
analysis is also written on tape 8 which may be saved for use at a 
later time. 

The input required to operate the program, the variable name, control 
options, and iteration intervals used, along with explanation where necessary, 
are contained in the following list. The numbers included in the parameter 
in the following list refer to engine station locations as shown in Figure 6. 


Input Parameter Identification 


I. Single Value Inputs 

A. Design Installation and Flight Variable 

ALT - Geometric altitude (ft) - must equal 0 for the design case. 

CV - Exhaust nozzle internal gross thrust coefficient (decimal value) 

DTEMP - Desired increment from standard day temperature (°C) - used in 
subroutine AT65SP (atmosphere). 

1. If DTEMP is input with standard atmospheric pressures (PBTAB), 
the calculated performance will be for a simple hot day; 
i.e., pressure variation with altitude is the same as stand- 
ard day (pressure altitude) while the ambient temperature 

is for a standard day +DTEMP atmosphere. 

2. If DTEMP is input with the correct PBTAB for that specific 
atmosphere - the calculated performance will be for the 
correct hydrostatic equilibrium day. 

EMACH - Flight Mach number - must equal 0 for the design case. 

ENN - Compressor design rotational speed ratio - equals 1.0 for 
compressor-map used in the program. 

ETA2D - Compressor design adiabatic efficiency (decimal value). 

ETA4 - Combustion efficiency (decimal value). 

ETAS - Turbine polytropic efficiency (decimal value). 

ETAS - Afterburner combustion efficiency (decimal value). 

HPEXT - Power extraction from the turbine (horsepower, HP). 

PT2T1D - Compressor design total pressure ratio. 
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PT4T2D - Combustor design total pressure ratio, (decimal value). 

PT65AD - - Design afterburner total pressure ratio (constant, decimal 

value) - If PT65AD < 1.0, the data produced will be for an engine 
equipped with an afterburner. If PT65AD = 1.0, there is no 
afterburner pressure loss, therefore, a non-afterburning case 
will provide data the same as if the engine did not have an 
afterburner. 

PT8PT6 - Exhaust nozzle to afterburner total pressure ratio, (decimal 
value). 

TT4N - Design maximum turbine inlet temperature, (°R). 

TT6MAX - Maximum afterburner temperature, (®R). 

WID - Design sea level static air flow rate, (Ibm/sec) - must be the 
same as the parameter WIDE. 

WB - Bleed flow rate from compressor discharge, (Ibm/sec) - air bleed 
quantity is fixed for all operating conditions. 

WCOWl - Turbine cooling flow factor (decimal value), such that the cooling 
flow ratio WC/Wl = WCOWl (TT4-2160°R) . 

NOTE: Turbine cooling flow is returned to the cycle downstream 

of the turbine. The associated pressure loss is fixed 
in the program. 

WIDE - Design sea level static air flow rate, (Ibm/sec).- must be the 
same as the parameter WID. 

X - Factor for lower limit of the compressor pressure ratio - Equals 

0.1 for the compressor use in the program. 

XMNSLU - Flight Mach number at which the compressor rotational speed is 
held constant at 100%. 

Y - Compressor surge margin at design - Equals PT2T1S/PT2T1D, Equals 

1.1 for the compressor used in the program. 

B. Standard Constants 

Ad - Mechanical equivalent of heat (778.156 ft-lbf/Btu). 

G - Acceleration due to gravity (32.174 ft/sec^). 

GAMO - Ratio of specific heats for air (1.4)' - assumed constant for all 

operating conditions. 
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PSLS - Atmospheric pressure at sea level static’, standard day conditions 
(2116.2381 lbf/ft2). 

RO - Gas constant for air (53.33 ft-ibf/lbm~°R). 

TSLS - Atmospheric temperature at sea level static, standard day 
conditions (518.67°R). 

C. Iteration Intervals and Miscellaneous 


AVAL - Interval used in the partial power routine to define the desired 
accuracy of matching the nozzle throat area ratio with the 
scheduled nozzle throat ratio (decimal value, .005 = interval used). 

ICASE - Case number - program has no provision for automatic indexing. 

NLIM - Maximum number of interations in the Iteration Subroutine (ITT). 

T5VAL - Iterative constant for turbine outlet total temperature (decimal 
value, .001 = constant used). 

T7VAL - Iterative constant for nozzle throat static temperature (decimal 
value, .01 = constant used). 

T8VAL - Iterative constant for nozzle exit static temperature (decimal 
value, .001 = constant used). 

VAL - Interval used to determine nozzle throat total to static pressure 
ratio (choke) when the nozzle exhaust Mach number is supersonic 
(decimal value, .025 = interval used). 

II. Tabular Inputs 

NOTE: The numbers in parenthesis following the variable name are the 

number of values that must be contained in the variable array. 

A. Atmospheric Pressure 

PBTAB (8) - Atmospheric pressure for a standard hydrostatic equilibrium 
day from the U.S. Standard Atmosphere 1%2 - with DTEMP = 0°C 
and the PBTAB table input to the program, atmospheric pressures 
and temperatures as a function of altitude will be generated 
for a standard hydrostatic equilibrium day. Should DTEMP f 0 
be input, then to produce data for a non-standard hydrostatic 
equilibrium day a new PBTAB table derived from the U.S. Standard 
Atmosphere 1962 and corresponding to the input DTEMP will have 
to replace the current PBTAB table. Should the PBTAB table not 
be replaced the atmospheric data generated will be for a simple 
non-standard day, where the temperatures v;ill be non-standard 
and the pressures will be for a standard day. 
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B. Air Inlet Total Pressure Recovery Figure 4 

TABM (19) - Mach number. 

TABPTl (19) - Inlet pressure recovery. 

C. Compressor Map - Figure 1 

•TABETA (99) - Adiabatic efficiency ratioed .to the surge value 

of adiabatic efficiency at each compressor rotational 
speed (N). 

TABNST (9) - Corrected Compressor rotational speed (N/Zejl ) ratioed 

to design corrected speed. 


TABPT2 (9) - Surge line pressure ratio - surge pressure ratio at each 

speed (N) is ratioed to the surge pressure ratio at 
design speed (100% N). 

TABWl (99) - Cor recte d air flow ratio - corrected air flov/ 

(Wa/0ji76xi } is ratioed to the surge corrected air 
flow at each rotational speed (N). 

TABZ (11) - Interval of pressure ratio along a corrected speed line. 

D. Air Thermodynamic Tables - Reference 2 

TABH (78) - Enthaply of dry air (H, Btu/lbm). 

TABPHI (78) - Constant pressure entropy of air (^, Btu/lbm °R). 

TABTEM (78) - Temperature (°R, from 0 to 5000). 

E. Part Power Speed/Nozzle Throat Area Control - Figure 3 

TABA7 (12) - Nozzle throat area ratio. 

TABN (12) - Engine rotational speed ratio. 

F. Exhaust Nozzle Sound Suppressor - Figure 5 

TABDLDB (17) - Sound level suppression (dB) 

TABVJ (17) - Exhaust nozzle jet velocity (ft/sec)' 


III. Program Operating Controls 


A. Design and Operating Code - IDES 

IDES = 0 Used to calculate the design values of the turbine flow function 
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{W4/0.T4/^T4) nozzle area ratio (A 7 /A 70 ), which are 

then used in the matching process for that particular engine 
design at off-design operating' conditions. IDES must be 
input as 0 before normal operating cases are computed. 

IDES = 1 Following an engine design case (IDES = 0), IDES = 1 must be 
input to fix the design matching parameters which are used 
in subsequent normal operation. Then IDES = 1 must be main- 
tained for all normal operation with variable inputs of Mach 
number, power setting, altitude, etc., for an engine design 
determined previously with IDES = 0. 


Power Setting Codes - IPC, ZTHRUST and IDDBS 


IDDBS = 0 Provides engine performance without suppression. 


IDDBS = 1 Estimates the available jet suppression (dB) based on the 
computed exhaust jet velocity and the tables TABDLDB and 
TABVJ. It also reduces the computed gross thrust by 0.5% 
for each dB of estimated available suppression. However, 
when the exhaust jet velocity is equal to or less than 1000 
ft/sec no jet suppression is available and no thrust loss is 
assessed. 


IPC = 1 Provides power points from maximum afterburning to maximum 
dry power. 

IPC = 2 Provides a maximum dry power point only. 


IPC = 3 Provides power points from maximum afterburning to 

maximum dry power and on through minimum dry partial power. 

IPC = 4 Provides power points from maximum dry power through minimum 
dry partial power. 


ZTHRUST This control is used when data is required at a specific 
thrust value. 


ZTHRUST=0 This control is inoperative and the engine operating mode is 
determined by the IDDBS and IPC controls. 


ZTHRUST>0 When this control is specified the IPC control must be set 
to 1, 3, or 4, since the ZTHRUST control causes tables of 
data from maximum through minimum partial power to be gen- 
erated. From these tables at the specified ZTHRUST are 
interpolated the specific engine parameters required for engine 
acoustic analysis. Use of this control overrides the normal 
output control and provides a printed output as shown in 
Example 6 of Appendix B. It also causes the program to write 
the parameters required for acoustics analysis on tape 8 
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which can be saved. The parameters written on tape 8 and the 
format in which they are written are given below: 

Parameters: CAPAS, WG, V8, FN, EMACH, T8 

(For parameter definitions, see output parameter 
identification) 

Format: (F10.4, F10.3, F10.2, FlO.l, F10.4, F10.2) 

C. Fuel Code - IFUEL 

The fuel code will cause the program to select the thermodynamic 
properties of the desired fuel with' which to compute the thermo- 
dynamic properties of combustion and subsequently, engine performance. 

IFUEL = 1 Selects JP4 (CH2) fuel. 

IFUEL = 2 Selects methane (CH4} fuel. 

IFUEL = 3 Selects hydrogen (LH2) fuel. 

D. Variable Geometry Turbine - A4A4DN 

The variable geometry turbine can only be employed for partial power 
engine performance calculations, therefore, use of this control can 
only be accomplished when IPC = 3 or 4 which causes entry into the 
part power routine. 

A4A4DN s 1.0 will produce part power performance for standard 
fixed area turbine engines. 

A4A4DN < 1.0 defines the lower limit of turbine area variation from 
the design area. The routine will divide the difference 
between 1.0 and the input value of A4A4DN into 0.05 incre- 
ments, from which the program will adjust the turbine inlet 
temperature to provide a flow match while maintaining a 
constant compressor pressure ratio and compressor corrected 
inlet flow. When the incremented value comes within 0.025 
of the input value of A4A4DN, the routine will iterate to 
define performance conditions at A4A4DN. 

A7A7DN > 1.0 with A4A4DN < 1.0. will provide partial power data at 
values lower than that permitted by the A4A4DN control 
alone. With A4A4DN i 1.0 will provide partial power data 
by standard routine for a fixed geometry turbine. 

A7A7DN < The partial power data will be determined by the A4A4DN 
control only. 
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E. Afterburner 


IPS - Number of afterburner points from maximum afterburning through 
minimum afterburning. 

IPS = 1 with IPC = 1 or 3 produces maximum dry power only. 

IPS = 2 with IPC = 1 or 3 produces no data and causes the 
program to dump. 

IPS = 3 with IPC = 1 produces maximum afterburning power, 
minimum afterburning and maximum dry power. With IPC = 3, in 
addition to the above, the program will produce data from 
maximum dry power through minimum partial power, 

IPS > 3 with IPC =1 or 3 will produce data similar to. the 
data when IPS = '3 except that there will be additional data 
points between the maximum and minimum afterburning data points. 

F. Nozzle Codes 

NFINAL - Number of partial power cases desired - must be equal to or less 
than the number of values contained in TABA7 or TABN. 

NSTART - Counter to determine the interval in the TABA7/TABN tables for 
the part power calculation: 

NSTART = 1 Read the first table value for the first calculation 

NSTART = 2 Read the second table value for the first calcu- 
lation. 

NUMBER - Sequence of TABA7/TABN lookup. 

NUMBER = 1 Use each value. 

NUMBER = 2 Use every other value. 

P80P8 - Exhaust nozzle exit static pressure ratio - equal to 1.0 

provides performance for a fully expanded convergent-divergent 
nozzle. 

NOTE: At the present- time this is the only option available. 

G. Output Code - IPUNCH 

IPUNCH = 0 Long form printed output only (Examples 1-4, Appendix B) 

IPUNCH = 1 Long form printed output plus punched cards in the missions 
analysis format. 
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Card output' parameters and mission format are:* 

Parameters; EMACH, ALT, POWER, FG, FD, WF, V8, CAPAS, WG 

(For parameter definitions, see output parameter 
- identification.) 

■ Mission Format: (F5.2, FlO.l, F5.0, 4F10.1 , F10.3, FlO.l ) 

IPUNCH = 2 Short form printed output only (Example 5, Appendix B) 

IPUNCH = 3 Short form printed output plus punched cards in the missions 
analysis format. 


(See also the explanation of the ZTHRUST for other modes of 
output. ) 


The output parameters along with their definitions are presented in the 
following list: 


A4/A4D 

A7/A7D 

A8 

A8/A7 

ALT 

CAPAS 

CASE = I CASE 
CV 

DFGS 


Output Parameter Identification 

- Turbine flow area ratio - actual to design. 

- Nozzle throat area ratio - actual to design. 

- Acoustic velocity of the exhaust gas at the nozzle exit, 
(ft/sec) 

- Nozzle area ratio - exit to throat. 

- Geometric altitude (ft) (input). 

- Nozzle exit area (ft ). 

- Case number - program has no provision for indexing 
(input). 

- Nozzle internal gross thrust coefficient (input). 

- Gross thrust loss due to exhaust jet noise suppression 
Obf). 


DH2-1 

DLDB 

DTEMP 

ETA2 


“ The actual enthalpy .change across the compressor minus 1. 
(Btu/lbm) (design case only). 

- Available exhaust jet noise suppression (dB). 

- Temperature variation from standard day temperature (°C) 
(input). 

- Compressor adiabatic efficiency (decimal value). 
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ETA5 

ETAC 

ETAO 

ETAP 

F4 

F6 

FD 

FG 

F61 

FN 

FNl 

’fn/wi 

FN/WID 

HPEXT 

HTl 

HT2I 

IDDBS 

M8 

MACH = 
N 

N/SRT 

ONSRT 


- Turbine polytropic efficiency (decimal value-input) 
(design case only). 

- Cycle efficiency (decimal value). 

- Overall propulsion efficiency (decimal value) 

- Propulsive efficiency (decimal value) 

- Primary combustor fuel air ratio 

- Afterburner fuel air ratio (if nonafterburning then 
F6 = F5A = fuel air ratio at afterburner inlet). 

- Ram drag (Ibf). 

. - Gross thrust (Ibf). 

- Gross thrust corrected for exhaust jet suppression loss 
(Ibf). 

- Net engine thrust (Ibf), 

■ ' - Net engine thrust corrected for exhaust jet suppression 

loss (Ibf). 

- Specific net thrust (Ibf/lbm/sec). 

- Specific net thrust based on sea level static design 
airflow (Ibf/lbm/sec). 

- Power extracted (horsepower) (input). 

- Enthalpy of the air entering the compressor (Btu/lbm) 
(design case only). 

- Ideal enthalpy of the air at the compressor discharge. 
(Btu/lbm) (design case only). 

- Exhaust jet suppressor control parameter, (input) 

- Nozzle exit Mach number. 

EMACH - Flight Mach number, (input) 

- Compressor rotational speed (%RPM). 

- Corrected compressor rotational speed - N/ /0j] (%RPM) 

- Corrected compressor rotational speed - N/Zejl (%RPM) 
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P3 = PS 

PC = IPC 

POWER = TT4 

PTIPTO 

PTl/PTO 

PT1PT7 

PT2T1D 

PT2/PT1 

PT4/PT2 

PT5/PT4 , 

PT5T5A 

PT6/PT5A 

PT6T5A 

PT7P7 

PT8PT6 

PT8/P8 

SFC 

SECT 

T8 

TTO 

TTl 

TT2 


(design case only). 

- Afterburner power point index - varies from 1 through IPS. 

For a nonafterburning case, it is equal to I = indeterminant. 

- Engine power setting code, (input) 

- Turbine inlet temperature (°R). 

- Inlet pressure recovery (design case' only). 

“ Inlet pressure recovery. 

- Engine pressure ratio. 

- Design compressor pressure ratio (input) (design case only). 

“ Compressor pressure ratio. 

- Combustor pressure ratio. 

- Turbine expansion ratio. 

- Pressure. ratio - total pressure before turbine cooling air 
is returned to the cycle to total pressure after the 
turbine cooling air is returned to the cycle - PT5/PT5A. 

- Afterburner total pressure ratio, 

- Afterburner total pressure ratio - PT6/PT5A. 

- Total to static pressure ratio at exhaust nozzle throat. 

- Exhaust nozzle exit to afterburner total pressure ratio. 

Total to static pressure ratio at the nozzle exit. 

- Specific fuel consumption (Ibm/hr/lbf). 

- Specific fuel consumption corrected for exhaust jet 
suppression loss (Ibm/hr/lbf). 

- Nozzle exit total temperature (°R). 

- Free stream total temperature (°R). 

- Total temperature at the compressor inlet (°R). 

- Compressor discharge total temperature (°R). 
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TT4 

TT5 

TT5A 

TT6 . 

VO 

V8 

W1 

WIDCOR 

WIK 

WB 

WBOWl 

WC 

WCOWl 

WF 

WG 

Z 

ZTHRUST 


- Combustor exit (turbine inlet) total temperature (°R). ■ 

- Turbine discharge total temperature (°R). 

- Total temperature at the afterburner inlet (°R). 
(Temperature after the turbine cooling air has been 

returned to the cycle.) 

- Afterburner exit total temperature (®R). 

“ Flight velocity (ft/sec). 

“ Nozzle exit velocity (ft/sec). 

- Actual compressor inlet air flow (Ibm/sec). 

- Corrected compressor air flow ratioed to the surge 
value at a given speed ratio (design case only). 

- Corrected compressor inlet airflow (W/efl/STl - Ibm/sec). 

- Service air bleed (Ibm/sec) (input). 

- Ratio of service air bleed to actual compressor air flow 
(WB/Wl). 

- Turbine cooling air flow (Ibm/sec). 

- Turbine cooling bleed factor [Turbine cooling airbleed 
ratio = WC/Wl = WCOWl (TT4~2160°R)] (input). 

- Total engine fuel flow = engine fuel flow plus afterburner 
fuel flow (Ibm/hr). 

“ Engine exhaust gas flow (Ibm/sec). 

- Compressor pressure ratio function. 

- Engine power control (input). 
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Surge line pressure ratio ratio, 

(CPR) surge surge at (N^orr^d 



Rotational speed ratio, /(Ngorr^d 


(a) Surge Line Pressure Ratio Ratio as a Function of Rotationol 
Speed Ratio 

Figure I - Nondimensional Compressor Map 
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Adiabatic efficiency ratio, 
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Corrected flow, Wcorr.l / < Wcorr, |)surge > 


(c) Efficiency as a Function of Corrected Flow 


Figure I. - Concluded 
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Figure 2- Variation of Compressor Corrected Inlet Flow w 
of Design Engines, 
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Figure 3- Schedule of Exhaust-Nozzle- 
Throat Area Ratio used for Partial 
Power Control. 



Moch number, M 

Figure 4- Schedule of Air- Inlet Total Pressure Recovery 
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Jet suppressor effectiveness— AdB 



Exhaust' Jet velocity- ft/sec 


Figure 5- Predicted Jet Noise Suppressor Effectiveness 
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Station 

0 Free stream 

1 Compressor Inlet 

2 Compressor outlet 

3 Combustor Inlet 

4 Combustor outlet {Turbine nozzle Inlet) 

5 Turbine outlet 

5 A Turbine cooling air returned to cycle (Afterburner Inlet) 

6 Afterburner outlet 

7. Ex ha list -nozzle throot 

8 Exhaust-nozzle exit 


Figure 6 - Engine Station Identification 



APPENDIX A 
ENGINEL 

PROGRAM LISTING 
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PROGRAM ENGINEL( INPUT»0UTPUT*TAPE5=INPUT«TAPE6=0UTPUT«PUNCH. TAPES) A 1 


COMMON /lEP/ IE A 2 

1E=0 A 3 

call main a 4 

END A 5— 


> 

I 

ro 


SUBROUTINE MAIN 81 

C MAIN PROGRAM B 2 

DIMENSION ALTSAV(15)o EMSAV(15)» NOSAV ( 1 5 > ♦ , FNS A V ( 1 5 )» W2SAV(15)« B 3 
ISFCSAVdS) B 4 

dimension TABTEM(78>* TA8H(78)» TABPHI(78)« TABPSH(78)» TABM(19)» B 5 
1TABPTI(19)« TABNST(9)« TABPT2{09)» TABZ ( II ) ♦ TABW1(99)» TABETA<99) B 6 
2» ANS(5). SFCC(25)» TABPSP ( 78 ) B ''7 

dimension TABPH5(78)» TABH4(78) B 8 

dimension TABN(I2)» TABA7(12)» RAT(2)» T4(3) B 9 

DIMENSION PBTAB(8)i HDTI20)* DTTABI20) B 10 

DIMENSION 2FN(50), YFN(50»10)« ZEE(IO) B 11 

dimension TABVj(17)» TABDLDB(17) B 12 

DIMENSION HJ<3)» HSJ<3)» TFJ(3)» FSJ(3)« AMP J < 3 ) B 13 

DIMENSION TABH4J(78.3) » TABPH5 J ( 78 ♦ 3 ) ♦ T ABPSH J ( 78 . 3 ) , TABPSPJ(78«3 B' 14 
1 ) B 15 

dimension ZZ(936) B 16 

EQUIVALENCE ( T ABH4 J ( 1 > » ZZ ( 1 ) ) » ( TABPH5 J ( 1 ) * ZZ ( 235 ) ) B 17 

EQUIVALENCE ( T ABPSH J ( 1 ) * ZZ ( 469 ) ) » ( TABPSP J ( 1 )» ZZ ( 703 ) ) B 18 

C B 19 

COMMON /ZPARAM/ ZFN « YFN ♦ NZCNT . ZTHRUST 8 20 

I NPUT B 2 1 

B 22 

COMMON /BKl/ ALT«EMACHtPT2Tl D»TT4N» WCOWl » W8* A7A7DN 0PT65AD *A4A4DN»P B 23 
180P8»W1D«H»ETA2D.ENN.X»Y.TF»PT8PT6»CV»TT6MAX»PT6T5A«ETA6 » IDES.VALo B 24 
2GAMO»RO»AJ.TSLS»PSLS,G»TABTEM,TABH»TABPHI dAePSH, TASPSP oTABM oTABNS B 25 
3T*TABPT2«TABZ»TABW1 *TABETA*.TABPT1 . I DDBS«TABVJ » TABDLDB, IFUEL B 26 

B 27 

COMPUTED values B 28 


B 29 

COMMON /BK2/ WBOw 1 » AOAS * ABCOR ♦ A80A7 . A7A7D ♦ ANS ♦ APO t A4A4D » A7 « A8 ♦ AO » A B 30 
180W1 .ABEM7«T0»T8,T7*T7I * T8 I ,TT0«TT2»TT1 »TT2I »TT2TTi »TTS,TT5A «TT7TT B 31 
21 ,TT69TT7TT6 iTT 4,P0»P80P0*PT7*HF»HBM«H8,H8I » H8 I M « H7M » H7 « H7M I tHAStH B 32 
3AS5<i H0» HT09HT29 HT2I « HT4M . HT4 * HT5M ♦ HT5 » HT505 ♦ HT5AM 9 HT5H 9 HT7M » HT7M I » B 33 
4HT6M9HT8M9HT5A05,PHI09PHIT0,PHIT2I »PHIT1 » PH I T4M » PH I T4 ♦ PH I T5A , PHT5A B 34 


OfllGjQ^AL PAGE IS 
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1 

CO 


PH !T8M, PH 10 *PH 12 I »PHlT6Mf PHI 7M i PH I T6 » PH I T7M ♦ PH0F6 « PH0H6 « PT B 35 
60P0,PT1PT0«PT2T1S*PT2PT1 ,PT4PT2 » PT5PT4 , PT5T5A , PT8P8 , PT 1 PT7*PT7P7 «P B 36 
7T6PT5»PT6T5D»PT5AM,PT55AD,PSIH4 ,PSP14»PSIH5»PSP15A,PSIH5A,PSPI8 tPS B 37 
8IH8«PSIH7»PSPI6»PSIPHI « PS I H6 » PS I H8 I » PSP I 7i EN ♦ ETA2 . ET5A , ETETS » ET A2S B 38 
9*ET2D«EM7»EM8« W1 COR o W1 DCOR ♦ W4C0R ♦ W4DC0R » W7COR ♦ W7DCOR * W4TOT4 « W4TEST B 39 
S«W1 «W1K»DH2M1 <STT0«DELTO9DELTOD*GAMA8«GAMA7»VO»V7»V3»GAM7Pi «GAM7Ml B 40 
$*0NSRT»SFC»SFCC»0F4»0F6»F0F6*F4«F5A»FAB,F6»FN0WI » FN * F6A » F0F4 » F05A » B 41 
SW8.RT5 9RTA»RT6»RT5A»CP7»CP8 9CAPA8»CAPA7»0WB»WC8* I'COMB* IPP» IS. ICOUN B 42 
$T» IKONT » JCONT. IC.KTEST .THETI .THETO.OWl .DVU2.DDHP.Z1 . Z2 » Z3« Z .ETAC .ET B 43 

SAO. FNOW1D.HT6.PHIH6.T5A.T4.RAT. WCOWIC.PBTAB.HDT.DTTAB. IPS.OTEMP.ET B 44 

SAP. alts A V.EMSAV. NOS AV.FNSAV.W2SAV.SFCSAV.no. 1N0Z2.TT2T1D.HT5A.PHT1 B 45 
5,TT4TT2.PHIT5.PSPI5.PHIT5!^»«HT1 .PT2T1L.PS, T ABA7 . T8 V AL » T7VAL ♦ NL I M . NU B 46 
SMBER.nSTART.NF INAL. AVAL .XMNSLU. T ABN* ETA4. ETA5. TABPH5.TABH4 .HPEXT . I B 47 


SCASE. IPC.PT4T2d.HS.T5VAL.WC B 48 

common /IEP/ IE B 49 

COMMON /NEW/ FD.FG.WF.DLDB.DFGS.'^GI »FN1 .SFCl B 50 

COMMON /RRR/ AMP28 . FS28MP . AMP B 51 

COMMON /PPP/ I PUNCH - B 52 

COMMON /WVUW/ wide B 53 

namelist /ENPUT/ ALT.EMACH.PT2T1D.TT4N.WC0W1 .WB.A7A7DN.PT65AD.A4A4 8 54 

lDN.PaOP8.WlD* ETA2D.ENN.X.Y. PT8PT6 . CV . T T 6MAX . ET A6 » I DES . VAL . GAM B 55 
20.R0.AJ.TSLS.PSLS.G.TABTEM.TABH.TABPHI ♦ T ABM ♦ TABNST « T B 56 

3A8PT?.TA8Z.TABW1 .TABETA.TABPTl .ETA4.ETA5. HPEXT. ICASE. B 57 

4IPC.PT4T2D. T5VAL.T8VAl.T7VAL.NLIM. number. NSTART.NFINAL. aval. XMN B 58 
5SLU.TABN.TABA7.DTEMP.PBTAB. IPS. WIDE* I PUNCH . ZTHRUST . 1 DOBS . TA B 59 

6BVJ.TABDLDB. IFUEL B 6° 

data HJ/20000. .23919,. 62000./ B 61 

DATA HS J/1600, .231 9, ♦ 3 0400o/ B 62 

data TFJ/536, .536. .536../ B 63 

DATA FSJ/.06763 . o 05801 . o0291 6/ B 64 

DATA AMPJ/28,907.27o74 1 . 24 .648/ B 65 


data (ZZ( I ). 1=1 .781/0.00000. 180.50000. 255,02000. 267, 44000. 279, 9400 B 66 
10 .292.52000. 305, 1 9000.317.93000 . 330 . 76000 , 343. 67000 . 356 . 68000 . 369, B 67 
276000.382.94000.396,21000 . 409 . 570DO . 423 . 04000 , 436. 6d000 , 450 , 28000 . B 68 
3464,06000.477.94 000.491.92000,506.08000.520.28000,534,58000.549,05 B 69'; 
4000,563.55000.578 » 1 5000.592,85000,607.70000.622,58000.637,56000.65 B 70 ' 
52. 64000. 667 .79000 ,683. 02000. 698 ,33000. 71 3.73000.729,20000.744 .7300 8 71 

60.760,34000.776.01000.791 . 76000 , 807 .67000 , 823 , 43000 , 839, 35000 .855, B 72 
734000,871 ,38000.887,48000,903,60000.919,82000,936.06000,952,36000. B 73 
8968.70000.985.08000, 1 001 ,51000. 101 7.98000, 1 034 .49000. 1051.04000.1 0 B 74 
967.63000, 1084,26000 , 1 100,92000. 1 1 17,62000, t 134.35000,1 151.1 1000 .11 B 75 
S67.91000. 1 184,74000, 1201 ,59000, 1218.47000. 1235.39000.1252.33000 .12 B 76 
$69,30000, 1286.29000, 1303,31000, 1320,35000,1337,42000,1354,51 000 ,13 B 77 
$71.62000,1543.89000,1717,92000/ 3 78 


O O 

sg 



DATA {Z2( I ) ♦ 1=79, 1S6 )/0 o 00000 ,239, 51 000 , 307 o 33000 « 320 « 30000 » 333 o35 
1 000,346*47000 ,359,66000,372 .92000 ,3S6« 25000 , 399,66000,413*16000,42 
26.73000,440.38000,454* 13000,467.97000,481.91 000 , 495. 94000 , 51 0 * 1 000 
30,524*34000,538.70000,553. 15000,567.81 000,582*48000,597.27000,612* 
425000,627.24000,642.33000 , 657.53000,672. 90000,688.29000, 703*78000* 
5719,39000,735.06000 ,750, 82000, 766 *67000 ,782. 60000, 798*6 1000, 81 4. 70 
6000,830*86000 , 847,10000,863.42000,879.79000,896.24000,912.74000,92 
79.32000,945*96000,962*65000,979,40000,996*21000,1013.07000,1029,99 
8000, 1046,96000, 1063*97000, 1081,04000*1098.15000,1 1 15.31000*1 132,51 
9000,1 149.55000 , 1 1 67*05000, 1 184,37000, 1201 .73000,1219.13000*1236*57 
$000,1254.05000 ,1271 .56000, 1289. 1 0000, 1 306.68000,1 324.28000,1341 .93 
$000,1359,59000 , 1 377.29000, 1 395,01000,1412.77000,1430.55000,1 448,35 
$000, 1466,18000, 1645.77000, 1827.33000/ 

DATA (Z2{ I ) , 1=1 57 ,234 1/0.00000,381 ,94000,470*1 4000,484.84000,499.5 
180 00,51 4.37000,529, 1 9000 , 544 .06000 , 558. 97000 * 573, 93000 , 588*95000*6 
204,02000,619.17000,634.40000,649.71000,665,1 1000,680.61000,696.240 
300,71 1 ,95000,727.77000,743.69000,759.92000,776*09000.792.36000*808 
4. 95000, 825. 5000, 842. 08000, 858.82000,875.81 000 ,892.77000,909.85000 
5,927,12000,944,42000 ,961 ,83000,979, 33000,996.95000 , 1 0 1 4 . 66000 , 1 032 
6*46000, 1050,35000,1068,33000 , 1 086,42000, i 104.58000 , 1 122.83000*1 141 
7,14000, 1 1 59.56000 , 1 1 78.05000 , 1 1 96.62000 , 1 2 1 5 . 26000 , 1 233 . 97000 * 1 252 
8.75000, 1271.60000 ,1290.52000,1309,50000, 1328.56000 » 1 347 , 67000 , 1 366 
9.84000 , i 386.06000 , 1405.34000 ,1424 .69000, 1444.08000 , 1463,52000*1 483 
$.01000, 1502.55000,1522.15000,1541 ,79000, 1561 ,47000,1581.19000,1600 
$.95000, 1620.77000 , 1 640.62000,1660,50000, 1680.42000 » 1 700. 39000* 1 720 
$.39000, 1740.42000 , 1 760.49000 .1962,86000,2167.87000/ 

DATA (ZZ( I ), 1=235,3121/0.00000,1 .29510, 1 ,49442,1 ,52764,1 .5571 1,1.5 
1 8363, 1.60778, 1 .62995,1 .65050,1 .66965* 1-.68766* 1 .70455*1 .72053*1.735 
271,1 .75017, 1 *76398, 1 .77723, 1 .78991 , 1 *80216* 1 .81398*1 .82540* 1 *83646 
3* 1 .84722*1 .85758, 1 *86770,1 *87754 » 1 .8871 1 , 1 .89645, 1 *90559* 1 *91448.1 
4.92316* 1 *93166*1 .93996* 1 .94809. 1 .95604* 1 .96383*1 *97147* 1 *97896* 1 *9 
58631 , 1 *99351 . 2 . 00059 *2 * 00754 , 2 . 0 I 43.6 * 2 , 02l 07 * 2 • 02766 * 2 * 034 1 4 * 2* 040 
652,2.04678,2*05896,2.05903*2.06501 ,2.07090,2.07670*2*08241 *2.088o4 
7,2*09359*2.09907,2* 1 0446,2.1 0978,2* 1 1503,2.12021,2*12531 ,2*13036,2 
'8.13533,2. 14024,2, 14509,2.1 4988,2. 15461,2. 15929,2. 16391 ,2* 16847.2. 1 
97298,2. 17743,2.18184,2.18619,2.1 9050,2.23107,2.26774/ 

DATA (ZZ( n , 1=313,390 1/0,00000,1 .34074, 1 .54894 , 1 *58364* 1 .61439. 1 .6 
14204,1 .66719,1 .69026,1 .71162,1 .73150*1.75020* 1 .76771 « 1 *78427*1*799 
299*1 .81496* 1 *82926. 1 .84296* 1 .85608. 1 *86905,1 *88097,1 *89279,1 *90423 
3,1 *91525,1 *92506,1 *93653,1 .94670, 1 .95660,1 .96625, 1 .97570,1 *98489, 1 
4.99387,2*00265,2.01 124,2.01 965,2*02788,2.03595,2.04385*2*05161 *2.0 
55922,2.06668*2.07401 .2.08121 , 2* 08829 *2 .09524 ♦ 2 * 1 0207 *2 . 1 0880 *2. 1 15 
641,2.12191,2.1 2832,2.1 3462,2*14083,2. 14694,2. 1 5296,2.15890*2*16475 
7,2*1 7052,2.17621 ,2.18181 ,2 . 1 8734 ,2 . 1 9280 , 2 . 1 98 1 9 , 2 *20350 , 2*20875,2 
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a o 21 392 -» 2 o21 903 o2o 22403 0 2 O 22907^2 0 53399, 2 <.23aS6«2o 24367 o2» 248435 2 0 2 B 123 
9531 1 »2.25776»2o26234 *2o2660a»2«27137,2<.3136752»35392/ B 1124 

DATA (ZZ ( n 0 I =391 ,468 )/OoOOOOO, 1 o48289« 1 o71 863 9 1 o75792» 1 o79266« 1 o8 B 125 
12383, 1 0 05210,1 ,87796,1 ,90134 , 1 ,92401 , 1 ,94483, 1 ,96429, 1 <>98265, 2,000 B 126 
204,2o01 660, 2 ,0 3240, 2 ,04753, 2 ,06204, 2 • 07601 ,2,08948,2,10249,2,11 SI 4 B 127 
3,2,12731 ,2,13918,2,15075,2,16194,2, 17285,2,18348,2,19390,2^20403,2 B 128 
4,21393,2,22360,2,23303,2,24237,2.25146,2,26038,2.26912,2,27771 ,2,2 B 129 
586 13, 2. 29439, 2, 30252, 2, 3 1050, 2, 3 1836, 2, 32607, 2, 33366, 2 ,341 13,2,348 B 130 
649,2,35560,2,36286,2 ,36987,2,37679, 2,38361 ,2,39033,2,39695,2,40349 B 131 
7,2.40993,2.41629, 2, 42256 ,‘2 ,42875, 2 9 43486, 2 9 44089, 2 ,44684 ,2,45272,2 B 1 32 
8.45852,2,46425,2.46992,2,47551 ,2,48104,2,48651 ,2,49191 ,2,49725,2,5 B 133 
90253,2.50775,2,51291 ,2,51801 ,2.52306,2,57072,2,61392/ B 134 

DATA (ZZ( I ), 1=469, 546)/0. 00000, 934,30000, 982, 30000, 990,30000, 999,1 B 135 
1 0000,1008,70000 ,1019 .30000, 1 030.60000, 1042,80000 i 1 055.800C0, 1070.3 B 1 36 
20000,1084.70000, 1099.80000, 1 1 15.90000, 1132,10000,1 148.80000, 1165.9 B 1 37 
30000,11 83.90000 » 1 202 ,00000, 1 220,40000, 1239,30000, 1 259, 60,000 , 1 279. 3 B I 38 

40000.1299.40000 . 1 320 .90000 . 1 341 .70000. 1363.00000. 1 384 . 60000 . 1 407 .3 .B 1 39 
50000,1429.60000, 1452.300CC , 1475.50000, 1 499.00000, 1 522,70000, 1546.8 B 140 
60000, 1571 *30000, 1596,00000, 1621 ,10000, 1646,40000, 1672,10000, 1698,0 B 141 

70000. 1724.20000. 1750.80000. 1777.40000. 1804.50000. 1831 .00000. 1859.4 8 142 

BOOOO, 1886.90000 , 1915.10000, 1943,40000, 1971 ,80000,2000,50000.2029,4 B 143 
90000,2058.50000 ,2037.70000 ,211 7,20000,2146.80000,2176.70000,2206.7 B 1 44 
£0000,2236.80000 ,2267.20000, 2297 .‘60000 ,2328.20000 ,2358.90000,2339,8 B 145 
$0000,2420.30000,2451 ,90000,2483, 10000,2514,50000,2545,90000,2577,4 B 146 
$0000,2609,00000 ,2640.70000, 2672.50000,2704,40000,2736,40000,3059,5 B 147 
$0000,3387,20000/ B 148 

DATA (ZZ( I) , 1=547, 6241/0,00000, 1972,50000,2088,90000,2108.30000,21 B 149 

128.40000.2149.30000.2171 .10000.2193.60000.2216.90000.2241 .00000. 22 B 150 
266,60000,2292,1 0000,2318,30000,2345,60000,2372,90000,2400 ,80000.24 B 151 

329.10000.2458.60000.2488.00000. 2518.00000.2548.40000.2581 0 1 0000.26 0 152 

4 12,60000,2644 .50000 ,2678,70000,27! 1 , 70000 , 2745 , 20000 , 2779 , 30000 .28 B 153 

514 .90000. 2850 .00000. 2885. 50 000.2922 .00000.2958. 60000. 2995 .70000 .30 B 154 
633.30000,3071 ,40000,31 09,90000,3148,90000,3188,30000,3228.20000 .32 B 155 
768,60000 ,3309.30000 ,3350.50000,3392 o 00000 ,3434 ,00000,3476,4 0000 ,35 B I 56 
8 19. 20000, 3562, 50000, 3605 ,80000, 3649, 70000,3693,90000,3738,50000 ,37 B 1 57 

983 .30000. 3828.60000. 3874.1 0000.3920.00000. 3966.10000.4012.50000.40 B 153 

$59, 20000, 4 106, 10000,41 53. 30000, 4200 ,80000, 4248 .40000, 4296, 30000, 43 B 159 
£44.50000,4392.30000 , 4441 .30000,4490,1 0000 . 4539 , OOOOO « 4580 , 10000.46 B 160 
$37,30000,4686.70000,4736.30000,4786, 1 0000,4835 ,90000 , 4836 ,00000 , 53 B 161 
£92.70000,5900.70000/ B 162 

DATA (ZZ( I ), 1=625, 702)/0, OOOOO, 9183, 30000, 9788, 10000, 9886, 60000,99 B 163 
185.40000, 10084,70000,10184,60000,10284,90000. 10385,70000, 10486,800 B 164 

200 . 1 0590.30000 . 10692.70000.1 0795.90000. 1 0901 .30000. 11006.00000. 111 B 1 65 

3! 1 ,60000, 1 1218. 10000, 1 1327,40000, 1 1435.90000, 1 1545,70000, 11656,300 3 165 


ORIGINAL PAGE IS 
OP POOR QUALITY 



A-6 


400.1 1775»00000, 1 I383o30000, 12003.. iOOOO, 121 1 9 » 50000 ♦ 1 2242o 90000« 12.3 
561 » 90000, 12482 .,80000, 1 2608 » 30000, 12731 o 80000, 1 2856 o 90000, I2935oao0 

600. 1 31 1 4« 1 0000, 1 3244 0 00000 , 1 3375 o 4 0000, 1 3508 o 30000 , 1 3642 o. 80000 » 1 37 

778.00000. 13916.20000. 14055.00000. 14195.50000. 14337.20000. 14480.600 
800 ,14624,70000 , 14771 ,20000, 14918,50000, 15067,40000,151 1 7,50000, 153 
968,80000, 15521 ,40000, 15675,30000, 15830,40000, 15986,60000, 16144,200 
•eOO, 16302,30000 , 16462,50000, 16623,10000, 16785,00000, 16947,80000, 171 
£11 ,60000, 17276.30000, 17441 ,90000, 17608,30000, 17775,70000, 17943,900 
$00 ,18113,00000, 13282 ,60000 , 1 8453 ,20000, 18624,50000 ,18796, 40000,1 89 
$69 9 00000, 19142,20000 ,19316,20000 , 1 9490,80000, 19665,80000, 19841,600 
£00,21626,50000 , 23452,30000/ 

DATA (Z2{ I ) , 1=703,780 )/0, 00000,-, 1 9680 ,06420 0421 0 02l 00 O 
1 0060 , ,0 1 930, ,03900 , , 05850, ,07800 , ,09730, , 1 1 650 , , 1 3550, ,15370 ,,.1713 
20, .18830, ,20520, ,22150, ,23730, ,25290, ,26850 < ,28420 , ,29900., , 3 1 370 » , 

332830. . 34260. .35640. .37020^ ,38380, , 39720 , ,4 1 04Q , , A 2320 , , 4 3590 , , 448 
450, ,46100, «/( 7340, ,48560 , ,49770 ,, 50970 ♦ ,521 40 , ,53310, ,54460, ,55600, 

5.56730. . 57850 . . 58950 .. 60040 ..61110. . 62 1 80 . . 63240 . .64280. .65320. .66 

6340. . 67350. .68350 . .69340. .70320. .71290. .72250. .73200. . 74 1 4 0 . . 75070 

7. . 75990 . .76900 . .7780c , ,78690 , ,79570 , ,80450, , 8 1 3 1 0 , » 82 1 70 , ,83020 , ®8 
33850, ,84680, ,8550''', ,36320, ,87120, ,94730, 1 o0l530> 

DATA (7Z( I ) , 1=781 , 853 ) /O , 00000 , ,60580, ,92020, ,97260,1 ,020 30,1 ,0 645 
1 0, r, 10590 , 1 , 14500, 1 , 18240, 1 ,21830 , 1 ,25370 , 1 , 28660® 1 ,31820 , 1 ,34990* 
£ 1 ,37950 , 1 ,40800 , 1 ,4 3590, 1,46230,1 .48370, 1 ,51410,1 ,53910 , i ,56430 , 1 , 
358330, 1,61130,1 ,63460, 1 ,65720,1 .67900 , 1 ,70070, 1 ,72210,1 ,74300,1 ,76 
4370, 1 ,78350, 1 ,80350, 1 ,82330, I ,84290, 1 ,86210, 1,88120*1 ,90000* 1,9185 
50,1 .93690, 1 ,95500 , 1 ,97290 , 1 ,99060,2,00810,2,02540 ,2,04260,2,05950* 
62,07620,2.09280*2*1 0920,2, 125^0,2, 14150,2. 1574 0,2, 17310,2, 18870 *2, 
7204 10, 2, 2 1930, 2.23440, 2. 24940,2,264 10, 2. 27860, 2 ,29330 *2, 30760 *2 ,32 
8180,2,33580,2,34980 « 2 , 36350 , 2 ,37720 , 2 ,39070 ,2 , 40400 ,2 ,4 1 720 * 2 ,4303 
90,2,44330,2,456! 0,2,46890 ,2,48150,2,60080,2,70950/ 

DATA <2Z ( ! ) , I =85« ,936 1/0,00000, 6, 1 8860,7,76960,8, 0331 0 « 8 , 26620 * 8° 4 
1 7560, 8. 66590, 8,84'' I 0,9, 00 140,9, 151 50 *9,2951 0*9,42740,9,552 10,9,672 
.260,9,78500,9,89400,9,99340, 10,09840* 10,19540, 10,28860, 10,37950* 10, 
34 7070, ! 0.56350, 10,63970, 10,72330 , 10,80340, 10,831 1 0 * 1.0,95770, 1 t ,033 
430,11 ,10700, 1 1 .17970,1 1 ,24910,1 1 ,31940, 11 ,38880,1 1 ,45700* I 1 ,52430* 
51 i ,59070, 1 1 ,65630,.! I ,72090, 1 1 ,73470,1 1 ,84700,1 1 ,91020,1 1,97180*12, 
603250* 12,09290,12,15250, 12,21140, 1 2 , 26890 , 1 2 , 32740 * 1 2 ,38440 , ! 2, 440 
790,12,49630, 12,55210,12,60630,12,661 10, 12,71470, 12,76790, 12,82020* 
812,37260, 12,92420 , 12,97530* 1 3,02580* 13,07590, 13,12550* 13, 17460* 13, 
922320,13,27140, 13,31940, 13,36640* 13,4131 0,13,45950*13,50540* 13,550 
£90,13,59590, 13,54050, 13,68470,14,10500, 14,48970/ 

CONTIN’UE 
!<TEST = 0 
IP=-1 
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> 
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3 

4 

5 

6 


7 

8 
9 


C 

10 
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IF ( IE.EQ.0.0 ) GO TO 5 

IF (ZTHRUST.EQ.O.O) GO TO S 

NLC=NZCNT-1 

IF (IDD8S.EQ.0) GO TO 3 

CALL MJLUP (ZTHRUST«ZEE»2«NLC*50»10« IP^ZEN^YFN) 

GO TO 4 

call MTLUP (ZTHRUST»ZEE»2«NLC»50» 10« IP»ZFN«YFN ) 

DLDB=0.0 

WRITE (6«19) ZTHRUST , (2EE{ I ) • 1 =1 « 1 O ) 

WRITE (8»20> ZEEC 3) »ZEE(1 ) »ZEE(2 ) »ZEE(5) «EMACH»ZEE < 10 ) 

READ (5«ENPUT) 

IF (E0F(5)) 18»5 

H=HJ ( IFUEL ) 

HS=HSvl( IFUEL ) 

TF=TFJ( IFUEL ) 

FS=FSJ( IFUEL 1 
AMP=AMPJ ( IFUEL ) 

DO 7 I =1 »78 

TABH4 ( I ) = TABH4J ( I t IFUEL > 

TABPH5C I ) =TABPH5J ( I » IFUEL) 

TABPSH ( I ) =TA8PSH J { I « IFUEL ) 

TABPSPI 1 )=TABPSPJ ( I « IFUEL) 

IF (IE) 8«8«9 
WRITE (6«ENPUT> 

1 E= IE+1 

N0=0 

NZCNT=1 

AMP28=AMP*28.97 

FS28MP= ( < I .+FS)*28.97-AMP)/FS 

WRITE ( 6 * 23 ) I CASE* I PC »PS ♦ ALT * EMACH * HPEXT » I DDBS » WCOWl * WB • C V » ZTHRUS 
IT . ' 

IPP=0 

EN=ENN 

P80P0=P80P8 

TT4=TT4N 

CALL FRES 

HNEW=ALT 

IC0DE=1 

IF (IDES-1) I0,ll9ll 



B 21 1 
B 212 
B 213 
B ‘214 
B 215 
B 216 
B 217 
B 218 
B 219 
B 220 
B 221 
B 222 
B 223‘ 
B 224 
B 225 
B 226 
B 227 
B 228 
B 229 
B 230 
B 231 
B 232 
B 233 
B 234 
B 235 
B 236 
B 237 
B 238 
B 239 
B 240 
B 241 
B 242 
B 243 
B 244 
B 245 
B 246 
B 247 
B 248 
B 249 
B 250 
B 251 
0 252 
B 253 
B 254 


PT2PT1 =PT2TID 
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c 

WIDE IS INPUT FOR WlD 

B 

255 



W1D=WIDE 

B 

256 



OELTOO=STTO/DELTO*1 ./PTIPTO 

B 

257 



EN=ENN 

B 

258 



A4A4D=1 • 

B 

259 



CALL COND 

B 

260 



CALL TURB 

B 

261 



CALL COOL 

B 

262 



HT6M=HT5AM 

B 

263 



PHIT6M=PHT5AM 

B 

264 



F6=F5A 

B 

265 



TT6=TT5A 

B 

266 



0F6=1 0+F6 

■B 

267 



PT6T5A=PT65AD 

8 

268 



CALL NOZZ 

B 

269 



CALL ENG 

B 

'270 



WRITE (6«2l ) PT2T1D*W1 DCOR « Z »ONSRT « ETA2 ♦ PT 1 PTO ♦ TTO ♦ TT2 ♦ TT4 » DH2M1 »E 

B 

271 



1TA5 

B 

272 



WRITE <6»22> HT1.HT2I 

B 

273 



GO TO 16 

B 

274 


c 


B 

275 

> 

c 

REGULAR CASE 

B 

276 

i 

Oo 

c 


B 

277 


1 1 

CALL MATCH 

B 

278 



CALL TURB 

B 

279 



CALL COOL 

B 

280 



IF <IPC-2) 13,12*17 

B 

281 


12 

HT6M=HT5AM 

B 

282 



PHIT6M=PHT5AM 

B 

203 



F6=F5A 

B 

284 


r 

TT.6=TT5A 

8 

285 



0F6=1 *+F6 

B 

286 



PT6T5A=PT65AD 

B 

287 



CALL NOZZ 

B 

288 



CALL ENG 

B 

289 



GO TO 16 

B 

290 


C 

NO AFTERBURNER 

B 

291 


C 

afterburner 

B 

292 


13 

CALL after 

B 

293 


14 

IF ( IPP-1 ) 2,15,15 

B 

294 


15 

CALL PARP 

8 

295 



GO TO 2 

8 

296 


1 6 

P0WER=TT4 

B 

297 



CALL OUTRUN (POWER) 

B 

298 



non 


GO TO 14 B 299 

17 IPP=1 B 300 

IF (IPC-3) 13*13«12 B 301' 

18 STOP B 302 

C B 303 

19 FORMAT (/30X»58H0UTPUT ENGINE PARAMETERS CORRESPONDING TO THRUST V B 304 

lALUE OF » lF8o 1 ♦5HLBS. /1X» lOHGASFLOW = »1F7«1«8H LBS/SEC * 2X 1 5H JET B 305 
2VEL0CITY = « lF7o2«7KFT/SEC»2X»l IHJET AREA = « 1 F7 . 3 » 7HSQ* FT o « 2X1 1 HR B 306 
3AM DRAG = » 1 F8 . 1 « 5HLBS * /1X13HNET THRUST = » 1 F8 * 1 « SHLBSo »2X19HTUR B 307 
4B1NE IN TEMPO = , 1 F7 o2 « 2HR • 2X 1 8HNO I SE REDUCTION = o 1 F7 o3 » 3HDB* » 2X 1 B 308 
S9HATM0SPHERE = STD + tlF3ol»3H C «/ l X28HSPEC I F I C FUEL CONSUMPTION = B 309 
6 * lF6»4o 1 IHLBM/HR/LBFo 9 I'xaOHTURB INE EXIT TEMP = «lF6.It2H R) B 310 

20 FORMAT < 1 F 1 0 « 4 » 1 F 1 0 « 3 » 1 F 1 0 » 2 9 1 FI 0 « 1 » 1 F 1 0 » 4 9 1 F 1 O » 2 > B 311 

21 FORMAT ( 1 H06X1 1 HDES I GN C ASE5X6HPT2T 1 DE 1 6 . 8 9 4X6HW 1 DCORE 1 6* 8 9 9X 1 HZEl B 312 

1 6»895X5H0NSRTE16.8/6X4HETA2E16.894X6HPTlPT0El6o8f7X3HTT0E16«897X3H B 313 
2TT2E16*897X3HTT4El6*8/6XFHDH2-lE16o896X4H£TA5EI6.8) B 314 

'22 FORMAT ( 1 H06X3HHT 1 5XE 1 6 0 8 9 SX4HHT2 I 4XE1 6 . 8/ ) B 315 

23 FORMAT ( 1 H06X5H I NPUT5X4HCASE4X I 4 9 BX2HPC6X I 4 9 5X2HP36XE 16.8/7X3HALTE B 316 

1 16«8 9 6X4HMACHE16»89 5X5HHPEXTE16.8 9 5X5HIDDBS4XI4/SX5HWC0W1E16.89 8X2 B - 317 
2HWBE16«898X2HCVEI6o8,3X7H2THRUSTE16.2/.) B 318 

END . 8 319- 


SUBROUTINE MATCH C 1 

DIMENSION ALTSAV(15)9 EMSAV<15)9 N0SAV(15)9 FNSAV(15>9 W2SAV(15)9 C 2 

ISFCSAV(IS) C 3 

DIMENSION TABTEM(78)9 TABH<78)9 TABPHI{78>9 TABPSH(78)» TABM{19>9 C 4 

1TABPT1<19>9 TABNST(9)9 TA3PT2<09>9 TABZ(11)9 TABW 1(99)9 TABETA(99) C 5 

2t ANS(5)9 SFCC(25)9 TABPSP(78> C 6 

DIMENSION PBTAB(8)9 HOT (20 >9 DTTAB(20) C 7 

DIMENSION TABPHS(78)9 TABH4 ( 78 ) ' C 8 

DIMENSION TABN(12)9 TABA7(12)9 RAT(2)9 T4(3) C 9 

DIMENSION TABVJ(17)9 TABDLDBC17) C 10 

C 1 I 

INPUT C 12 


C 13 

COMMON /BKl/ ALT 9 EMACH 9 PT 2 TID 9 TT 4 N 9 WCOWI 9WB«A7A7DN9PT65AD9A4A4DN9P C 14 
ia0Pa9WlD»H9ETA209ENN9X9Y9Tr9PT8PT69CV9TT6MAX9PT6TSA9ETA69 1DES*VAL9 C 15 
2GAMO 9 RO 9 A J 9 TSLS 9 PSLS 9 G 9 TA3TEM 9 TA 3 H 9 TABPH I 9 TABPSH 9 T ASPSP « TA8M 9 TABNS C 1 6 
3T9TABPT29TAB29TASW1 9TABETA9TA0PT1 9lDOBS*TA3VJ9TAODL08 C 17 

C 18 


o o 

8| 

sg; 

I? 



c 



n o 


o 


computed values c 

c 

COMMON /BK2/ WBOWI « A0AS« ABC0R»A80A7» A7A7Dt ANS« APO* A4A4D« A7» AStAOf a C 
180W1 .ABEM7»TO»T8»T7«T7I » T8 I »TT0*TT2»TTl «TT2U TT2TT1 tXT5» TT5A « TT7TT C 
21 »TT6»TT7TT6.TT4,P0*P8OP0»PT7»HF«H8M,H8,HSI » H8 I M ♦ H7M ♦ H7 ♦ H7M I C 
3AB5» HO » HTO« HT2 » Hj2 I » HT4M , HT4 » HT5M » HT5 » HT505 * HT5AM . HT5H * HT7M < HT7M I . C 
4HT6M« HT8M»HT5A05, PH IO»PH ITO .PH I T2 I ♦ PH I T 1 » PH I T4M . PH I T4 . PH I T5A . PHTSA C 
5M.PHI8M.PHITBM.PHl8iPHI2I .PHIT6M.PHI7M.PHIT6.PHIT7M.PH0F6.PH0H6.PT C 
60P0.PT1PT0.PT2T1 S.PT2PT1 .PT4PT2 . PTSPT4 . PT5T5A . PT8P8 . PT 1 PT7 »PT7P7. P C 
7T6PT5.PT6T5D.PT5AM .PT55AD.PS IH4.PSPI4.PS1H5.PSP1 5A.PSIH5A.PSPI8.PS C 
8IH8.PS1H7.PSP16.PSIPHI . PS I H6 .PS I HB I » PSP 1 7 ♦ EN . ETA2 . ET5A « ETETS . ET A2S C 
9.ET2D.EM7.EM8.W1C0R.W1 DCOR . W4C0R . W4DCOR . W7C0R . W70COR. W4T0T4 « W4TEST C 
S.Wl .W1K.DH2M1 .STTO.DELTO.DELTOD. GAMAS. GAMA7.V0.V7.V8.GAM7P1 .GAM7M1 C 
$ . ONSRT . SFC . SFCC ♦ 0F4,. 0F6 . F0F6 . F4 . F5A . FAB . F6 . FNOWl .FN . F6A ♦ F0F4 . F05A ♦ C 
SW8.RT5.RTA.RT6.RT5A.CP7.CP8.CAPA8.CAPA7.0WB. WCB. I COMB. IPP. IS. ICOUN C 
ST. IKONT. JCONT. IC.KTE'ST.THETI .THETO.DWl .0W2.DDHP.Z1 . 22 . Z3 . Z . ETAC . ET C 

SAO. FN0W1D.HT6.PH1 H6 . T5 A . T4 . R AT . WCOW 1 C . PBT AS . HDT . 0TTA8 . IPS.DTEMP.ET C 

SAP . ALTSAV .EMSAV .NOSAV . FNSA V . W2SAV . SFCSAV . NO . I NOZZ . TT2T1D. HT5A .PHT 1 C 
S.TT4TT2.PHIT5.P5PI5.PHIT5M.HT1 . PT2T 1 L .PS . TABA7 . T8VAL . T7VAL . NL I M .NU .C 
SMBER.NSTART . NF I NAL .AVAL . XMNSLU . TABN ♦ ETA4 . ETAS. TABPH5 . TA0H4 . HPEXT. I C 


SCASE.IPC.PT4T2D.HS.T5VAL.WC ' C 

ICOUNT=0 C 

JC0UNT=1 C. 

THET1=THET0 C 

0NSRT=EN/SQRT (THeTI ) C 

Z=1 . C 

CALL 01 SCOT (ONSRT. ONSRT. TABNST. TABPT2. TABPT2. -1 1 .09.0. PT2T1L) C 

PT2TIL=PT2T1L*PT2T1 S C 

CALL DISCOT (ONSRT . Z . TABNST . TABW 1 . TABZ . 1 1 . 99 « 1 1 « W 1 COR ) C 

CALL DISCOT (ONSRT . Z , TABNST . TABETA . TABZ . 1 1 . 99 . 1 1 .ETETS > C 

PT2PT1= (Z* ( 1 .-X ) + X)* (PT2T1L-1 » )+l . C 

ETA2=ETETS1<-ETA2S - C 

IF (PT2PT1 .LT.O ) GO TO 1 1 C 

PHIT2I=PHT1+R0/AJ#AL0G(PT2PTI ) C 

CALL DISCOT (PH I T2 I . PH I T2 I . TABPH 1 . TABTEM . TABTEM . - 1 1 .78 . 0 . TT2 1 ) C 

CALL DISCOT ( TT2 I .TT2 I . TABTEM .TABH. TABH . - 1 1 . 76 . 0 . HT2 I ) C 

DH2M1= (HT2I-HT1 )/ETA2 C 

HT2=HT1+DH2M1 C 

CALL DI SCOT ( HT2 . HT2 .TABH . TABTEM . TABTEM .-1 I . 78 .0 .TT2 ) C 

TT2TT1=TT2/TT1 C 

WlK=WiCOR/Wl DC0R^^W1D■^(•DELT0D C 

W1=W1K*DELTO/STTO*PT1PTO C 

W30WI=W8/W1 C 

0WS=Io-WB0Wl C 


19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 



u ro 


I 



5 

6 


7 

a 

9 

10 


1 1 


c 

12 

13 

14 


CALL COMB 

DW2=W4DC0R“W4C0R»1 «/A4A4D 
ABCOR=ABS{DW2 ) 

W4TEST= <.OOl«W4DCOR 

IF (JCOUNT-2) 2»8»7 

IF (DW2) 4»5*3 

IF ( A8C0R-W4TEST > 5»5«6 

WRITE (6«12> W4DCOR*W4COR 

IC0UNT=2 

Z = 22 

RETURN 

JCOUNT= JCOUNT+1 
DWl =DW2 
22=0 
21=2 
2 = 22 
GO TO 1 

IF (ABCOR-W4TEST) 5»5«10 
IF (DW2) 9»5«4 
IF (ABCOR-W4TEST) l0,5»10 
23=21“ ( ( (21-22)*DW1 >/(DWl-DW2) ) 
IF (23«LToO) 23=o5'«-22 
IF (Z3oGT«l.) 23 =o5-»Z1 
21=22 
22=23 
2 = 22 
DW1=DW2 

JC0UNT=JC0UNT+1 

IF (JCOUNT-25) 1,1,4 

WRITE (6«13) PT2PT1 ,X»2,PT2T1L 

WRITE (6,14) 21 ,Z2,23,DW1 ,DW2 

STOP 


O O 
O 

Sg 

tr' 

i> 
Q 
tg 


jO 

I 


CO 


FORMAT (1H06X11HNO SOLUT I ON6X6HW4DCOR2XE 1 6o8 , 5X5HW4C0R3XE 1 6« 8/ ) 
FORMAT nH06X4E15o8) 

FORMAT ( lH05X2H2lE15.8,2X2H22ElS«a,2X2HZ3E15o8,2X3HDWlE15eB,2X3HDW 
12E15,8/ > 

END 


C 63 
C 64 
C 65 
C 66 
C 67 
C 68 
C 69 
C 70 
C 71- 
C 72 
C 73 
C 74 
C 75 
C 76 
C 77 
C 78 
C 79 
C 80 
■C 81 
C , 82 
C 83 
C 84 
C 85 
C 86 
C 87 
C 88 
C 89 
C 90 
C 91 
C 92 
C 93 
C 94 
C 95 
C . g -6 ' 
C 97 
C 98 
C 99 
C lOO- 



non o n o 
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SUBROUTINE NOZZ D 

dimension ALTSaV< 15)« EMSAVU5)* NOSAVdS). FNSAV(15)« W2SAV(!5)« D 
ISFCSAVdS) D 

dimension TABTEM(78)» TABH<78)i TABPHn78)» TABPSH(78)* TABM(19)t D 
lTABPTld9). TABNST(9). TA8PT2(09)« TABZ(ll). TABW1{99)< TA8ETA(99) O 
a, ANS<5). SFCC(25)« TABPSP{78) D 

DIMENSION P8TAB(8)» HDT(20>* DTTAB(20> O 

DIMENSION TABPH5(78>» TABH4(78) D 

DIMENSION TABN(12)» TABA7(12)» RAT(2)» T4(3) ’ O 

DIMENSION TABVJ(17)» TABDLDBd?) D 

D 

INPUT D. 

O 

COMMON /BKl/ ALT*EMACH*PT2tlD«TT4N*WC0Wl ♦WB»A7A7DN«PT65AD«A4A4DN»P D 
180P8«W10«H.ETA2D»ENN»X»Y,TF»PT8PT6»CV»TT6MAX.PT6T5A«ETA6» lOEStVAL. D 
2GAM0«R0»AJ«TSLS«PSLS»G»TA8TEMtTA8H«TABPHl , TA8PSH » TABPSP » TA8M f TABNS D 
3T« TABPT2»TABZ»TABW1 * TABETA » TABPT 1 , IDDB£«TABVJ» TABDLDB D 

D 

computed VALUES D 

D 

COMMON /BK2/ WBOWI » AOAS»A8COR«A80A7»A7A7D»ANS«Ap 6»A4A4D»A7»A8»AO»A D 
1 8OWl‘«ABEM7.T0d8,T7.T7I ♦ T8 I »TT0.TT2»TT1 tTT2I .TT2TT1 dT5, TT5A * TT7TT . D 
21 f TT6»TT7TT6»TT4»P0»P80P0»PT7*HF<H8Mf HBtHSI » H8 I M i H7M » H7 » H7M I .HABrH D 
3AB5«HO»HTO»HT2»HT2I ♦ HT4M »HT4 » HT5M » HT5 • HT505 » HT5AM » HT5H « HT7M ♦ HT7M 1 « D 


4HT6Mt HTBM.HT5A05,PHl0»PHIT0t PH I T21 t PH IT 1 « PH I T4M » PH I T4 ♦ PH I T5A « PHT5A D 
5M.PHiaM.PHIT8M»PHia»PHI2 I «PHIT6MtPHI7M»PHIT6<PHIT7M»PH0F6«PH0H6«PT ' D 
60PO«PTIPTO«PT2T1S»PT2PT1 .PT4PT2 ♦ PT5PT4 * PT5T5A » PT8P8 ♦ PT 1 PT7< PT7P7» P D 
7T6PT5«PT6T5DfPT5AM,PT55AD»PS IH4»PSPI4»PSIH5*PSP15A»PSIH5A»PSPI8*P‘= D 
8IH8»PSIH7»PSPI6«PSIPHI *PS I H6 »PS I H8 I , PSP I 7 »EN i ETA2 »ET5A ,ETETS fETAES D 
9 « ET2D » EM7 » EM8 ♦ W 1 COR * VJ 1 DCOR ♦ W4C0R ♦ W4DC0R » W7COR » W7DC0R » W4T0T4 » W4TEST D 
S,W-1 »W1K»DH2M1 tSTTOfOELTOiOELTODt GAM A8 . GAMA7 « VO « V7 « V8 » GAM7P1 »GAM7M1 D 
$»ONSRTiSFC»SFCC«OF4 »0F6«F0F6»F4tF5A»FAB»F6»FN0Wl » FN « F6A « F0F4 ♦ F05A » D 
$W8«RT5»RTA»RT6»RT5A«CP7^CP8»CAPA8»CAPA7»0WB»WCB« ICOMB* I PP d S « ICOUN D 
$TdKONT» JCONT* I C » KTEST » THET 1 »THET0»DW1 »DW2«DDHP.Z1 ♦ Z2 » Z3 ♦ Z » ETAC »ET D 
$AO »FNOW 1 D,HT6 ♦ PH I H6 » T5A d4» RAT » WCOW 1 C»PBTAB»HDT»DTTAB» I PS » DTEMP .ET D 
$AP*ALTSAV,EMSAV»N0 SAViFNSAV»W 2SAV«SFCSAV»N0» I NOZZ» TT2T1D«HT5A,PHT1 D 
$ » TT4TT2 * PH I T5 iPSp I 5 »PH I T5M . HTl *PT2T 1 L »PS . TABA7 < T8VAL ♦ T7VAL » NL I M tNU D 
$MBER»NSTART»NFINAL»AVAL»XMNSLU*TABN»ETA4»ETA5*TABPH5»TABH4»HPEXT* I D 


SCASEvIPC»PT4T2D«HS,T5VAL D 

COMMON /RRR/ AMP28»FS28MP, AMP D 

PHIT7M=PHIT6M D 

IIK0NT=0 D 

F0F6 = F6/d .+F6) D 


PT8P8 = PT0P0i(-PTlPT0*PT2PTl*PT4PT2’ifrPT5PT4*l • /PT5T5 A*PT8PT6*PT6T5A'* ( 1 D 
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CD -J O' U1 -£!■ u ro 


1 


U) 


9 


10 


1 1 
12 

13 


1 O /P80P8 ) 

IF (PT8P8-1 .005 ) 1 .1.2 
WRITE (6.31 ) 

GO TO 8 

IF (F6-F5A) 3.4.3 

RT6= 1545 .*F0F6* (FS28MP+AMP/F6 ) /AMP28 

GO TO 5 

RT6=RT5A 

TT6=TT5A 

PHIT8M=PHIT6M 

PHI 8 M = PHIT 8 M- (RT6/AJWAL0G(PT8P8 ) ) 

IROU=8 

call ITT (PHI8M.TABPH5.TABPHI .TABTEM.TABPSP.T8VAL.F0F6. I R0U.T8. PH I 
1 8 . PSP i 8 ) 

call DISCOT (T8.T8.TABTEM.TABH.TABH.-1 1 ,78,0. H8) 

CALL DISCOT ( TS . T8 . TABTEM . T ABPSH . TABPSH . - 1 1 . 78 . O .PS I H 8 ) 

H 8 M=H 8 +PS IH8^<-F0F6 
HT 8 M=HT 6 M 

IF (HT 8 M-H 8 M) 6.6.9 
WRITE (6.32) 

IF (.IDES) 7.7.8 

STOP 

A4A4D= 1 * 

CALL MAIN 

V8=SQRT (2.*G*AJ* CHT8M-H8M> ) 

T8I=T8+2. 

call DISCOT (T8I .T81 .TABTEM. TABH.TA8H.-1 1 ,78.0. H8I ) 

call DISCOT (T8I .T8I .TABTEM, TABPSH. TABPSH.-l 1 ,78.0. PSIH8I ) 

H8 1 M=H8 I +PS I H8 I ^F0F6 
CP 8 = (H 8 IM-H 8 M )/ (T 8 I-T 8 ) 

W8=W1*0F6#0WB 

CAPA 8 = ( W 8 *RT 6 *T 8 >/ < V8*P80P0*P0 ) 

GAMA 8 =CP 8 / ( CP 8 -RT 6 /A J ) 

A 8 =SQRT (GAMA 8 *RT 6 *T 8 ^G ) 

PT1PT7= (PT2PT1-»PT4PT2*PT5PT4*1 . /PT5T5A*PT6T5A»PT8PT6 ) #* ( - 1 ) 

TT7TTl=TT2TTl*TT4/TT2*TT5/TT4#TT5A/TT5-!<-TT6/TT5A 

IF { IDES-1 >10.11,11 

W7DC0R = W1DC0R*( 1 ,-WBOWl > #0F6*SQRT (TT7TT 1 li^PTlPTT 
A7A7D=1 • 

GO TO 12 

W7C0R = WIC0R*{ 1 .-WBOIVI } *0F6*SQRT ( TT7TT1 )»PTlPT7 

EM8=V8/A8 

IF (EM 8 -I 0 ) 13,13.16 
A80A7=1 » 
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D 

48 

D 

49 

D 

50 

D 

51 

D 

52 

D 

53 

D 

'54 . 

D 

55 ■ 

D 

56 

D 

57 

D 

■58 

0 

59 

D 

60 

D 

61 

D 

62 

D 

63 ■ 

0 

64 

D 

65 

D 

66 

D 

67 

D 

68 

0 

69 

D 

70 

0 

71 

D 

72 

D 

73 

D 

74 

D 

75 

D 

76 

D 

77 

D 

78 

D 

79 

D 

80 

D 

81 

0 

82 

D 

83 

D 

84 

0 

85 

D 

86 

D 

87 

D 
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EM7=EM8 D 89 

GAMA7=sGAMA8 D 90 

V7=V8 O 91 

GAM7P1=GAMA7+1 . D 92 

GAM7M1=GAMA7-1 • 0 93 

A0AS= (2./GAM7P1 )»« ( GAM7P 1 / < 2 o*GAM7M 1 ) )*1 •/EM7*U .+GAM7M1 /2 • *EM7**2 D 94 

1 )**(GAM7P1/(2.*GAM7M1 n D 95 

IF (lDES-1) 14«15»15 O- 96 

14 AOASO=AOAS D 97 

RETURN D 98 

15 A7A7D=(W7C0R*A0AS)/(W7DC0R*A0ASD) D 99 

RETURN D 100 

16 A0AS=1. D 101 

PT7P71=1.8 D 102 

PT7P7=PT7P71 D 103 

17 CONTINUE D 104 

IF <PT7P7) 18*19,19 D 105 

18 WRITE (6*33) PT7P7 D 1 06 

CALL MAIN D 107 

19 CONTINUE D 108 

PHI7M=PH1T7M~RT6/AJ#AL0G (PT7P7 ) D 109 

IR0U = 7 D 1 1 0 

CALL ITT {PHI7M*TABPH5*TABPH1 ,TABTEM,TABPSP.T7VAL«F0F6, IR0U*T7,PHI Dill 

17,PSP17) D 112 

CALL OISCOT (T7«T7,TABTEM,TABH,TABH*-1 1 *78*0*H7) D 113 

CALL DISCOT < T7 * T7 * T ABTEM , TABPSH * TABPSH * - 1 l * 78 * 0 * PS I H7 ) O 114 

H7M=H7+F0F6*PS IH7 D 115 

HT7M=HT6M O 116 

IF (HT7M-H7M) 20,21,21 O 117 

20 WRITE (6,34) D 118 

PRINT 35, H7M,HT7M D 119 

IF (IDES) 7,7,8 D 120 

21 V7=SQRT (2,»G*AJ-» (HT7M~H7M) ) D 121 

T7I=T7+2« D 122 

call DISCOT (T7I ,T7 I ,TABTEM ,TABH,TABH,-1 1 ,78,0 ,H7I ) D 123 

CALL DISCOT < T7 I , T7 I , TABTEM , TABPSH, TABPSH , -1 1 , 78 , 0 , PS IH7 I ) D, 124 

H7IM=H7I+F0F6*PSIH7I D- 125 

CP7=<H7M-H7IM)/(T7-T7I ) D 126 

GAMA7=CP7/(CP7-RT6/A J) . D 127 

A7=SQRT (GAMA7*RT6*T7-»-G ) D 128 

£M7=V7/A7 D 129 

I IKONT= I IKONT+1 D 130 

IF (II!<0NT--2) 22,27,27 D 131 

EM71=EM7 D 132 


22 



A-15 


IF (EM7.GE*o999«AND*EM7oLE*l .001 ) GO TO 30 
IF <EM7~1.) 23*30.26 

23 PT7P72=PT7P71 +VAL 

24 I 1K0NT= 1 IKONT+I 
PT7P7=PT7P72 

25 IF (IIKONT-10) 17*17*29 

26 PT7P72=PT7P71-VAL 
GO TO 24 

27 IF (EM7*GEo .999oAND.EM7.LE. 1 »00l > GO TO 30 
IF (EM7 .lt. .999. OR. EM7.GT. 1 .001 ) GO TO 28 


28 EM72=EM7 

PT7P73=PT7P71- ( (EM71~1 . )/ ( EM71 -EM72 )» (PT7P71 -PT7P72 ) ) 

PT7P71=PT7P72 
PT7P72=PT7P73 
EM71 =EM72 
PT7P7=PT7P73 
I I k0NT= I IKONT + 1 
GO TO 25 

29 WRITE (6*361 

PRINT 37* T7*H7*PSIH7*H7M*V7*HT7M*PHI7.PSPI7*F0F6*T7I *H7 I * PS I H7 I » H 
17IM,.A7»CP7*GAMA7*EM7*PT7PT3« I IK0NT»PT7P7 

30 CAPA7= (Wl *0F6^^RT6^^T7 > / C V7* ( PT8P8* I ./PT8PT61M . /PT7P711'P80PO-»P0 ) )1«-0W8 
A80A7=CAPA8/CAPA7 

IF (IDES) 14*14*15 


C 

31 FORMAT (1 

32 FORMAT ( 1 

33 FORMAT { 1 

34 FORMAT (1 

35 FORMAT (1 

36 FORMAT (1 

37 FORMAT ( 1 

1F12.5*2X* 
2*2X*3HH7I 
3*3HCP7*F1 
41 IKONT* 14 

END 


H06X37HNOZZLE PRESSURE RATIO LESS THAN 1.005) 

HO , 6X ♦ *HT8M I SLESSTHANH8M# ) 

HO *6X**PT7P7I SNEGAT IVE=-S-*EI 4.6 ) 

HO * 6X * *HT7M I SLES3THANH7M1<- ) 

H * 2X * 3HH7M , F 1 2 o 5 ? 2 X * 4HHT7M . F 1 2 . 5 ) 

H06X24HSTOPPED IN M7 ITTERATION) 

H *2X*2HT7«F12o5*3X«2HH7«F12.5*3X*5HPSIH7*Fl2<i5» 3X.3HH7M* 
2HV7«F12«5*2X*4HHT7M*F12c5*2X* 4HPH I 7 * F 1 2 o 5/ 1 X* 3HT7 I *F12.5 
*F12.5*2X*SHPSIH7I * FI 2 . 5 * 2X * 4HH7 I M * FI 2 . 5 ♦ 2X *2HA7 * F 1 2 .5 * 3X 
2 o 5 . 1 X * 5HGA MA7 * F 1 2 . 5/ 1 X * 3HEM7 * F 1 2 . 5 * 1 X * 4HF0F6 « F 1 2 . 5 * 2X * 5H 
* 1 X * 5HPSP I7*F12o5*lX* 5HPT7P7 * F 1 2 . 5 « 1 X ♦ 6HPT7PT3 » F 1 2 . 5 ) 


D 133 
D 134 
D 135 
D 136 
D 137 
D 138 
D 139 
D 140 
D I4l 
D 142 
D 143 
D 144 
D 145 
D 146 
D 147 
D 148 
D 149 
D 150 
D 151 
D 152 
D 153 
D 154 
D 155 
D 156 
D 157 
D 158 
D 159 
D 160 
D 161 
D 162 
D 163 
D 164 
D 165 
D 166 
D 167 
D 168 
O 169- 




subroutine COOL E 1 

dimension AlTSAVUS)* EMSAVdS)* N0SAV(15>* FNSAVdS)* W23AV(15)* E 2 





C 

C 

C 


C 

c 

c 


1SFCSAV( 15) 

DIMENSION TABTEM(78)* TABH(78>9 TABPHI(78>9 TABPSH<78)» 
1TABPT1(19)9 TA0NST(9)9 TABPT2(09)« TABZ(11)« TABW1(99)» 
2 9 ANS(5)9 SFCC(25>9 TABPSP(78) 

DIMENSION P3TAB <8 ) 9 HDT(20)9 DTTAB(20) 

DIMENSION TABPH5(78)9 TA8H4(78) 

DIMENSION TABN(12)9 TABA7(12)9 RAT (2) 9 T4(3) 

DIMENSION TABVJ(17)9 TABDLDB(17) 


TABM(19)9'.: 

TABETA(99) 


INPUT 

COMMON /BKl/ ALT 9 EMACH 9 PT 2 TI D9TT4N9WC0WI 9 W8 9 A7A7DN 9 PT65AD 9 A4 A4DN9 P 
I80P89W1D9H9ETA2D9ENN9X9Y9TF9PT8PT69CV9TT6MAX9PT6T5A9ETA69 IDES 9 VAL 9 
2 GAMO 9 RO 9 AJ 9 TSLS 9 PSLS 9 G 9 TABTEM 9 TABH 9 TABPHI 9 TABPSH 9 TABPSP 9 TABM 9 TABNS 
3T9TABPT29TAB29TABW1 9 TABETA 9 T ABPT 1 9 I DDBS 9 TABVJ 9 TABDLDB 


COMPUTED VALUES 

COMMON /BK2/ WBOW 1 9 AOAS 9 A 8 COR 9 A 8 OA 79 A7A7D 9 ANS 9 APO 9 A 4 A 4 D 9 A7 9 A8 9 AO 9 « 
180W1 9 ABEM 79 TO 9 T 89 T 79 T 7 I 9T8I 9 TTO 9 TT 29 TTI 9 TT 2 I 9TT2TT1 9 TTS 9 TT5A 9 TT7TT 
21 9TT-69TT7TT69TT49P09P80P09PT79HF9H8M9H89H8I 9 H8 I M 9 H7M 9 H7 9 H7M I 9 HAB,H 
• 3 AS 59 HO 9 HTO 1 HT 29 HT 2 I 9HT4M9HT4 9HT5M9HT5 9HT505 9HT5AM9HT5H9HT7M9HT7MI 9 
4 HT 6 M 9 HT 8 M 9 HT 5 A 059 PH IO 9 PH ITO 9 PHI T2I 9 PHITI 9 PH I T4M 9 PH I T4 9 PH I T5A 9PHT5A 
5M 9 PH 1 8M I PH I T8M 9 PH I 8 9 PH I 2 I 9 PH I T6M 9 PH I 7M 9 PH I T6 9 PH I T7M 9 PH0F6 9 PH0H6 9 PT 
6 OPO 9 PTIPTO 9 PT 2 TIS 9 PT 2 PTI ♦PT4PT29PT5PT49PTST5A9PT8P89PT1PT79PT7P79P 
7T6PT59PT6T5D9PT5AM9PT55AD9PSIH49PSPI49PSIH59PSP15A9PSIHSA9PSPI89PS 
8IH89PS1H79PSPI69PSIPHI 9 PS I H6 9 PS I HB I 9 PSP I 7 9 EN 9 ETA2 9 ET5A 9 ETETS 9 ET A2S 
9 9 ET2D 9 EM7 9 EM8 9 W 1 COR 9 W1 DCOR 9 W4C0R 9 W4DC0R 9 W7C0R 9 W7DC0R 9 W4T0T4 9 W4TEST 
S9 W1 9W1K9DH2M1 9 STTO 9 DELTO 9 DELTOD 9 GAMA 89 GAMA 79 VO 9 V 79 V89GAM7P1 9GAM7M1 
$90NSRT9SFC9SFCC90F490F69P0F69F49F5A9FAB9F69FN0W1 9 FN 9 FftA 9 FOF4 9F05A 9 
$W89RT59RTA9RT69RT5A9CP7iCP89CAPA89CAPA790WB9WCB9 I COMB 9 I PP 9 IS 9 ICOUN 
$T9 IKONT 9 JCONT 9 IC 9 KTEST 9 THETT 9 THETO 9 DWI 9 DW 29 DDHP 921 9229239Z9ETAC9ET 
$ AO 9 FNOW 1 D 9 HT6 9 PH I H6 9 T5A 9 T4 9 RAT 9 WCOW I C 9 PSTAB 9 HD T 9 QTTAB 9 I PS 9 DTEMP 9 ET 
$AP9AlTSAV9EMSAV9N0SAV9FNSAV9W2SAV9SFCSAV9N0 9 IN0ZZ9TT2T1D9HT5A 9PHT1 
S,TT4TT2 9 PHIT 5 9 PSPI 5 9 PHIT 5 M 9 HT 1 9 PT2T 1 L 9 PS 9 TABA7 9 T 8 VAL 9 T7VAL 9 NL I M 9 NU 
SMBER 9 NSTART 9 NF I NAL 9 AVAL 9 XMNSLU 9 TABN 9 ETA4 9 ETA5 9 TA8PH5 9 TABH4 9 HPEXT 9 I 
SCASE9 IPC9PT4T2D9HS9T5VAL 
COMMON /RRR/ AMP 289 FS 28 MP 9 AMP 

HT5AM= ( HT5M*WCB*0F4-)-HT2*WC0W1 C ) / (0F4»WCB+WC0W1 C ) 

F5A=F4*WCB/0WB 
F05A=F5A/ { 1 .+F5A ) 

IROU=6 

CALL ITT (HT 5 AM 9 TABH 49 TABH 9 TA 8 TEM 9 TA 8 PSH 9 T 5 VAL 9 FO 5 A 9 IROU 9 TT 5 A 9 HT 5 A 

1 « oc t A \ 


E 3 
E 4 
• E ' 5 

E 6 
E ■ 7 

E 8 
E .9 
E 10 
E 1 1 
E 12 
E 13 
E 14 
E 15 
E 16 
E 17 
E 18 
E '19 
E 20 
E 21 
E 22 
E 23 
E 24 
E 25 
E 26 
E 27 
E 28 
E 29 
E 30 
E 31 
E 32 
E 33 
E 34 
E 35 
E 36 
E 37 
E 38 
E 39 
E 40 
E 41 
E 42 
E 43 
E 44 
E 45 
F 46 



o n o 


I CALL DISCOT ( TT5A . TT.5A » TABTEM*TABPHI ♦ TABPH! » -11 » 78 * 0 »PH I TSA > 

call DISCOT (TT5A»TT5A*TA8TEM»TABPSP*TABPSP«~1n1 ♦78»0»PSPI5A ) 
RT5A=1545,-»F05A* (FS28MP+AMP/F5A )/AMP28 
PHT5AM=PH I T5A+PSP ISA-K-FOSA 
PTST5A=1. 

return 

END 


■E 47 
E 48 
E 49 
E 50 
E 51 
E 52 
E 53- 


5 = 

I 


C 

c 

C 


SUBROUTINE ITT (PHOM » H45TAB » PHOT AB« TEMTAB 
I »PHO«PSIO) 

DIMENSION ALTSAV(15)» EMSAV<15)« NOSAVU5 
l.SFCSAVI 15) 

DIMENSION TA0TEM(78)» TABH(78)» TABPHI(78 
1TABPT1(19>« TA8NST(9)« TABPT2<09)» TABZ ( 1 
2« ANS(5)f SFCC(25)» TABPSP(78) 

DIMENSION PBTAB(8). HDT ( 20 ) ♦ DTTAB(20) 
DIMENSION TABPH5(VS>, TASH4(78) 

DIMENSION TABN(12)« TABA7(12)» RAT(2)» T4 
DIMENSION TABVJ(17)» TABDLDB<17) 


♦PSITAB.TVAL 

) « FNSAV ( 15 ) 

)♦ TABPSH(78 
1 > « TABWl (99 


(3) 


,FOFF« I ROU «T3 

, W2SAV(15>* 

) • TABM ( 19 ) « 
)» TABETA(99) 


INPUT 

COMMON /BKl/ ALT*EMACH»PT2T1D»TT4N» WCOWl » WB » A7A7DN ♦ PT65AD ^ A4 A4DN * P 
1 80P8* Wl D«H»ETA20»ENN<X« Y»TF»PT8PT6»CV«TT6 MaX»PT 5T5A«ETA6» IDES.VALt 
2GAM0 . RO ♦ A J * TSLS » PSLS < G « TABTEM tTABH , JABPH I « TABPSH « TABPSP » T ABM « TABNS 
3T*TABPT2»TA8Z»TABW1 ♦ TABETA * T ABPT 1 » IDDBS«TABVJ?TABDLDB 


COMPUTED values 


common /BK2/ WBOWl . AOAS« ABCOR ♦ A80A7» A7A7D ♦ ANS i APO« A4A4D» A7 » A8 . AO , A 
180W1 , ABEM7»T0»TS«T7»T7I ♦ T8 I ,TT0»TT2»TT1 .TT2I »TT2TT1 ♦ TT5 ♦ TT5A , TT7TT 
2 1 » TT6 » TT7T T6 « TT4 « PO ♦ PBOPO « PT7 ♦ HF * H8M * H8 « H8 I f H8 I M « H7M » H7 ♦ H7M I ♦ HAS » H 
3A85»H0»HT0»HT2«HT2I ♦ HT4M » HT4 » HT5M » NTS » HT505 » HTSAM ♦ HT5H* HT7M ♦ HT7M I « 
4HT6M t HT8M , HT5A05 , PH I O * PH I TO t PH I T2 I ♦ PH I T 1 * PH I T4 M i PH I T4 * PH I T5A , PHT5A 
5M«PHI8M»PHIT8M»PHIS*PHI2I »PHlT6MtPHI7MtPHIT6«PHIT7M»PHOF6»PHOH6»PT 
60P0»PT1PT0 iPT 2T1S«PT2PT1 »PT4PT2 t PTSPT4 * »T5T5A , PT8P8» PT 1 PT7 *PT7P7 * P 
7T6PT5 » PT6T5D ^ PT5AM « PT55AD* PS I H4 » PSP I 4 »PS I H5 » PSPl SA » PS I H5A » PSp I 8 »PS 
3IH8«PSIH7 qPSPI6«PSIPHI o PS I H 6 « PS I HS I ♦ PSp I 7 » EN » ETA2 » ET5A ♦ ETETS » ETA2S 
9 1 ET2D1 EM7 OEMS o W 1 COR o W 1 DCOR o W4C0R o W4DC0R » W7C0R o W7DC0R* W4T0T4 o W4TEST 
$0 WI o'tf.'!KoDH2Ml oSTTOoDELTOoDELTODoGAMASoGAMA 7» VOoV7» V3,GAM7P1 o.GAMTMl 


F 1 
F -2 
F 3 
F 4 
F 5 
F ' 6' 
F 7 
F -8 
F 9 
-F 10 
E 1 L. 
■F.. 12 
F 13 
F 14 
F 1 5 
F 16 
F 17 
F 1 8 
F 19 
F 20 
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F 22 
F 23 
F 24 
F 25 

F 26 
F 27 
F 28 
F 29 
F 30 
F 31 
F 32 
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S.ONSRT»SFC»SFCCiOF4iOF6,FOF6»F4,F5A»FABtF6«FNOWl ♦ FN ♦ F6A ♦ F0F4 » F05A . 
SW8 »RT5 1 RTA «RT6 «RT5A » CP7 » CPS * CAPA8 , CAPA7 ♦ O WB ♦ WCB « I COMB , I PP « I S « I COUN 
$T» IKONT* JC0NT» I C » KTEST » THET 1 ♦THETOtDWl »0W2«DDHP.Z1 »Z2»Z3« ZtETAC «ET 
$A0»FN0WlD»HT6»PHIH6«T5A»T4«RATtWC0WlC»PBTAB*HDT»DTTAB» IPS*DTEMP«ET 
$AP.AlTSAV«EMSAV»N 0SAV»FNSAV»W2SAV»SFCSAV»N0» IN0ZZ»TT2T1D»HT5A*PHT1 
S.TT4TT2»PHIT5«PSPI5*PHIT5M»HT1 . PT2T 1 L • PS , TAB A7 » T8VAL ♦ T7VAL» NL I M «NU 
$M8ER ♦ NSTART » NF I NAL « A VAL « XMNSLU. TABN » ETA4 ♦ ETA5« TABPH5 » TABH4 tHPEXT t I 
SCASEi IPC»PT4T2D«HS»T5VAL 

dimension PH0TAB(1)» TEMTABCII* PSITAB(1)» H45TAB(1) 

CALL DISCOT (PHOM tPHOM » H45TAB . TEMTAB » TEMTAB » -1 I ♦ 78 »0 *T1 ) 

KONT=0 
T0 = T1 

CALL DISCOT ( TO ♦ TO * TEMTAB t PH0TA8 ♦ PHOTAB « *- II » 78 « 0 ♦ PHO ) 
call DISCOT (T0»T0»TEMTAB.PS!TAB»PSITAB»-1 1 «78«0,PSlO) 
RATlO=PHOM/(PHO+PS I O^FOFF ) 

kont=kont+i 

IF (KONT-2) 2»3»5 
Rl=RATIO 
T2=T1*R1 • 

TO=T2 
GO TO 1 
R2=RATI0 

T3=( 1 «-R2 )/(R2-Rl )» {T2-T1 )+T2 

TO = T3 

GO TO 1 

R3=RATI O 

TPLUS=1 *+TVAL 

TM1N=1 *-TVAL 

IF (R3-TPLUS) ,6«9«7 

IF (R3-TMIN) 7.9*9 

R1 -R2 

R2=R3 

T1 =T2 

T2=T3 

IF (KONT-NLIM) 4.8.8 
WRITE (6* 10) IROU 
RETURN 

FORMAT (1H06X16HDID NOT C0NVERGE3X I 4 ) 

END 
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SUBROUTINE OUTRUN (POWER) 

dimension ALTSAV( 15) t EMSAV(15)» N0SAV(I5). FNSAV(15)» W2SAV(15)* 
1SFCSAV( 15 ) 

DIMENSION TABTEM( 78) * TABH C 78 ) » TABPH[(73)» TABPSH(78)» TABM ( 19 ) » 
1TABPT1U9)» TABNST(9)* TABPT2(09)» TABZ(11)» TABW1(99)* TABETA(99) 
2» ANS(5)» SFCC<25)* TABPSP C 78 ) 

DIMENSION TABPH5(78)f TABH4 ( 78 ) 

DIMENSION TABN(12)« TABA7(12)» RAT(2). T4(3) 

DIMENSION PBTA8<8)» HDT<20)« DTTAB(20) 

DIMENSION ZFN<50), YFN ( 50 . 1 0 ) 

DIMENSION TABVJ(17)» TABDLDBUT) 

COMMON /ZPARAM/ ZFN f YFN i NZCNT ♦ ZTHRUST 

Input 

COMMON /BKl / ALT«EMACH»PT2T 1 D « TT4N » WCOW I » WB » A7A7DN » PT65AD « A4A4DN ♦ P 
1 80P8 « W 1 D » H « ETA2D , ENN « X « Y « TF « PT8PT6 * C V ♦ TT6MAX » PT6T5 A » ET A6 ♦ I DES ♦ VAL » 
2GAM0iR0tAJ.TSLS*PSLS.G»TABTEM*TASH*TABPHI » TABPSH » TABPSP ♦ T ABM » TASNS 
3T»TABPT2»TABZ»TA3W1 * TABETA t T ABPT 1 , I DDBS , T ABVJ » TABDLDB 

COMPUTED values 

COMMON /BK2/ WBOWI , AOAS < ABCOR « A80A7 « A7A7D » ANS * APO « A4A4D ♦ A 7» A 8 * A 0 1 A 
1 80W1 » ABEM7» TO«T8, T7 » T7 I , T8 I ♦ TTO i TT2 » TTl , TT2 I » TT2TT 1 »TT5»TT5A ,TT7TT 
21 .TT6»TT7TT6<TT4,P0»P80P0»PT7«HF,H8M,H8,HBI « H8 I M « H7M t H7 « H7M I ,HAB«H 
3A65»H0*HT0»HT2»HT2 I » HT4M ♦ HT4 * HT5M i HT5 . HT505 * HT5AM t HT5H » HT7M , HT7M I , 
4HT6M , HT8M ♦ HT5AQ5 , PH I O ♦ PH I TO , PH I T2 I , PH I T 1 i PH 1 T4M , PH I T4 ♦ PH I T5A , PHT5A 
5M»PHI8MiPHlT8M*PHI8»PHI2I , PH I T6M ♦ PH I 7M , PH I T6 * PH I T7M , PHOF6 » PH0H6 « PT 
60PO<PT1PTO^PT2T1 S *PT2PT1 . PT4PT2 , PT5PT4 , PT5T5A , PT8P8 « PT 1 PT7 . PT7P7 , P 
7T6PT5«PT6T5D.PT5AM«PT55AD«PSIH4,PSPI4tPSIH5iPSPl 5A,PSIH5A»PSPI8»PS 
alH6.PSIH7»PSP16«PSIPHI «PSIH6 iPSIH8I » PSP I 7 ,EN » ETA2 , ET5A , ETETS . ET A2S 
9 «ET2D« EM7 »EM8« W 1 COR » W 1 DCOR ♦ W4C0R ♦ W4DC0R » W7C0R « W7DC0R , W4T0T4 » W4TE3T 
$. W1 »W1K«DH2M1 «STTO*DELTO .DELTOD* GAM A 8 1 GAM A7 « VO » V7« VB,GAM7P1 .GAM7M1 
$t0NSRT«SFC»SFCC*0F4 *OF6 iFOF6'F4iF 5A»FAB.F6»FNOW1 . FN» F6A » F0F4 » F05A » 
$Wa»RTS.RTA«RT6,RT5A<.CP7*CP8»CAPA8.CAPA7*0W8»WCB. I COMB » I PP » I 3 « ICOUN 
•6T* II<ONT« JCONT* IC,KTEST»THET1 »THETO»OWl DW2 ^ DDHP ♦ 7 1 ♦ Z2 ♦ Z3 « Z« ETAC » ET 
SA0»FN0W1D«HT6»PH1 H6»T5A* T4 »RAT*VJC0W1 C«PBTAB»HDT»DTTA8* IPS,DTEMP »ET 
$AP»ALTSAV»EMSAV»N0SAV»FNSAV»W2SAV»SFCSAV»N0* IN0ZZ«TT2T1D« HT5A »PHTl 
$» TT4TT2 »PHIT5»PSP 15 ^ PH I T5M»HT1 » PT2T 1 L *PS , TABA7 » T8VAL » T7VAL » NL I M « NU 
SMBER*NSTART«NFINAL« AVAL«XMNSLU»TABN«ETA4 9ETA5*TABPH5»TA8H4»HPEXT» 1 
SCASE* IPC«PT4T2D9HGijT5VAL9WC 
COMMON /NEW/ FD^FG^WF^DLDBiDFGS^FGl »FN1 »SFC1 
COMMON /PPP/ I PUNCH 
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WG=W1*( 1 «+F6 ) 

IF (ZTHRUST.GT.OoO ) GO TO 2 
IF ( IPUNCH*GT.O ) GO TO 3 

WRITE <6»7) FNtSFCtFNOWl . PT 1 PTO ♦ PT2PT 1 »ETA2 » W1 K ♦ ONSRT . TTO « TT2 « TT4 < 
1PT4PT2»TT5»PT5PT4«TT5A *TT6«PT6T5A.F6»T8»PTeP8» V0tV8«EM8»EN»A4A40»A 
27A7D«A80A7 i A 8»CAPA8»W1 » ETAP » ETAC » ETAO » FNOW 1 D t WF » WG » PT7P7 » PTl PT7 »PT 
35T5A‘«PT6T5A»PT8PT6tTTl »WB»WBOWl » DLDB « DFGS tFC 1 »FN1 ♦Spci *DTEMP«WC»F4 
IF ( IPUNCH.EO.O ) GO TO 6 
IF ( IPUNCH.EQ*! ) GO TO 4 
ZFNINZCNT )=FNI 
YFN(NZCNT« I ) =WG 
YFN(NZCNT,2 )=V8 
YFN < NZCNT , 3 ) =CAPA8 
YFN(NZCNT,4 )=FD 
YFN(NZCNT,5 }=FN 
YFN(NZCNT«6)=TT4 
YFN ( NZCNT « 7 ) =DLDB 
YFN ( NZCNT . 8 ) =DTEMP 
YFN(NZCNT,9)=SFC 
YFN(NZCNT«10)=T8 
NZCNT=NZCNT+1 
RETURN 

IF { IPUNCH.EO.l ) GO TO I 
IF ( IPUNCH*EQ«2 > GO TO 5 
IF (IPUNCH.EQ. 3) GO TO 5 

PUNCH 8« EMACH.ALT.P0WER»FGI .FD.WF. VS.CAPABtWG 
RETURN 

WRITE '(6»9 ) FN.FG.FD.SFC.WF. W1 .WG.FN1.FG1 .EN.SFCl .WlK.VB.CAPA8.TT4 
1 , A4A4D.TT1 

IF ( IPUNCH.EQ. 3 ) GO TO 4 

CONTINUE 

RETURN 

FORMAT ( iH07X2HFNFl2.5.6X3HSFCFI2.5»4X5HFN/WlFl2.S.2X7HPTl /PTOF12. 
15.2X7HPT2/PT1F12.5.5X4HETA2F12.5/6X3HW1KF12.5. 4X5HN/SRTF12.5.6X3HT 
2T0F12.5.6X3HTT2F12.5 .6X3HTT4F12.5.2X7HPT4/PT2F12.5/6X3HTT5F12.5.2X 
37HPT5/PT4F12.5.5X4HTT5AF12.5.6X3HTT6F12.5. 1 X8HP [ 6/P T 5AF 1 2 • 5 . 7X2HF6 
4F12.5/7X2HT8F12.5.3X6HPT8/P8F12.5.7X2HVOF12.5.7X2HV8F12.5.7X2HM8F1 
52.5.8X1HNF12.5/3X6HA4/A4DF12.5.3X6HA7/A7DF12 .5.4X5HA8/A7F12.5.7X2H 
6A8F12.5.4X5HCAPAaFi2.5»7X2HWlFI2o5/5X4HETAPF12.5.5X4HETACF12.5. 5X4 
7HETA0F12o5.3X6HFN/Wl DFI2.5.7X2HWFFI2.5.7X2HWGF12.5/4X5HPT7P7F12 .5. 
83X6HPT1PT7F1 2.5 . 3X6HPT5T5AF 1 2 « 5 . 3X6HOT6T5AF 1 2 * 5 ■> 3X6HPTBPT6F 1 2 o5 , 5X 
93HTT 1 F 1 2 o 5/7X2HW3F 1 2 .5 . 4X5HWB0W 1 F 1 2 . 5 . 5X4HDLD3F 12.5. 5X4HDFGSF 12.5. 
S6X3HFG1 F12o5 .6X3HFN1 F12.5X5X4HSFC1F12.5.4X5HDTEMPF 12 .5. 7X2HWCF1 2.5 
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00 O' 


$»7X2HF4F12o5) 

FORMAT (F5«2»F10.1 iF5.0*4F10ol «F10.3»F10.1 ) 

FORMAT < 1H03X2HFNF1 1 .4»5X2HFGF1 1 .AoeXaUPDFl 1 .4 »6X3HSFCF 1 0 . 5 -.7X2HWF 
1F10.4»7X2HW1F10«4»7X2HWGF10.4/3X3HFN1F1 1 .4»4X3HFG1F1 i ,4,6X1 HNFl 1 .5 
2 ,5X4HSFC1 FI 0.5 ,6X3HW IKFl 0»4 , 7X2HV8F 1 0 • 4 , 4 X5HC APA8F 1 0 . 4/3X3HTT4F U • . 
35, 1X6HA4/A4DF1 1 • 5 «SX3HTT 1 F 1 1 .5) 

END 
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SUBROUTINE TURB 

DIMENSION AlTSAV( 15) , EMSAV(15), NOSAVdS), FNSAV(15), 
ISFCSAVI 15) 

DIMENSION TABTEM(78), TABH{78), TABPHI{78>» 


1TABPT1(19>, TABNST(9), TABPT2(09), TABZ { 1 I ) , 
2i ANS(5), SFCC(25), TABPSP ( 78 ) 

DIMENSION PBTABC8), HDT(20), DTTA8(20) 
dimension TASPH5(78>, TABH4 ( 78 ) 

DIMENSION TA8N<12), TABA7{12)» RAT(2), T4(3) 
DIMENSION TABVJ(17), TA8DLDB<I7) 


TABPSH (78 ) 
TABWl (99 ) 


W2SAV ( 15 ) ♦ 

, TABM< 19 ) , ■ 
, TABETA(99) 


I NPUT 

COMMON /BKl/ ALT,EMACH,PT2T1D,TT4N, WCOWl , WB » A7A7DN , PT65AD , A4A4DN , P 
180P8, W1 D,H,ETA2D,ENN,X, V, TF,PT8PT6,CV,TT6MAX,PT6T5A ,ETA6« IDES,VAL , 
2GAMO,RO,AJ,TSLS*PSLS,G,TABTEM,TABH,TABPHI , TABPSH, TABPSP,T ABM, TABNS 
3T, TABPT2 , TABZ, TABWl , TABETA , T ABPT 1 , I DDBS ♦ T ABV J , TABDLDB 


computed values 


COMMON /BK2/ WbOWI , AOAS , ABCOR , A 80A7 , A7A7D , ANS , APO , A4A4D » A7 , A8 , AO , A 
laoWl ,ABEM7,TO,T8,T7,T7I , T8I ,TT0,TT2,TT1 ,TT2I ,TT2TT1 » TT5 , TT5A , TT7TT 
2 1 , TT6 , TT7T T6 , TT4 , PO , P80P0 , PT7 , HF ♦ H8M , H8 , H8 I , H8 I M , H7M , H7 , H7M I , HA3 , H 
3AB5,H0,HT0,HT2,HT2I , HT4M , HT4 , HT5M ♦ HTS » HT505 , HT5AM , HT5H , HT7M , HT7M I , 
4HT6M , HT8M , HT5A05 , PH I O , PH I TO , PH I T2 I , PH I T I , PH I T4 M , PH I T4 , PH I T5A , PHT5A 
5M, PH I 8M, PH IT8M, PH 1 8,PH1 2 I , PH I T6M , PH 1 7M , PH I T6 , PH I T7M • PH0F6 , PH0H6 , PT 
60P0»PT1PT0,PT2T1S,PT2PT1 ,PT4PT2 , PT5PT4 , P T5T5A , P I 8P8,Pr 1 PI 7,P I 7P7,P 
7T6PT5«PT6T5D,PT5AM,PT55AD»PSIH4,PSPI4,PS1H5,PSP15A,PSIH5A ,PSPI8»PS 

8 I H8 , PS I H7 , PSP I 6 »PS I PH 1 ,PS I H& , PS 1 H8 1 , PSP \1 , EN , ETA2 , ET5A , ETETS , ET A2S 

9 ,ET2D<>EM7 ,EM8, W1 cor , W1 DCOR , U/4C0R , W40C0R , W7C0R , W7DC0R , W4T0T4 , W4TEST 
S»W1 »W1!<,DH2M1 , STTO,DELTO,DELTOD, GAMAS, GAMA7, VO, V7 ,V8,GAM7P1 ,GAM7M1 
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$«ONSRT»SFC»SFCC»OF4»OF6»FOF6»F4*F5A«FAB*F6tFNOWl ♦ FN » F6A , F0F4 « F05A » H 32 
$WBtRT5»RTA»RT6»RT5A^CP7»CP8f CAPAa«CAPA7.0WB»WCB. ICOMB. I PP t IS» ICOUN H 33 
$T« IKONT* JCONT» IC«KTEST iTHETl »THET0»DW1 t d'W 2 ♦ DDHP ♦ Z 1 ♦ Z2 , Z3 1 Z » ETAC » ET H 34 
SAOfFNOWl D*HT6»PHlH6i T5A* T4 » RAT ♦ WCOW 1 C « PBTAB * HOT » DTTAB • I PS .DTEMP »ET H 35 
£AP»AUTSAV»EMSAV»N0 SAViFNSAV»W 2SAV»SFCSAV«N0« INOZZ*TT2TlD«HT5AtPHTl H 36 
$»TT4TT2tPHIT5»PSPI5»PHIT5M«HTl » PT2T 1 L ♦ PS » TABA7 » T8VAL » T7VAL » NL I M »NU H 37 
$MBER » NSTART » NF I NAL < A VAL i XMNSLU» TABN « ETA4 » ETA5 « TABPH5 » TABH4 » HPEXT » I H 38 


SCASEf IPC»PT4T2D*HS, T5VAL H 39 

COMMON /RRR/ AMP28 ♦ FS28MP » AMP H 40 

F0F4=F4/0F4 H 41 

HT4M=HT4+PS IH4*F0F4 H 42 

CALL DI SCOT I TT4 f TT4 « TABTEM , TABPH I , TABPHI ♦- 1 1 « 78 » 0 ♦ PH I T4 ) H 43 

CALL D I SCOT ( TT4 ♦ TT4 ♦ TABTEM ♦ TABPSP , TABPSP . - I 1 . 78 . 0 » PSP I 4 ) H 44 

PHIT4M=PHIT4+PSPI4*F0F4 H 45 

DDHP=*7068*HPEXT/W1 H 46 

DH2M1=DH2M1 +DDHP H 47 

HT5M=HT4M-DH2MI/(0F4*WCB) H 48 

IR0U=5 H 49 

CALL ITT (HT5M»TABH4 »TABH»TABTEM»TABPSH»T5VAL»F0F4 < IR0U*TT5»HT5«PS H 50 

11H5) H 51 

CALL. D I SCOT < TT5 » TT5 ♦ T ABTEM . TABPH I , TABPH I » - 1 1 < 78 » 0 » PH I T5 ) H 52 

CALL DISCOT ( TT5 . TT5 » TABTEM , TABPSP » TABPSP . - 11 ♦ 78 ♦ 0 « PSP I 5 > H 53 

PHIT5M=PHIT5+PSPI5*F0F4 H 54 

RT5=1545**F0F4» (FS28MP+AMP/F4)/AMP28 H 55 

PT5PT4 = EXP( ( (PHIT5M-PHIT4M)*AJ)/(ETA5*RT5n H 56 

RETURN H 57 

END ' H 58- 


SUBROUTINE COND I I 

DIMENSION ALTSAV(15)« EMSAV(15)» NOSAVdSli FNSAv'<15), W2SAV(15)* I 2 

iSFCSAVdS) I 3 

DIMENSION TABTEM(78)» TABH(78)» TABPHI(78)* TABPSH(78)» TABM(19)» I 4 

1TABPTI(19). TABNST<9). TABPT2(09), TABZdl)* TABW1(99), TABETA(99) I 5 
2» ANS(5). SFCC(25 )i TABPSP ( 78 ) I 6 

DIMENSION P8TAB<8)* HDT(20)t DTTAB(20) I 7 

DIMENSION TABPH5(78)« TABH4(78) I 8 

DIMENSION TAeN(12)» TABA7(12)« RAT(2)» T4(3) I 9 

DIMENSION TABVJd7)» TABDLDB(17) I 10 

I 1 1 

INPUT 1-12 
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COMMON /BKl/ ALT,EMACH«PT2T1D«TT4N» WCOWl t WB « A7 A7DN » PT65 AD « A4 A4DN « P 
180P8.W1D»H»ETA2D«ENN.X» YtTF»PT8PT6tCViTT6MAX«PT6T5A9ETA6» IDES«VAL» 
2GAM0»R0« AJ«TSLS9PSLS9G»TABTEM*TABH»TABPHI «TA8PSH» TABPSP»TABM » TABNS 
3T*TABPT29TABZ»TABW1 9 TABETA ♦ TABPT 1 9 IDDBS 9 TABVJ 9 TABDLDB 

COMPUTED VALUES 

COMMON /BK2/ WBOWl 9 AOAS 9 ABCOR 9 A80A7 9 A7A7D 9 ANS 9 APO 9 A4A4D 9 A7 9 AS 9 AO 9 A 
180W1 9 ABEM 79 TO 9 T 89 T 79 T 7 I 9 T 8 I 9 TTO 9 TT 29 TTI 9 TT 2 I 9TT2TT1 9TT59TT5A 9TT7TT 
21 9TT69TT7TT69TT49P09P80P09PT79HF9HaM9H89H8I 9 H8 I M 9 H7M 9 H7 9 H7M I 9 HAB 9 H 
3 AB 59 HO 9 HTO 9 HT 29 HT 2 I 9HT4M9HT49HT5M9HT59HT5059HT5AM9HT5H9HT7M9HT7MI 9 
4 HT 6 M 9 HT 8 M 9 HT 5 AO 59 PH IO 9 PHITO 9 PHI T21 9 PHITI 9 PHI T4M9PHI T4 9 PHI T5A9PHT5A 
5 M 9 PH I 8M 9PH I T8M 9 PH 18 9PH I 2 1 9 PH I T6M 9 PH I 7M 9 PH I T6 9 PH I T 7M 9 PH0F6 9 PHOH 6 9 PT 
60PO9PT1PT0 9PT2T1S9PT2PT1 , PT4PT2 9 PT5PT4 , PT5T5 A 9 PT8P8 9 PT 1 PT7 9 PT7P7 9 P 
7T6PT5 9PT6T5D 9PT5AM9 PT 55 AD 9 PS IH49 PSPI4 9PS 1H5 9 PSPI 5A 9PSIH5A 9 PSP 18 9PS 

8 I H8 9 PS I H 79 PSPI 6 9 PSI PHI 9 PSI H6 9 PSIH 8 I 9PSPI79EN9ETA29E15A 9E TETS 9ETA2S 

9 9ET2D9EM7 9 EM 89 Wi COR 9 Wl DCOR 9 W4C0R 9 W4DC0R 9 W7COR 9 W /■UCOR 9 w4 t O l4 9 W4 I Eb I 
$ 9 W 1 9 WIK 9DH2M1 9 STTO 9 DELTO 9 DELTOC 9 GAMA 89 GAMA 79 VO 9 V 79 V 89 GAM 7 PI 9GAM7M1 
$9 0NSRT9SFC9SFCC9 0F4 9 OF 69 FOF 69 F 49 F 5 A 9 FAB 9 F 69 FNOWI 9 FN 9 F 6 A 9 FOF 4 9 F 05 A 9 
£W89RT59RTA9RT69RT5A9CP79CP89CAPA89CAPA790WB9 WCB 9 I COMB, I PP 9 I S 9 ICOUN 
£T9 IKONT 9 JCONT 9 IC 9 KTEST 9 THETI 9THET09DW1 9DW29DDHP9Z1 9 Z2 9 Z3 9 Z 9 ETAC 9 ET 
$A0 9FNOW1D9HT6 9PHIH6 9 T 5 A 9 T 49 RAT 9 WCOWI C 9 PBTAB 9 HDT 9 DTTAB 9 IPS 9 DTEMP 9 ET 
£AP9ALTSAV9EMSAV9NOSAV9FNSAV9W2SAV9SFCSAV9N0 9 INOZZ 9 TT2T1D9 HI5A 9 PHTI 
$,TT4TT2 9PHIT5 9PSPI59PHIT5M9HT1 9 PT 2 TIL 9 PS 9 T AB A7 9 T8 V AL 9 T7 VAL 9 NL I M 9 NU 
$MBER9 NSTART 9 NF I NAL 9 A VAL 9 XMNSLU 9 TABN 9 ETA4 9 ETA5 9 TABPH5 9 TABH4 9 HPEXT 9 I 
SCASE9 IPC9PT4T2D9HS9T5VAL 

. Wl D=W1D*PT1PT0 
W1=W10 

PT2T1S=Y*PT2T1D 
ONSRT= 1 * 

Z=( <PT2T1D-1 • )/<PT2TlS-l , ) -X ) * 1 9 / ( 1 , -X ) 

CALL DISCOT ( ONSRT 9 Z 9 T ABNST 9 TABW 1 9 TABZ 9 1 1 9p9 9 1 1 9 W 1 DCOR ) 

CALL DISCOT ( ONSRT 9 Z 9 TABNST 9 TABET A 9 T ABZ 9 1 1 9 99 9 1 1 VETETS ) 

ETA2S=ETA2D/ETETS 

ETA2=ETETS*ETA2S 

PHTl=PHITO 

PHIT2 I =PHT1 +RO/A J*ALOG (PT2T1 D) 

call DISCOT (PH1T2I 9PHIT2I »TASPHI 9 TABTEM 9 TABTEM 9 -I 1 9789 O 9 TT 2 I ) 

CALL DISCOT < TT2 1 9 TT2 I 9 TABTeM 9 T ABH 9 T ABH 9 - 1 1 9 78 9 0 , HT2 I 3 
DH2M1 = <MT2! -HTl 3/ETA2 
HT2=HT1 -f-DH2Ml 

CALL DI SCOT ( HT2 9 HT2 9 TABH 9 TABTEM 9 TABTEM 9 - 1 1 9 73 9 0 9 T T2 ) 
TT2T1D=TT2/TT1 


I 1 3 
I 1 4 
I 15 
I 16 
I 17 
I 18 
I 19 
I 20 
I' 21 . 

I 22 
I 23 
I 24 
■I 25 •. 
I 26 
I 27 
I 28 
I 29 
I 30 
I 31 
I 32 
I 33 
I 34 
I 35 
I 36 
1 37 

1 38 
I 39 
I 40 
I 41 
I 42 . 

I 43 
I 44 
I 45 
I 46 
I 47 
I 48 
I 49 
I 50 
I 51 
I 52 
I 53 
I 54 
I 55 
I 36 


ORIGINAL PAGE IS 
OP POOR QUALITY 



TT2TT1=TT2T1 D 

WB0W1=WB/W1 

0WB=1 .-WBOWl 

WC0W1C=WC0W1 

WCB= < 1 .-WCOWl -WBOWl ) 

CALL COMB 

RETURN 

END 


I 57 
.1 ' 58 
■ I 59 
.1 60 
I 61 
I 62 
1 63 

I 64- 
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SUBROUTINE FRES J 1' 

DIMENSION ALTSaVUB)* EMSAV(15)« N0SAV(15)» FNSAV<15)t W2SAV(15)* J 2' 
ISFCSAVdS) J 3 

dimension TABTEM(78>» TABH { 78 ) » TABPHI(78>, TABPSH(78)« TABM<19)c J 4 
1TABPT1U9)* TABNST(9)f TABPT2<09), TABZ ( 11 ) « TABW1C99)» TABETA(99) J 5 
2« ANS<5). SFCC(25)» TABPSP(78) J 6 

DIMENSION PBTAB(8)» HDT(20)» DTTAB(20) J 7 

DIMENSION TABPH5<78), TABH4(78) J 8 

DIMENSION TABN(12)» TABA7(12)» RAT (2)* T4 ( 3 1 J. 9 

DIMENSION TABVJ<17)t TABDLDB(17) J 10 

C • J 1 1 

C INPUT J 12 

C ■ J 13 

COMMON /BKI/ ALTvEMACH»PT2Tl D*TT4N« WCOWl » WB» A7A7DN»PT65AD tAAAADN.P J 14 
180P8»WlD»H»ETA2D»ENN»X»Y»TF»PT8PT6»CVtTT6MAX»PT6T5A»ETA6» IDES-.VAL* J 15 
2GAMO»RO«AJ»TSLS»PSLS,G»TABTEM»TABH*TABPHI ♦TABPSHdABPSP.TABMfTABNS J 16 
3T»TABPT2»TABZ»TABW1 »TABETA*TABPT1 , IDD8S«TARVJ«TABDLD8 J 17 

' C , • J 16 

C COMPUTED values J 19 

C J 20’ 


COMMON /BK2/ WBOWl ♦ AOAS » ABCOR * A80A7 ♦ A7A7D » ANS « APO « A4A4D. A 7 » AS ♦ AO » A J 21 
laoWl *ABEM7*T0»T8,T7»T7I »T8I fTT0«TT2»TTl «TT2I »TT2TTl » TT5 , TT5A , TT7TT J 22 
21 «TT6«TT7TT6»TT4»PO»P8OP0.PT7iHF,H 8M.H8.H8I « H8 I M i H7M ♦ H7 , H7M I «HAB,H J 23 
3AB5»H0.HT0»HT2»HT2I « HT4M « HT4 » HT5M » HT5 ♦ HT505 » HT5AM » HT5H » HT7M . HT7M I . J 24 
4HT6M»HT8M«HT5A05»PHI0»PHIT0»PHIT2I iPHITl »PH I T4M« PH I T4 , PH I T5A »PHT5A J 25 
5M»PHI8M»PHIT8M»PHI8«PHI2I « PH I T6M .PH I 7M» PH I T6 9 PH I T7M» PH0F6 *PH0H6 »PT J 26 
60P0.PT1 PT0.PT2T1 5 .PT2PT1 .PT4PT2 . PT5PT4 . PT5T5A .PT8P8 »PT 1 PT7.PT7P7. P J 27 
7T6PT5.PT6T5D.PT5AM.PT55AD.PS I H4» PSP 14 «PS IH5.PSP1 5A . PS IH5A .PSP 18 ,PS J 28 
8 I H8.PSIH7.PSP16.PS I PH I . PS I H6 ♦ PS I H8 I . PSP I 7 . EN . ETA2 . ET5A . ETETS . ET A2S J 29 
9 . ET2D » EM7 »EM8 » W 1 COR . Wl OCOR . W4C0R ♦ W4DC0R . W7C0R . W7DC0R , W4T0T4 . W4TEST J 30 
S.Wl .W1K.DH2M1 .STTO. DELTO.DELTOD. GAMAS. GAMA7. VO »V7. V8.GAM7P1 .GAM7MI J 31 



u o 


1 = 

1 

rv> 

an 


C 

C 

C 

1 


C 


$ « ONSRT ♦ SF C«SFCC«0F49 0F6 «F0F6»F49F5A» FAd » r 6 « FNO W 1 « FN o F6A « F OF 4 * F05A ♦ 
$W8»RT59RTA«RT69RT5A9CP7 cCP8»CAPA8»CAPA7«0WB9 WC8» ICOMB» I PP ? IS» ICOUN 
$T« IKONT* JCONT* IC»KTEST9THET1 »THET0» D'V 1 . DW2 « DDHP» 2 1 » Z2 ♦ Z3 9 2 9 ETAC ■>E7 
$A0»FN0W1D*HT6»PHIH6 • T5A • T4 « R AT » WCOW 1 C « PBTA8 « HOT « DT TAB i IPS tDTEMP «ET 
$AP»ALTSAV«EMSAV9N0SAV«FNSAV9W23AV 9SFCSA\/9N0« I N022 9 TT2T1D9HT5A 9 PHTI 
S,TT4TT2 9PHIT5*PSPI5 9PHIT5IV)9HT1 9 PT 2 TIL 9 PS 9 TAdA7 9T8VAL9T7VAL9N’LIM»NU 
SMBER 9 NSTART 9 NF I NAL 9 AVAL 9 XMNSLU 9 TABN 9 ETA4 9 ETA5 tTABPHB 9 TABH4 9 HPEXT 9 I 
$CASE9 1PC9PT4T2D9HS9T5VAL 
HOT( 1 )=-l a 

call AT62SP (ALT9ANS9DTEMP9PBTAB9HDT9DTTAB) 

P0=ANS ( 2 ) 

AO=ANS ( 4 ) 

T0=ANSI3)*1 «8 
TT4=TT4N 

STANDARD DAY 

call DISCOT (T09T09TABTEM9TABH9TABH9-I1 97 B 9 O 9 HO) 

HTO=HO+ (GAM0-»-EMACH#*2»T0*R0 )/(2a*AJ ) 

HTl =HTO 

CALL DISCOT (HTO 9 HTO 9 TABH 9 TABTEM 9 TABTEM 9 -I 1 9789 O 9 TTO) 
call DISCOT (TO 9 TO 9 TABTEM 9 TABPHI .TABPHl 9 - 1 1 9789 O 9 PHIO) 

CALL DISCOT < TTO 9 TTO 9 TABTEM 9 TABPH 1 9 TABPH I 9 - I 1 9 78 9 0 ♦ PH I TO ) 
P70P0=EXP( AJ/RO* (PHI TO-PHIO ) ) 

STTO=SQRt < TTO/TO*TO/TSLS ) 

THET0=STT0**2 

TTl=TTO 

PHTl=PHITO 

DELTO=PTOPO*PO/PSLS 

INLET CONDITIONS 

CALL DISCOT ( EM ACH 9 EMACH 9 TABM 9 TABPT 1 9 TASPT 1 9 - 1 1 9 I 9 9 0 9 PT 1 PTO ) 

return 

END 
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J 33 
J 34 
J 35 
J 36 
J 37 
J 38 
J 39 
J 40 
J 41 
J 42 
J 43 
J 44 
J 45 
J 46 
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J 48 
J 49 
J 5 O 
J 51 
J 52 
J 53 
J 54 
J 55 
J 56 
J 57 
j 58 
J 59 
J 60 
J 61 
J 62 
J 63 
J 64 
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SUBROUTINE COMB 

DIMENSION TABTEM(78)9 TABH(78)9 TABPHl ( 78)9 TABPSH<78)9 
1TABPT1(19)9 TA3NST(9)» TABPT2(09)9 TABZ(11)9 TASW1<99}9 
29 ANS<5)9 SFCC(25>9 TASPSP ( 78 ) 


TABM <19)0 
TASETA(99> 


l< 

l< 

l< 

K 


1 

2 

3 

4 



uro*-- non non 


DIMENSION PBTAB(8)« HDT(20)» DTTAB(20) K 5 

DIMENSION ALTSAVUS)^ EMSAVdS)? N0SAVU5)* FNSAV(IS), W2SAV(15)* K 6 

1SFCSAVU5) K 7 

DIMENSION TABPH5(78)* TABH4(78) K ’ 8 

DIMENSION TABN(12)» TABA7<12)» RAT(2)« T4<3) K 9 

DIMENSION TABVJ(17)« TABDLDB(17) -K. 10 

K 1 1 

INPUT K 12 

K 13 


COMMON /BKl/ ALT«EMACH«PT2T1 D«TT4N«WC0W1 »WB* A7A7 DNiPT65AD oA 4A4DN»P K 14 
180P8»W1D»H*ETA2D,ENN»X«Y«TF.PT8PT6»CV»TT6MAX<PT6T5A«ETA6* IDES»VAL« K 15 
2GAM0fR0«AJ? TSLS«PSLStG« TABTEM»TABH»TABPHI « TABPSH » TABPSP , TABM « TABNS K 1 6 
3T.TABPT2«TABZ»TABW1 *TABETA»TABPT1 dDDBS«TABVJ«TABDLD3 -K 17 

•K 18 '■ 

COMPUTED VALUES K 19 

•K 20 

COMMON /BK2/ WBOW 1 ♦ A0AS» ABC0R» A80A7 « A7A7D » ANS » APO » A4 A4D» A7 ♦ A8 , A O *A K 21 
180W1 »ABEM7»T0«T8,T7»T7I » T8 I ♦TT0»TT2.TTl .TT2I »TT2TT1 « TT5 . TT5A ♦ TT7TT K 22 
21 »TT6«TT7TT6»TT4«P0«P80P0»PT7*HF»H8M»H8«H8I « H8 I M ♦ H7M ♦ H7 * H7M I ,HAB»H K 23 
3AB5*H0*HT0*HT2iHT2I »HT4M«HT4»HT5M.HT5«HT505»HT5AH»HT5H»HT7M»HT7MI * K 24 
4HT6M*Hf8M»HT5A05»PH10»PHlTO*PHIT2I ,PHIT1 . PH I T4M » PH I T4 , PH I T5A » PHT5A K 25 
5M,PHiaM»PHIT8M*PHI8«PHI2I ♦ PH I T5M « PH I 7M » PH I T6 * PH I T7M »PH0F6 »PH0H6 »PT K 26 
60PO»PT1PTO«PT2T1S*PT2PT1 , PT4PT2 « PT5PT4 * PT5T5 A ♦ PT8P8 « PT 1 PT 7 » PT7P7 ♦ P K 27 
7T6PT5»PT6T5D»PT5AM»PT55AD.PSlH4»PSPI4fPSIH5«PSP15A«PSIH5A»PSPI8tPS K -28 . 
8 I H8 « PS IH7 t PSP I 6 *PS I PH I ,PSI H6 «PS I H8I . PSP I 7 «EN ♦ ETA2 * ET5A , ETETS » ETA2S K 29 
9 ♦ ET2D * EM7 « EM8 » W i COR ♦ W1 DCOR * W4C0R ♦ W4 DC OR « W7C0R » W7DC0R i W4T0T4 t W4TEST K 30 
S»W1 «WlKtDH2Ml » STTO » DELTO, DELT0D» GAMAS »GAMA7t VO iV 7»V8»GAM7P1 »GAM7M1 K 31 " 
$»ONSRT»SFC»SFCC«OF4«OF6»FOF6tF4»F5A.FAB»F6*FNOWl *FN »F6A , F0F4 *F05A ♦ K 32 
$W8«RT5»RTA*RT6*RT5A»CP7»CP8»CAPA8.CAPA7*0WB* WCBf ICCMB» 1PP»IS^ ICOUN K 33 
ST» IKONT* JCONT« IC«KTEST«THET1 »THET0»DW1 ♦DW2»DDHPiZ 1 *Z2«Z3'Z*ETAC»ET K 34 , 
$AO»FNOW1D»HT6«PHIH6»T5A»T4»RAT»WCOW1C*PBTAB*HDT»DTTAB* IPStDTEMP»ET K 35 
$AP»ALTSAVtEMSAV»NOSAV»FNSAV»W2SAV,SFCSAV» NO t I NOZZ * TT2T 1 D » HT5A «PHTl K 36 
$,TT4TT2»PHIT5«PSP 15.PHIT5M,HT1 . PT2T 1 L ♦ PS » TAB A7 , T8V AL t T7VAL * NL I M , NU K 37 
$MBER»NST/^RT»NFINAL»AVAL.XMNSLU«TABN*ETA4»ETA5 iTABPH5»TABH4»HPEXT, I K 38 . 


SCASE* IPC»PT4T2D*HS, TSVALtWC K 39 

IF (TT4€LT.TT2) go to 2 K 40 

IF (TT4oGT«TT4N ) GO TO 1 K, 41 

GO TO 3 K 42 

TT4=TT4N K 43 

GO TO 3 , K 44 

TT4 = TT2+50o l< 45 

IF {TT4oLEo2160» 5 GO TO 8 K 46 

WCOWlC = WC0WdMTT4-2l60 9 ) K 47 

WC=WCOWlCi(-Wl l< 47A 



p'CO ou. u ouo 


4 WCB= < I •-WC0W1C-WB0W1 ) K 

HF=o5*(TF-536» ) K 49 

•CALL DISCOT ( TT4 » TT4 < T ABTEM » TABH , T ABH . - 1 1 , 73 » 0 » HT4 ) K 50 

call di scot (TT4» TT4 »TABTEM» TABPSH«TABPSH«-1 1 * 78*0 »PSIH4 ) K 5l 

F4= (HT4-HT2 )/(ETA4* ( H-HS I-HT4-PS IH4+HS-fHF ) K 52 

0F4=l*-fF4 •< 53 

IF (IDES-1) 5,6»6 K 54 

5 W4DC0R=W1DC0R*SQRT(TT4/TT2*TT2T1 D)*l •/(PT2T1D*PT4T2D)*0F4*WCB K 5S 

PT4PT2 = PT4T2D , K •• 56 ' ' 

GO TO 7 K 57 

6 PT4PT2=1*-( (1 •-PT4T2D)*W1C0R**2*W1DC0R**C-2)*(PT2T10/'PT2PT1 )**2*TT K 58 

12TT1/TT2T1D) K 59 

IF (PT4PT2.lt. .7) PT4PT2=.7 K 60 

W4C0R=W1C0R»SQRT(TT4/TT2*TT2TT1 >*1 ./(PT2PT1*PT4PT2)*0F4*WCB K 61 

RETURN « 62 

WC0WIC=0 K 63 

GO TO 4 K 64 

END K 65- 


subroutine PARP L ' 1 

DIMENSION ALTSAVdS)* EMSAV(15>. N0SAV(15)» FNSAV<15). W2SAV(1S)« L 2 
ISFCSAV(I5) L 3 

DIMENSION TA8TEM(7a). TABH ( 78 ) » ‘ TABPH I ( 78 ) » TABPSH{78). TABM(I9)« L 4 
1TABPT1(I9)» TA8NST(9). TABPT2(09 m TABZ(11)» TABWM99). TABETA(99) L 5 
2. ANS(5>« SFCC(25)« TABPSP(78) L 6 

DIMENSION PBTAB<8)« HDT(20)» DTTAB(20) L 7 

DIMENSION TABPH5(78)* TABH4<78) L 8 

DIMENSION TABN(12)» TABA7U2)’ RAT(2)» T4(3) L 9 

DIMENSION TABVJ(17)» TABDLDB(17) L lO 

L 1 I 

I NPUT t- 1 2 

L 13 

COMMON /8K1/ ALT,EMACH»PT2T1D«TT4N»WC0W1 « W8.A7A7DN«PT65AD»A4A4DN«P L 14 
ie0P8.WlD»H.ETA2D»ENN.X.Y«TFiPT8PT6»CV^TT6MAX«PT6T5A»ETA6» IDEStVAL. L 15 
2GAM0»R0<AJ»TSLS»PSLS,G«TABTEM«TABH»TA8PHI 9TABPSH«TABPSP,TA8M«TABNS L 16 
3T«TABPT2»TABZ»TABW1 «TABETA*TABPT1 « IDDBS*TABVJ*TABDLDB L 17 

L 18 

COMPUTED values L 19 

L 20 

common /BX2/ WbOWI ? A0AS<> ABC0R9A30A7-SA7A7D9 ANS?AP0? A4A4D?A7oA3 9AO?A L 21 



ro o n 


> 

ro 

00 


1 

c 


3 

4 

5 


6 

7 


180W1 »ABEM7tT0*T8*T7«T71 »T8I ,TT0»TT2»TT1 .TT21 »TT2TT1 « TT5-.TT5A ♦TT7TT 
21 »TT6 9TT7TT6»TT4«P0»P80P0»PT7»HF 9H8M«H8»H8I » H8 I M » H7M » H7 » H7M I »HAe»H 
3AB5*H0tHT0«HT2»HT2I « HT4M »HT4 » HT5M» HT5 ♦ HT505 » HT5AM « HT5H« HT7M * HT7MI>, 
4HT6M ♦ HT8M « HT 5A05 » PH I 0 * PH I TO « PH I T2 I » PH I T 1 . PH I T4 M « PH I T4 » PH I T5A t PHT5A 
5MfPHI8M«PHlT8M«PHl8»PHI2I *PHIT6M «PH I 7M. PHI T6» PHI T7M,PH0F6»PH0H6»PT 
60PO»PT1PTO»PT2T1S<PT2PT1 ,PT 4PT2 « PT5PT4 . PT5TSA ♦ PT8P8 ♦ PT i PT7 * PT7P7» P 
7T6PT5»PT6T5D*PT5AM»PT55AD«PS IH4«PSPI4»PSIH5«PSPI5A»PSIH5A iPSPie tPS 
8 I H8 • PS I H7 ♦ PSP I 6 »PS I PH I »PS I H6 ♦ PS I H8 1 « PSP I 7 « EN « ETA2 « ET5A * ETETS « ET A2S 
9 ♦ ET2D » EM7 » EM8 » W 1 COR ♦ W 1 DCOR » W4C0R * W4DC0R » W7C0R » W7DC0R « W4T0T4 » W4T EST 
$ » W 1 . W I K » DH2M 1 • STTO « DEL TO * DELTOD » GAMA8 . GAMA7 ♦ VO ♦ V7 » V8 » GAM7P 1 ♦ GAM 7M 1 
$*ONSRTtSFC»SFCC»OF4«OF6tFOF6«F4»F5A»FAB»F6»FNOWl i FN »F6A . FOF4 ♦ F05A ♦ 
$W8,*RT5»RTA»RT6»RT5A»CP7«CP8tCAPA8»CAPA7.0WB«WCB» ICOMBt I PP » I S . ICOUN 
ST» IK0NT»JC0NT» IC^KTeST.THETI »THET0»DW1 »DW2*DDHP»Z1 »22 *23 * Z»ETAC »ET 
$ AO » FNOW 1 D » HT6 » PH I H6 « T5 A ♦ T4 * R AT « WCOW 1 C » PBT AB » HD T * DT TAB * I pS , DTEMP * ET 
$AP * ALTSAV *EMSAV .NOSAV * FNSAV ♦ W2SAV ♦SFCSAV « NO » I NOZZ « TT2T 1 D » HT5A .PHT 1 
$*TT4TT2«PHIT5*PSP 1 5 * PH I T5M » HT 1 » PT2T 1 L * PS , TA3A7 » T8VAL » T7VAL * NL IM * NU 
SMBER *NSTART « NF I NAL * AVAL » XMNSLU* T ABN . ETA4’* ETA5 * TABPH5 t TABH4 * HPEXT * I 
SCASE. IPC»PT4T2D*HS*T5VAL 
ITJ = 0 
IPP=0 
EN=ENN 
1TT4=0 
IM=0 

W4C1=W4C0R' 

IF (A4A4DN-1*) 2*18*18 
GO TO 35 

VARIABLE A4/A4D ROUTINE 

TT4=TT4N*A4A4DN**2 

IF (TT4.lt. (TT2+1 00. > ) GO TO 3 

GO TO 4 

TT4=TT2+100. 

A4A4DN=SQRT ( TT4/TT4N ) 

DELTAT= (TT4N-TT4 )/ ( < 1 . -A4A4DN ) /. 05 ) 

IF (IM) 5.5*6 

TT4=TT4N-DELTAT 

IM=1 

GO TO 7 

TT4=TT4“DELTAT 
T4 ( 1 ) = TT4 
CALL COMB 
W4C2=W4C0R 
A4A4D=W4C2/W4C1 
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22 


P 

-23 


L 

24 


,L- 

25 


L 

26 


L 

27 


L 

28 


L 

29 


L 

30' 


l' 

31' 


L 

32 


L 

33 


L 

34 


L 



•L' 

36 


.L 

37 


L' 

30 


L 

39 


L, 

40 


L 

41 , 


L 

42 


L 

43 


L 

44 


L 

45 


L 

46.. 


L 

47 


L 

48 


L 

49 


L 

50 . 


L 

51 


L 

52 


L 

53 


L 

54 


L 

L 

L 

L 

L 

‘55 , 

56 

57 

58 

59 

If 

L 

L 

L 

60 
6 1 
62 


L 

63 


L 

64 

L 

65 




vO D 


IF UM~1 ) a« a »9 L 66 

IF (A4A4D-(A4A4DN+.025) ) 1 1 »9»9 L 67 

call TURB L 68 

CALL COOL L ■ 69 

PT6T5A=PT65AD L. 70 

HT6M=HT5AM L 71 

PHIT6M=PHT5AM L 72 

F6=F5A ■ L 73 

0F6=le+F6 *- 74 

F0F6=F6/0F6 L. 75 

CALL NOZ2 L 76 

CALL ENG L 77 

EN=ENN L 78 

IF ( IM-1 ) 1 0 • 1 O « 1 1 L 79 

1 0 P0WER=TT4 C 80 

CALL OUTPUN (POWER) L 81 

IF (FNOWl .LT.O . ) GO TO 35 L 82 

IF (A4A4D- ( A4A4DN+.025 ) ) L 83 

11 ABSA4 = ABS ( A4A4D-A4A40N ) L 84 

ITJ=ITJ + 1 L. 85 

IF ( ITJ-1 ) 1 3« ) 3, 12 L .86' 

12 IF (ABSA4-,001) 17,17»13 L 87 

13 IF (lM-2) 14*15«15 L 88 

14 T4 (? )=T4 ( 1 >* (A4A4DN**2 )/< A4A4D**2 ) L 89 

A4C1=A4A4D *- 90’ 

1M=2 , L 91 

TT4=T4(2) E 92 

GO TO 7 , L 93 

15 A4C2=A4A4D L 94 

IF ( A4C1 .EQ»A4C2 ) GO TO 35 L 95 

T4(3)=T4<1 )“( (A4C1**2-A4A4DN^<-*2 ) / ( A4C 1 **2-A4C2**2 ) )* < T4 { 1 ) -T4 (2 ) ) L 96 

T4 { 1 )=T4 (2 ) L 97 

T4(2)”T4(3) t- 98 

TT4=T4(3) t- 99 

A4C1=A4C2 t- 100 

IM=IM+1 L 101 

IF (lM-25) 7»7«16 L 102 

16 WRITE (6*36) L. 103 

GO TO 35 L 104 

17 P0WER=TT4 L l05 

CALL OUTPUN (POWER) L 106 

IF (FNOWl oLToOo ) GO TO 35 L 107 

IF (A7A7DN“1o) 35«l8oia L 108 

C PPCNA L 1 09 



A-30 


18 

TPLUS= 1 .+AVAL 

L 

110 


TMIN=1 9-AVAL 

L 

1 1 1 


A7A7DS=A7A7D 

■ L 

1 12 


T4 < 1 >=TT4-1 OOo 

L 

113 


DO 34 JT=NSTART»NFINAL»NUMBER 

L 

1 14 


TTP=TT4 

. L 

1 15 


1NUM=0 

L 

1 16 


11 = 1 

L 

1 17 


EN=TABN ( JT ) 

L 

1 18 


A7A7DL = TABA7 ( JT ) 

L 

1 1 9 


A7A7DO=A7A7DS*A7A70L 

L 

120 


IF (EMACH-XMNSLU) 20,19»19 

L 

121 

19 

EN=1 . 

L 

•122 

20 

TT4 = T4 ( 1 ) 

L 

123 

21 

CALL MATCH 

L 

124 


IF (I COUNT-2) 22,1,22 

L 

1 25 

22 

IF (ONSRT-.4) 34,34,23 

L 

126 

23 

CALL TURB 

L 

127 


CALL COOL 

L 

128 


PT6T5A=PT65AD 

L 

129 


HT6M=HT5AM 

L 

130 


PH 1T6M=PHT5AM 

L 

1 31 


F6=F5A 

L 

132 


0F6=1 •^-Fe 

L 

133 


F0F6=F6/0F6 

L 

134 


CALL N0Z2 

L 

135 


RAT( I I ) =A7A7D0/A7A7D 

L 

136 


INUM= INUM+1 

L 

137 


IF (RAT { I 1 )-TPLUS ) 24*31,28 

L 

138 

24 

IF <RAT{ 1 n-TMIN) 25,31,31 

L 

1 39 

25 

IF, (INUM-2) 27,30,26 

L 

1 40 

26 

IF .{INUM-25) 30,30,32 

L 

141 

27 

CONTINUE 

L 

142 


T4 (2 )=T4 ( 1 )+ (TTP-T4 ( 1 ) 

L 

1 43 


TT4=T4 (2 ) 

L 

144 


I 1=2 

L 

145 


GO TO 21 

L 

146 

28 

IF (INUM-2) 29,30,26 

L 

147 

29 

T4 <2 )=T4 ( 1 )- (TTP-T4 Cl ) )#.5 

L 

148 


TT4=T4 (2 ) 

L 

149 


11=2 

L 

150 


GO TO 21 

L 

151 

30 

CONTINUE 

L 

152 


IF (RATCi )bEQoRAT(2) ) RAT < 2 ) =RAT ( 1 >-» 1 , 0 1 

L 

153 



T4 (3 ) = ( < 1 o-RAT (2 ) > / <=5AT(2)-RAT ( 1 ) )«-<T4 <2 )-T4 < 1 ) ) ) + T4 (2 ) 
IF (T4(3)oGT»TTP) T4 {3 ) =TTP-25» 

IF (T4(3)oLT.TT2) T4 ( 3 ) =TT2-«-50 « 

TT4 = T4 < 3 ) 

T4 <1 )=T4(2) 

T4 (2 >=T4 (3 ) 

RAT < 1 )=RAT < 2 ) 

I 1=2 

GO TO 21 

31 CALL ENG 
POWER=TT4 

CALL OUTRUN (POWER) 

GO TO 33 

32 WRITE (6.36) 

33 T4 ( 1 >=TT4-1 00* 

IF (FNOWl .LT.O, ) GO TO 35 

34 CONTINUE 

35 A4A4D=1. 

RETURN 

C 

36 format ( 1H06X26HI TTERATI on GREATER THAN 25! 

^ END- 

CO 


L 154 
L 155 
L 156 
L 157 
L 158 
L 159 
L 160 
L 161 
L 162 
L 163 
L 164 
L 165 
L 166 
L 167 
L 168 
L 169 
L 170 
L 171 
L 172 
L 1 73 
L 174 
L 175- 


C 

C 

c 


SUBROUTINE ENG M 

DIMENSION ALTSAV(15). EMSAVdS). NOSAVdS), FNSAV(15), W2SAV(15). M 
ISFCSAVdS) M 

DIMENSION TABTEM(78). TABH(78). TA8PHI(78). TABPSH(78). TABM(19). M 
1TA8PT1(19)» TABNST(9). TABPT2(09). TABZdl). TABW1(99). TABETA(99) M 
2. ANS(5). SFCC(25>« TASPSP(78) M 

DIMENSION PBTAB(8). HDT(20>, DTTAB(20) M 

DIMENSION TABPH5(78). TABH4(78) M 

DIMENSION TABN(12). TABA7(12)« RAT (2). T4(3) M 

DIMENSION TABVJ(17). TA8DLDB(17) M 

M 

INPUT M 

M 

COMMON /BKl/ ALT.EMACH.PT2T1 D.TT4N. WCOWl .WB.A7A7DN.PT65AD.A4A4DN.P M 
180P3. W 1 D.H.ETA2D.ENN.X. Y.TF.PTaPT6eCV.TT6MAX«PT6T5A.ETA6o IDES.VALo M 
2GAM0.RO. AJ.TSLS.PSLS .GiTASTEM.TASH.TABPH! « TA3P3H. T ABPSP » T ASM . TABNS M 
3T.TABPT2.TABZ.TA3W1 . TABETA. TASPT 1 , I DD3S . TASVJ . TA3DL0B M 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 
1 5 
16 
17 


ORIGINAE PAGE iS 
OF POOR QUALITY 



non 


1 

w 

r\5 


M 18 

COMPUTED values M 19 

M 20 

COMMON /BK2/ WBOW 1 o AOAS « ABCOR « A80A7 ♦ A7A7D « ANS « APO ♦ A4A4D 9 A7 » AS . AO 9 A M 21 
180VM1 9ABEM79T09T89T7»T7l9T81 9TTO.TT2»TTl 9 TT 2 I 9TT2TT1 9 TT5 9 TT5A 9 TT7TT M 22 
21 9 TT 6 9 TT7TT6 9 TT4. PO 9 P 8 OPO 9 PT 7 9 HP 9 HBM 9 H8 9 H8 I 9 H 8 I M 9 H7M 9 H 7 9 H 7 M I 9 HAB,H M 23 
3 AB 59 HO 9 HTO 9 HT 29 HT 2 I 9HT4M9HT49HT5M9HT59HT5059HT5AM9HT5H9HT7M9HT7MI 9 M 24 
4HT6M9 HT8M 9 HT5A059 PH I 0» PH I TO 9 PH 1T2 I 9 PHITI 9 PH I T4M 9 PH I T4 9 PH I T5A 9 PHT5A M 25 
5 M 9 PHI 8 M 9 PHI T 8 M 9 PH 18 9 PH 12 I 9PH IT 6 M 9 PHI 7 M 9 PH I T 69 PH IT7M9PH0F6 9PH0H6*PT M 26 
60P09PT1PT09PT2T1S9PT2PT1 9 PT4PT2 9 PT5PT4 9 PT5T5A 9 PT8P8 9 PT 1 PT7 9 PT7P79 P M 27 
7T6PT5 9PT6T5D9PT5AM9PT55AD9P5 IH49 PSP I 49 PS IH59PSP15A 9PSIH5A 9 PSPI 8 9 PS M 28 
8 IH 89 PSIH 79 PSPI 69 PSIPHI 9 PSIH 6 9PS I H8I 9 PSPI 79 EN 9 ETA 29 ET 5 A 9 ETETS 9 ETA 2 S M 29 
9 9ET2D9EM7 9EM8 9W1 COR 9 W I DCOR 9 W4COR 9 W4DC0R 9 W7C0R 9 W7DCOR 9 W4T0T4 9 W4TEST M 30 
S9W1 9WIK9DH2M1 9 STTO 9 DELTO 9 DELTOD 9 GAMA 89 GAMA 79 VO 9 V 79 V89GAM7P1 9GAM7M1 M 31 
$9 0 NSRT 9 SFC 9 SFCC 9 0F4 9 OF 69 FOF 6 9 F 4 9 F5A 9 FAB 9 F 69 FNOVJI 9 FN 9 F 6 A 9 FOF4 9 F05A 9 M 32 
SW89RT59RTA9RT69RT5A9CP7»CP89CAPA89CAPA790WB9WCB9 ICOMB 9 I PP 9 iS9 ICOUN M 33 
ST 9 IKONT 9 JCONT 9 IC 9 KTEST 9 THETI 9 THETO 9 DWI 9 DW 29 DDHP 9 ZI 9Z29Z3«Z9ETAC 9 ET M 34 
$AO 9 FNOW 1 D 9 HT6 9 PH I H6 9 T5A 9 T4 9 RAT 9 WCOW 1 C 9 PBT AB 9 HOT 9 DTTAB 9 I PS 9 DTEMP 9 ET M 35 
$AP9ALTSAV9EMSAV9N0SAV9FNSAV9W2SAV9SFCSAV9N0 9 INOZZ 9 TT 2 TID 9 HTSA 9 PHTI M 36 
$9TT4TT2 9PHIT59PSP I59PHIT5M9HT1 9 PT2T 1 L 9PS 9 T A8 A7 9 T8VAL 9 T7VAL 9 NL I M 9 NU M 37 
SMBER 9 NSTART 9 NF I NAL 9 AVAL 9 XMNSLU9 T ABN 9 ETA4 9 ETA 5 9 TABPH5 9 TABH4 9 HPEXT 9 I M 38 


SCASE9 IPC9PT4T2D9HS9T5VAL M 39 

COMMON /NEW/ FD 9 FG 9 WF 9 OLDS 9 DFGS 9 FGl 9 FN 1 9 SFC 1 M 40 

COMMON /WWW/ WIDE M 41 

VO=EMACH*AO M 42 

A80W1= (0WB*0F6*RT6*T8 )/<P80P0*PO*V8 ) M 43 

FN0W1=0WB*0F6*(V8*CV/G)-V0/G+(P80P0-1 • )*A80W1*PO M 44’ 

FN=FNOWl#Wl M 45 

SFC=F6*0WB/FN0W1 *3600# M 46 

FD=Wl*VO/G M 47 

FG=FN+FD M 48 

WF=SFC*FN M 49 

IF UDDBS.EQoO) GO TO 1 M 50 

IF (V8.LT*1000) GO TO 2 M 51' 

CALL FTLUP { V8 9 DLDB 9 89 1 79 TABVJ 9 TABDLDB ) M 52 

DFGS=1 o-«005*DLDB M 53 

FG1=FG»DFGS M 54 

FN1=FG1-FD M 55 

SFC1=WF/FN1 M 56 

1 IF { IDDBS.GT 9 O ) GO TO 3 ■ M 57 

2 DLDB=0<» M 58 

FG1=FG M 59 

FN1=FN M 60 

SFC1=SFC M 61 



oon ooo O'Ui 4>cj 


IF (VO) 595»4 

VJE=( <FNOW1+VO/G>-»G)/(OW3-S-OF6) 

ETAP= (2o*G*FNOWl*\/0 ) / (0WB-»-0F6* ( VJE#*2-V0-»-i«-2 > ) 

ETAC= (Oli/B«-OF6-*^ ( V JE*-«-2- V0-^*2 ) ) / < 2 • *G*A J»OWB'»'F6^'‘- (H-HS) ) 

FNOWlD=FNt/WIDE 

ETAO=ETAP»ETAC 

GO TO 6 

ETAC=0 

ETApcO 

etao=o 

FNOWJD=FN/W10E 

RETURN 

END 


M 62 
M 63 
M 64 
M 65 
M 66 
M 67 
M 68 
M 69 
M 70 
M 71 
M 72 
M 73 
M 74- 


SUBROUT INE AFTER 
DIMENSION ALTSAV(15)« 
ISFCSAV ( 15 ) 
dimension TABTEM(78). 
1TABPT1U9)* TABNST(9) 


EMSAV(15)« NOSAV( 15 ) » FNSAVI 15) ♦ W2SAV(15)* 


TABH<78)« TABPHI(78) 
TABPT2(09)» TAB2(11 


TABPSH(78)< TABM(19)« 
TABW1<99)» TA8ETA(99) 


2 * ANS(5)« 
D I MENS I ON 
D I MENS I ON 
DIMENSION 
DIMENSION 


SFCC(2S)» TABPSP(78) 
PBTAB(8)» HDT(20)« DTTAB(20) 
TABPH5 ( 78 ) » TABH4 ( 78 ) 
TABN(12)» TABA7(12)* RAT<2)» 
TABVJ(17)» TABDLDB(17) 


T4 


3 ) 


INPUT 

COMMON /BKl/ ALT*EMACH«PT2TI D»TT4N»WC0W1 » WB ♦ A7 A7DN f PT65AD . A4A4DN » P 
iBOPStWl D«H?ETA2D.ENN«X»Y«TF»PT8PT6»CV*TT6MAX.PT6T5A.ETA6» IDES«VAL* 
2GAM0»R0*AJ.TSLS »PSLS,G«TABTEM«TABH»TABPHI .TABPSH.TABPSP.TABM* TABNS 
3T«TABPT2*TABZiTABWl 9 TABETA ♦ TABPT 1 9 1 ODBS 9 TABV J 9 TABDLDB 


COMPUTED VALUES 

COMMON /BK2/ WBOWl 9 AOAS 9 ABCOR 9 A80A7 9 A7A70 9 ANS» APO 9 A 4 A 4 D 9 A7» A8 9 AO 9 A 
180W1 tABEM7»T09T89T79T7I 9T3I 9 TTO 9 TT 29 TTI 9 TT 2 I 9 TT 2 TTI 9 TTS 9 TT5A 9 TT7TT 
2 1 9 TT 69 TT7TT6 9 TT 4 9 PO 9 P 8 OPO 9 PT7 9 HF 9 H8M 9H8 9 HS I 9 H8 I M 9 H7M 9 H7 9 H7M I 9 HAS 9 H 
3 AS 59 H 0 9 HT09HT2 9 HT 2 I 9 HT4M 9 HT4 9 HT5M 9 HT5 9 HT505 9 HT5AM 9 HT5H 9 HT 7M 9 HT7M I 9 
4HT6M9HT8M->HT5A059PH i 0 9 PH I TO 9 PH 1 T2 I 9 PH I T 1 9 PM I T4 M 9 PH I T4 9 PH I T5A 9 PH T5A 
5M 9 PH ! 8M 9 PH I T8M 9 PH ! 8 o PH I 2 I 9 PH I T6M 9 PH ! 7M 9 PH I T6 9 PH I T 7M 9 PH0F6 9 PH0H6 o PT 


N 

1 

N 

2 

N 

3 

N 

4 

N 

5 

N 

6 

N 

7 

N 

8 

N 

9 

N 

10 

N 

1 1 

N 

12 

N 

13 

N 

14 

N 

15 

N 

16 

N 

,17 

N' 

18 

N 

19 

N 

20 

N 

21 

N 

22 

N 

23 

N 

24 

N 

25 

N 

26 


Q S 
Hi a 

Ip 




A-34 


I 


2 


3 


4 


60P0«PT1PT0»PT2T1S«PT2PT1 « PT4PT2 « PT5PT4 » P r5T5A g PT8P3« PTl PT7« PT7P7» P 
7T6PT5 « PT6T5D »PT5AM « PT55AD, PS I H4 » PSP 1 4 » PS I H5 » PSPl 5 A « PS I H5A i PSP 1 8 « PS 
aiH8»PSIH7*PSP16«PSIPHl ?PS I H6 «PS I H8 1 « PSP 1 7 » EN » ETA2 < ET5A »ETETS g ET A2S 
9 » ET2D * EM7 « EM8 ♦ W 1 COR ♦ W1 DCOR « W4C0R » W4DC0R ♦ W7C0R g W7DCOR g W4T0T4 g W4TEST 
$g W1 ♦ WiK»DH2Ml » STTO * DELTO g OELTOD g GAMAS g GAMA7 g VO g V7 g V8 g GAM7P 1 gGAM7Ml 
$gONSRTgSFCgSFCCgOF4gOF6gFOF6gF4gF5AgFA6gF6gFNOWl g FN g F6A gF0F4 g F05A g 
£WSgRT5gRTAgRT6gRT5AgCP7gCP8gCAPA8gCAPA7g0WBgWCBg ICOMBg 1 PP g ISg iCOUN 
STg IKONJg JCONTg I Cg KTEST g THE! 1 gTHETOgOWl g DV\/2 g DDHP g Z 1 g 22 g Z3 g Z g ETAC g ET 
$AOgFNOWlDgHT6gPHlH6g T5Ag T4 gRATg WCOWl CgPBTABgHDT gOTTASg IPS gDTEMP gET 
SAPgALTSAVgEMSAVgNOSAVgFNSAVg WaSAVgSFCSAVgNOg I NOZZ g TT2T 1 D g HT5A gpHTl 
S,TT4TT2gPHIT5gPSPI5gPHIT5MgHTl gPT2TlLgPSgTA6A7gT8VALgT7VALgNLlMgNU 
SMBERgNSTART g NF I NAL g AVAL g XMNSLUg TABN g ETA4 g ETAS g TABPHS g TA8H4 g HPEXT g I 
SCASEg IPCgPT4T2DgHSgT5VAL 
CPS= IPS 
DO 4 I =1 g IPS 
PS= I 

IF (l-IPS) ig3g3 

TT 6 = (TT5A+50. ) + <TT6MAX-TT5A-S0. )*(CPS- 1 . -PS ) / < CPS-2 . ) 
PT6T5A=PT65AD-2.*10g**(-5)* ( TT6MAX-TT5A-50 , )*(CPS-1 . -PS )/( CPS-2 • ) 
CALL DlSCOT ( TT6 g TT6 g T ABTEM g.TABH g TABH g - 1 1 g 78 g 0 g HT6 ) 

CALL DI SCOT ( TT6 g TT6 g T ABTEM g TABPSH g T ABPSH g - 1 1 g 78 g 0 g PS I H6 ) 

FAB= (HT6-HT5A )/ <ETA6* (H-HS ) +HS-HT6-PS I H6+HF ) * ( 1 #+F5A*( 1 ,+ CPS IH5-PS 
1 IH5A)/(HT6-HT5A)) ) 

F6=F5A+FAB 
OF6=l 0+F6 

HT6M=HT6+PS I H6*F6/ ( 0F6 ) 

CALL DISCOT ( TT6 g TT6 g T ABTEM g TABPSP g TABPSP g - 1 1 g 78 g 0 g PSP I 6 ) 
call DISCOT (TT6gTT6gTABTEMgTASPHI gTABPHI g-1 1 g78g0gPHIT6) 

PH I T6M=PH I T6+PSP I 6*F6/0F6 
CALL NOZZ 
CALL ENG 
POWER=TT6 

CALL OOTPUN (POWER) 

GO TO 4 

TT6=TT5A 

PH IT6M=PHT5AM 

HT6M=HT5AM 

FAB=0 

PT6T5A=PT65AD 
F6=F5A 
0F6=1 g+F6 
GO TO 2 
CONTINUE 
return 


N 

27 

N 

28 

N 

29 

N 

30 

N 

31 

N 

32 

N 

33 

N 

34 

N 

35 

N 

36 

N 

37 

N 

38 

N 

39 

N 

40 

N 

41 

N 

42 

N 

43 

N 

44 

N 

45 

N 

46 

N 

47 

N 

48 

N 

49 

N 

50 

N 

51 

N 

52 

N 

53 

N 

54 

N 

55‘ 

N 

56 

N 

57 

N 

58 

N 

59 

N 

60 

N 

61 

N 

62 

N’ 

63 

N 

64 

N 

65 

N 

66 

N 

67 

N 

'68 

•N 

69 

N 

70 
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END 


71 


N 


> 

I 

CO 

ui 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE AT62SP (H«ANSeDTEMP»PBTAB »HDT« DTTAB ) 
answers compatible with those of AT62 
AT62SP CONVERTED TO CDC 1~67 

SUBROUTINE AT62SP WILL COMPUTE THE STANDARD DAY«1962 ATMOSPHERIC 
PROPERTIES (PRESSURE (.DENSITY »TEMP»SONIC SPEEDoAND REYNOLDS NO.PER 
MACH“FTo)AS A FUNCTION OF AN INPUT ALTITUDE (H«GEOMETRIC FT« ) 

FROM S.L. TO 200000 FT«o 

IN ADDITION A NON-STANDARD DAY MAY BE COMPUTED IN SEVERAL OPTIONAL 
WAYS# 

BY INPUTING A TEMPERATURE I NCREMENT ( DTEMP« DEG«KELV I N ) A HOT DAY 
UTILIZING THE STANDARD DAY TEMPERATURE GRADIENTS MAY BE COMPUTED# 
THIS ATMOS .WILL USE THE HYDROSTATIC EQUA.AND THE EQUA.OF STATE# 

BY INPUTING A TABLE OF TEMPERATURE INCREMENTS AS A 

FUNCTION OF ALTITUDE CDTTAB VS HDT) A HOT DAY MAY BE COMPUTED 

WHICH WILL OBEY THE EOUA.OF STATE BUT NOT THE HYDROSTATIC EQUA# # 

COMBINATIONS OF THESE TWO OPTIONS MAY BE USED TO COMPUTE 
A WIDE VARIETY OF NON-STANDARD ATMOSPHERES# 

FOR EITHER ATMOS# DESIRED# A TABLE OF PRESSURES AT SPECIFIED 
ALTITUDES MUST BE INPUT (PBTABCB))# THIS TABLE WILL BE SUPPLIED 
BY ACG FOR SELECTED DTEMP VALUES. 

(THE BASIC TEMPS.AND GRADIENTS ARE BUILT IN > 

DIMENSION ANS(5># HMTAB(9)» HBTAB(0)# ELMTAB<0)# TMBTAB{8) 
DIMENSION PBTAB(8)# HDT(20)# DTTABC20) 

DIMENSION ZDPRCl)# HDPR < 1 ) 

DOUBLE PRECISION ZDPR # HDPR « Z6 « Z6P « TEMP2 

DOUBLE PRECISION TEM2 « TEM3 « TEM4 , # TEM5 # TEM9 

DOUBLE PRECISION S 1 NZ # DSIN # DLOG # DCOS 

DATA (HMTAB ( I ) # 1 = 1 #9 ) # CHBTA3 ( I ) # 1 = 1 # 8 ) « (ELMTA 8 ( I ) # 1 = 1 # 8 ) # (TMBTAB ( I 
1 ) # 1 = 1 #8 >/ -99999999# # 1 1000# #20000# # 32000 ## 47 OOO , #52000# #61000# #7900 
20# #99999999# #0# # 1 1000# # 20000 ## 32000 # # 47000 ® # 52000 ## 61 000# #79000# «- 
3o0065# »0» oOOl q s0o2S# oO#—« 002 # — o004# 0 Q v23B » 15#216865#216# 65 #228#65o 
4270 #65 >270 #65,252 #65 9 ISO #65/ 


O 1 
O 2 
' 0 3 

O 4 
O 5 
0 6 
O 7 
O 8 
0 9 

O 1 0 
O I 1 
O 12 
O 1 3 
O 1 4 
O 15 
0 16 
0 17 

O 18 
O 19 
O 20 
0 21 
0 22 
O 23 
0 24 

0 25 

O 26 
O 27 
O 28 
O 29 
0 30 

O 31 
O 32 
0 33 

O 34 
O 35 
O 36 
O 37 
0 38 


II 

a 

^65 



(jfo— n o nnonooo oooooo 


C 


EMOR=®00343368 
GAM=1 «4 
G0=9o80665 
CONVl = 0 3048 
C0NV2= » 45359237 
ZFT = H 

convert z ft. to z meters 

Z=ZFT*C0NV1 

COMPUTE H DOUBLE PRECISION 

1111 CALL STORE<ZDPR.ZDPR{ 1 ) .1 .2 ) 

ZDPP=Z 

CALL STORECZDPR .ZDPR (2 ) . 1 .0.0 ) 

Z6 = ZDPR < I )*. 000001 DO 
Z6=ZDPR*.00000 1 

210 Z6P=Z6+6.3675708D0 
Z6P=Z6+6. 3675708 

21 1 TEMP2=, 001 05591 79D0*Z6+1 .5640943DO 
TEMP2=. 00 10559179*26+1 .5640943 
TEM2=Z6P*Z6P 

TEM3=TEM2#Z6P 

T£M4=TEM2*TEM2 

TEM5=TEM2*TEM3 

TEM9=TEM4*TEM5 

S I NZ=DS I N < . 00052795893D0*Z6+ .782047 1 3D0 ) 

SINZ=S INI .00052795893*26+. 7820471 3) 

217 HDPR ( 1 >=-40648635. D0/(Z6P 1-894899. 28DO/TEM3-16302410.DO/TEM9-7 
1 3724D0/TEM4-1 35. 55760D0*Z6*Z6-1 724.849000*26+ (2684697. 8D0/TEM3 )* ( 
2 INZ*S INZ > + (70S .70509D0/TEM2-256757.85D0 

3Z6- 1 634923 . 8D0 + . 000790 1 7959DO*DLOG ( ZDPR ( 1 > + 6367570 .6D0 > ) * ( DS I N ( TE 
4P2> )+( .7483343600/ { Z6P ) -243 1 60820 .DO )*( DCOS ( TEMP2 > >+9646530* 

5284232D0 

HDPR=-40648635./(Z6P >-894899. 28/TEM3-1 63024 1 0./TEM9-7408.3724/TEM4 
1 -135. 55760*Z6*Z6- 1724. 8490*26+ (2684697.a/TEM3)*<S INZ*SINZ )+ (708.70 
2509/TEM2-256757 .85*26-1 634923.8+. 000790 1 7959*AL0G ( ZDPR ( 1 >+6367570. 
38 ) >* (SIN (TEMP2 ) )+ (. 74833436/ ( Z6P > -243 1 60820 . >* (C0S(TEMP2 > >+9646530 
4.7284232 

H =PART (HDPR.HDPR ( I > . 1) 

H=HDPR 

IF (Z“90000.0) 1.1,13 

Z LESS THAN 90000. M FIND H IN HMTAB TABLE 
N=1 

IF (H-HMTAB(N)> 5.4.3 
N = N+1 
GO TO 2 
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O 42 
O • 43 
O 44 
O 45 
O 46 
O 47 
O 48 
O 49 
O 50 
O 51 
O 52 
O 53 
O 54 
O 55 
O 56 
0 57 

O . 58 
O 59 
0 60 
0 61 
O 62 
O 63 
0 64 

O 65 
O 66 
O 67 
0 68 
O 69 
O 70 
0 71 

O 72 
O 73 
O 74 
O 75 
O 76 
O 77 
0 78 

0 79 

O 80 
0 81 
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vO CD O O' O U1 


> 

■ 

to 


4 


1 0 
C 
1 1 

C 

12 

13 
C 

14 
C 

15 

16 
17 
C 

C 

c 

c 

c 


J = N 

O 

83 

GO TO 6 

0 

84 

J = N-1 

0 

85 

FIND H8«P8»ELM»TMBfl TM 

0 

86 

HB=H3TAB ( J ) 

0 

87 

PB=PBTAB( J) 

0 

88 

ELM=ELMTAB ( J ) 

0 

89 

TM8=TMBTAB ( J )+DTEMP 

• ' 0 

90 

TM=TMB+ELM* ( H-H8 ) 

0 

91 

IF <HDT ( 1 ) oLToO* ) GO TO 7 

o 

92 

call discot (2FT»ZFT,HDT,DTTAB»DTTAB*-01 1 ♦20»00*OT) 

0 

93 

TM IN KELVIN 

0 

94 ' 

TM=TM+DT 

• 0 

95 

TMB=TMB+DT 

o 

96 . 

IF (22640o-PBTAB( 2 ) ) 10»10tS 

0 

97 

IF (PBTAB(2 )-22630* ) 10tl0»9 

• 0 

98 

DTSAVE = TMB-TMBTAB ( J ) 

0 

99 • 

TN=TM-DTSAVE 

0. 

100 

TMB=TM3-DTSAVE 

0 

101 

IF CELM) 11«12»I1 

0 

102 

IF ELM NOT =0 COMPUTE P 

0 

103 

P = E><:P< ALOG (PB )- (GO'«'EMOR/ELM) *AL0G(TM/TM8 ) ) 

0 

104 ■ 

GO TO 14 

0 

105 

IF ELM=0 COMPUTE P 

o 

1 06 

P = EXP ( ALOG ( PB ) - ( GO*EMOR ) * < < H-HB ) /TMB ) ) 

0 

1 07 

2 = Z 

0 

108 

convert P NT/M2 TO PSF 

- O 

109 ■ 

PRESS=P/47 o880 1 83 

o 

1 10 

T degrees KELVIN 

0 

111 

IF (22640. -PBTA8<2 ) ) 17*17.15 

0 

112’ 

IF. (PBTAB(2>-22630, ) 17*17*16 

. o 

1 13 

tm=tm+dtsave 

0 

114 

T = TM 

0 

US 

RHO KG/M3 

. 0 

1 16 

RH0=EM0R*(P/TM ) 

0 

M7 

CS M/SEC 

0 

1 1 8 

CS=SQRT <GAM*TM/EM0R ) 

0 

1 19 

CONVERT RHO KG/M3 TO LB/FT3 

0 

120 

RH0=RH0*C0NV 1 ^CONVl *CONV 1 /CONV2 • 

0 

121 

CS FT/SEC 

0 

122 

CS=CS/CONVi 

0 

123 

VISCOS IN LB/FT-SEC 

0 

1 24 

C0NST= ( CONV 1 /C0NV2 ) ^ 1 * 458/ 1 0 o o 

0 

125 

VlSCOS=CONST«T^ti«'l .5/ (T + l 10o4 ) 

o 

1 26 
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C RNOML IN RN/MACH-FT 

RNOML=RHO*CS/V I SCOS 
C convert RHO LB/FT3 TO SLUG/FT3 

RHO=RHO/,32174O405E+O2 
AnS( 1 )=RHO 
ANS<2 ) =PRESS 
ANS(3)=T 
ANS(4 )=CS 
ANS(5)=RN0ML 
H=2FT 
RETURN 
END 
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O 133 
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O 135 
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O 137 
O 138- 
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2 

3 

4 


7 
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SUBROUTINE FTLUP (X » Y » M*N» VAR I . VARD ) 

P 

1 

DIMENSION VARMl)* VARD(l)* V(3)« YY(2) 

P 

2 

DIMENSION 1 I (43) 

P 

3 

DATA ( I I { J) * J=1 »43)/43*-l/ 

P 

4 

MA= TABS (M ) 

P 

5 

LI=MOD(LOCF(VARI (1 ) ) ,43)+l 

P ■ 

• 6 

1 = I I (LI ) 

p. 

7 

IF (I »GE*0 ) GO TO 6 

' P 

8 

IF (N.LT*2) GO TO 6 

• P 

9 

IF (VARl (21-VARI ( I ) ) 2«2.4 

P 

1 0 

K = L0CF< VARI ( 1 ) ) 

P 

1 1 

PRINT 17* J*K* (VaRI ( J) . J=1 *N) » <VARD( J) *U=1 *N) 

P 

12 

STOP 

P 

1 3 

DO 3 J=2«N 

P 

14 

IF (VARI ( J)-VARI { J-1 ) ) 3*1*1 

P 

15 

CONTINUE 

P 

16 

GO TO 6 

P 

17 

DO 5 J=2*N 

P 

18 

IF (VARI ( J)-VARI (U-1 ) > 1*1*5 

P 

19 

CONTINUE 

P 

20 

IF (I.LEoO) 1=1 

P 

21 

IF (I.GEoN) I=N-1 

P 

22 

IF (N.LE* 1 ) GO TO 7 

P 

23 

IF (MAoNEoO) GO TO 8 

P 

24 

Y=VARD( 1 ) 

P 

25 

GO TO 16 

P 

26 

IF ( (VARI ( n-Xl^^IVARl ( I + l )-X) ) 11*11*9 

P 

27 



A-39 


9 

IN = S IGN < 1 .0» (VAR I (1 + 1 )-VARI ( I ) )* (X-VARI ( 1 ) ) ) 

, P 

28 

10 

IF ( ( i+in)»le.o> go to H 

’ P 

29 


IF ( ( I+IN).GE*N) GO TO U 

P 

30 


I=T+IN ' • 

P 

31 


IF ( (VARI ( I )-X )*( VARK I + l >-X > ) 11»11»10 

P 

32 

1 I 

IF (MA.EQo2 ) GO TO 12 

P 

33 


Y= (VARD( I )*(VARI (I + l )-X)-VARD( I + i )^MVARI ( I )-X) )/(VARl (I + l )-VARI ( I ) 

P 

34 


1 ) 

P 

35 


GO TO 16 

P- 

•36 

12 

IF (N.EQ.2) GO TO 1 

P 

37 


IF (IoEQ»(N-l)) GO TO 14 

P 

38 


IF ( I oEQ. 1 ) GO TO 13 

P 

39 


SK = VARI ( 1 + 1 )-VARl { I ) 

P 

40 


IF ( (SK*(X-VAR! ( I-l > ) )oLTo (SK*(VARI ( I+2)-X) ) ) GOTO 14 

P 

41 

13 

L=I 

P 

42 


GO TO 15 ■ 

P 

43 

14 

L= I-l 

P 

44 

15 

V ( 1 ) =VARI (L )“X 

P 

45 . 


V (2 >=VAR1 (L+1 )-X 

P 

46 


V(3)=VAR! (L+2 )-X 

P 

47 


YY(1 )=(VARD(L)*V(2)-VARD(L+1 )*V( 1 ) )/(VARI (L+1 )-VARl (L> ) 

P 

48 


YY (2 )= ( VARD (L+i )* V ( 3 )- VARD ( L+2 > ( 2 ) )/ ( VARI (L+2 ) -VAR I (L+1 ) ) 

P 

49 


Y=CYY( 1 )*V(3)-YY(2)*V(1 ) )/(VARI (L+2) -VAR I (L) ) 

P 

50 

16 

I I (LI ) = I 

R 

51 


RETURN 

P 

52 

C 


P ■ 

53.- 

17 

FORMAT { 1 HI < *TABLEBEL0W0UT0F0RDERF0RFTLUPATP0S 1 T I ON* » 1 5 » /+XTABLE I S 

P 

54 


lSTORE0INL0CATI0N*,O6«//(8G15.8) ) 

P 

55 


END 

P 

56- 


SUBROUTINE DlSCOT < XA« ZA, TABX ♦TABY » TABZ« NC *NY ♦ NZ» ANS ) 

C THE DIMENSIONS IN THIS SUBROUTINE ARE ONLY DUMMY DIMENSIONSo 

DIMENSION TABX(2)» TABY(2)« TABZ(2)« NPX(8)» NPY(8)« YY ( 8 ) 

C DIMENSION TABX (2 ) »TABY(2 ) ♦TABZia ) ^NPX(8) tNPYIS ) »YY (8) 

CALL UNS (NC» IA» IDX» IDZ« IMS) 

I F (NZ-1 ) 1 1 1 <2 

I CALL DISSER ( XA •. TABX { 1 > ♦ 1 t NY » I DX . NN ) 

NNN= IDX+1 

CALL LAGRAN ( XA -» TABX t NN ) ^ T ABY (NN ) oNNN « ANS ) 

GO TO 12 


Q 1 
Q 2 
Q 3 
0 4 
O 5 
Q 6 
Q 7 
Q 8 
Q 9 
Q 10 



n 


2 


ZARG=ZA 
1P1X=I0X+1 
IP12=1D2+1 
IF (lA) 3»5*3 

3 IF (ZARG~TABZ(N2) ) 5.5»4 

4 ZARG=TABZ (NZ ) 

5 call DISSER <2ARG*TABZ( 1 ) . 1 «NZ» IDZtNPZ ) 

NX=NY/NZ 

NPZL=NPZ+IDZ 
1 = 1 

IF (IMS) 6»6«8 

6 CALL DISSER (XA »TABX ( 1 ) « 1 »NX * I DX » MPX ( 1 > ) 

DO 7 JJ=NPZ,NPZL 

NPY ( I ) = < J J-1 ) *NX+NPX ( 1 ) 

NPX( I )=NPX( I ) 

7 1=1+1 

GO TO 10 

8 DO 9 JJ=NPZ»NPZL 
IS= ( JJ-1 )*NX+1 

call DISSER (XAtTABXd ) ♦ !S^NX» IDX«NPX(I ) ) 

NPY( I ) = NPX( I ) 

^ 9 • 1=1+1 

10 DO 11 LL=I»IP1Z 

o NL0C=NPX(LL) 

NL0CY=NPY(LL ) 

11 CALL LAGRAN ( XA ♦ TABX (NLOC ) « TABY ( NLOC Y ) d P 1 X » YY < LL ) ) 
CALL LAGRAN ( ZARG dABZ ( NPZ ) ♦ YY < 1 ) i I P 1 Z i ANS ) 

12 RETURN 
END 


Q 

1 1 

, Q 

12 

0 

13 

Q 

14 

0 

15 

Q 

16 

.0 

,-r, 1 7 

o 

18 

Q‘ 

19 

Q 

20 

. 0 

21 

Q 

22 

Q 

23 

Q 

24 

Q 

25 

Q 

26 

■0 

27 

Q 

■ 28 

Q 

29 

Q 

30 

Q 

31 

Q 

32 

Q 

33 

Q 

34 

0 

35 

Q 

36 

Q 

37 

’q 

38 

■Q 

39' 


SUBROUTINE UNS ( I C ♦ I A ♦ I DX d DZ » I MS ) R 1 

IF ( IC) I d »2 R 2 

1 1MS=1 R 3 

NC=-IC R A 

GO TO 3 R 5 

2 IMS=0 R 6 

NC=IC R 7 

IF (NC-lOO) 4»5»5 R 8 

IA=0 R 9 

GO TO 6 RIO 



O' 


5 IA=1 
NC=NC-1 00 

6 IDX=NC/10 
ID2=NC-IDX*l0 
RETURN 

END 


R 1 1 
R 12 
R 13 
R 14 
R 15 
R 16- 


SUBROUTINE DISSER (XA«TAB» I »NX» ID,NPX) 
DIMENSION TAB<2) 

C DIMENSION TAB (2) 

NPT= ID+1 
NPB=NPT/2 
NPU=NPT-NP8 
IF (NX-NPT) 2*1.2 

1 NPX=I 
RETURN 

2 NLOW=I+NPB 
NUPP=I+NX-(NPU+1 ) 

DO 3 II=NLOW.NUPP 
NLOC=I I 

IF (TAB(II)-XA) 3.4.4 

3 CONTINUE 
NPX=NUPP-NPB+1 
RETURN 

4 NL=NLOC-NPB 
NU=NL+ID 

' ■ DO. 5 JJ = NL.NU 

NDIS=3J 

IF (TAB ( JJ )-TAB { JJ+1 > ) 5.6.5 

5 CONTINUE 
NPX=NL 
RETURN 

IF (TAB(NDIS )-XA ) 8.7.7 
NPX=NDIS-ID 
RETURN 

8 NPX=NDIS+1 

RETURN 
END 


S 1 
S 2 
S 3 
S 4 • 
S 5 
S 6 
S 7 
S 8' 
S 9 
S 10 
S 1 1 
S 12 
S 13 
S 14 ' 
S 15 . 
S 16 
S 17 
S 18 
F 19' 

S 20 
S 21 
S 22 
S 23 ■ 
S 24 
S 25 
S 26 
S 27 
S 28 
S 29 
S 30 
S 31- 
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SUBROUTINE LAGRAN ( XA t X ♦ Y » N ♦ ANS ) 
DIMENSION X(2) ♦ Y(2) 

C DIMENSION X(2) ,Y(2) 

SUM=0.0 
DO 3 I =1 »N 
PROD=Y( I ) 

DO 2 J= 1 *N 
A=X( I )-X( J) 

IF (A) 1*2«1 

1 B= (XA-X( J ) >/A 

prod=prod*b 

CONTINUE 
SUM=SUM+PROD 
ANS=SUM 
RETURN . 

END 


T 1 
T 2 
T 3 
T 4 
T 5 
T 6 
T 7 
T 8 
T 9 
T 10 
T 1 1 
T .12 
T 13 
T 14 
T 15 
T 16- 


> 

I 

ro 


SUBROUTINE MTLUP (X » Y » M » N» MAX »NTAB » I . VAR I » VARD ) 

DIMENSION VARI<1)» VARD(MAX»1)» Y(l)» V(3)» YY(2) 

LOGICAL EX 

IF (M.EQ*0) GO TO 17 

IF (N.LE* 1 ) GO TO 17 

EX=.F. 

IF ( I.GE.O) 60 TO 6 

IF {N.LT*2) GO TO 6 

IF (VARI {2)-VARIU ) ) 2*2»4 

1 K=LOCF( VARI < 1 ) ) 

PRINT 19. J.K. (VARI ( J) . J=1 .N) 

STOP 

2 DO 3 J=2.N 

IF (VARI ( J)-VARI ( J-1 ) ) 3.1.1 

3 CONTINUE 
GO TO 6 

4 DO 5 J=2.N 

IF (VARI (J)-VARI (J-1 ) ) 1.1.5 

5 CONTINUE 

ft IF (loLEaO) 1=1 


U 1 
U 2 
U 3 
U 4 
U 5 
U 6 
U 7 
U 8 
U 9 
U 1 d 
U 1 1 
U 1 2 
U 1 3 
U 14 
U 15 
U 1 6 
U 17 
U 18 
U 19 
U 20 



• 3 > 

I 

■?=• 

W 


7 

8 


9 
1 0 

1 1 


12 


13 

14 

15 


16 


17 

18 

C 

19 


IF (I«GE«N) I=N-1 

IF ( (VARI n >-X>* ( VARI ( I +1 )-X ) ) 10.10f7 
IN=SIGN< 1 .0» (VAR I (I + I )-VAR 1(1) )* (X-VARI ( 1 ) ) ) 

IF ( ( I+IN)*LE#0 > GO TO 9 
IF { ( 1+ IN ) .GE.N ) GO TO 9 
I=I + IN 

IF < (VARI (I )-X)*(VARI C I+l )-X> ) 10tl0«8 
EX=.T. 

IF (M.EQ.2) GO TO 12 . 

DO 11 NT=1»NTAB 

Y(NT )= ( VARD { I »NT )* { VARI <1 + 1 )~X )-VARD ( I + l .NT )* (VAR I ( I )~X) )/<VARI ( 1 + 
1 1 l-VARI < I ) ) 

IF (EX ) 1 = 1+ IN 

RETURN 

IF (N.EQ.2) GO TO 1 
IF (I*EQ«(N-1)) GO TO 14 
IF (I .EQ. 1 ) GO TO 13 
SK-VARl ( I+l )-VARl ( I > 

IF ( (SK*(X-VARI ( I -I ) ) >.l-T« (SK*(VARI ( 1+2 )-X ) ) ) GO TO 14 
L= I 

GO TO 15 
L=I^1 

V ( 1 )=VARI (L)-X 
V (2 )-VARI (L+1 )-X 
V<3)=VARI (L+2 )-X 
DO 16 NT=1.NTAB 

YY ( 1 ) = ( VARDIL.NT )*V <2 ) -VARD (L+1 .NT l-^Vl 1 ) >/( VARl (L+l )-VARI (L) ) 

YY(2) = (VARD(L+1 .NT ) * V ( 3 ) -VARD CL+2 . NT ) *V ( 2 ) )/(VARI (L+2)-VARl (L + 1 ) ) 

Y (NT )= ( YY< 1 )*V(3)-YY (2)*V ( 1 ) )/ ( VARI (L+2 )-VARI (L) ) 

IF (EX) I=I+IN 
RETURN 

DO 1 8 NT=1 .NTAB 
Y(NT )=VARD( 1 .NT ) 

RETURN 

FORMAT ( IHl .*TABLEBELOWOUTOFORDERFORMTLUPATPOSITION#. 15. /*XTABLE I S 
1 STORED INLOCATION*. 06.// (86 15 .8) ) 

END 


u 

21 

u 

22 

u 

23 

'u 

24 

u, 

,••25 

M 

26 

u 

27 

u 

28 

u 

29 

u 

30 

u 

'31 ■ 
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APPENDIX B, 
SAMPLE PROBLEMS 


B-1 



B-2 


The examples, which follow, were chosen to show the proper sequence for 
establishing an engine design and to show samples of the three modes 
of printed output. 

Example 1 (Case 1) - The following input (ENPUT) data is required to 
establish an engine design. 

SENPUT 


ALT=0. »CV=.985»OTEMP=0, «EMACH=0. »ENN=1 . , ETA2D= .875 • ETA4 = .98 » ET A5= . 9« ETA6= • 9 » 
HPEXT = 0. «PT2T1D=15. , PT4T2D= *95 1 PT65A0= 1 . «PT8PT6=1 • .TT4N = 3060. » TT6MAX=4060. » 
W1D =700. ♦ WB = 0. . WCOWl =.0001 . WIDE = 700. »X=. 1 tXMNSLU=l .5. Y= 1 .1 * 

AJ=778. 156.G=32. 1 74 « GAMO= I . 4 , PSLS=2 1 16.2381 ,R0=53. 33032 1 . TSLS=518. 66999. 

AVAL=.005. ICASE=1 »NL I M=20 ,T5VAL= .001 .T7VAL=.0l «T8VAL=.000l . VAL=.025. 

PBTAB=101325..22632.025.5474,8705.868,01404. 1 10,90547. 
59.000452.18,209893.1.0376983. 


T ABM = 0 . ..2* .4. .6. .6.1 , .1 .2. 1 ,4. 1 .6. 1 ,8.2. .2.2. 2, 4.2, 6. 2, 8. 3, .3, 2. 3. 4. 3, 6. 
TABPTI=19*1 . . 

TABETA=,7609. ,8045. ,86. .901 . .87. .797. ,777. .721 . ,695 ., 7538. ,81 39 . ,876 . .927 ., 909. 
.83. .807. .753. .735. .7409. .8245. .89. .951 . . 942. .886. .871 . .807. ,783. .7232. .8292. 
,901 . ,968. ,974. ,933. ,926. ,87. ,815. .7091 . .8316. .912. ,984. 1,001. ,973., 965. , 903. 

, 839 . , 6867 . .8316. .922. .991 . 1 ,0l6. 1 .005. .989. «919* *846. ,6714. .8 292 ♦ .923. .995. 

1 ,025. 1 .025. 1 ,005. .931 . . 848 . ,6478 . ,828 . , 91 9. ,994. 1 ,028. 1,031.1 . 0l 3 . , 939 . , 847 . 
,6207. ,8174. ,91 . ,989.1 ,025. 1 ,029, 1 , 0 1 3. , 942 . ,837 , ,5971 . ,8045. , 895 . ,978. 1,015. 
1,021.1 .009. .937. .‘806. . 5689 .. 7845 . ,878. , 961 .1,014.1.007.1,. .919. .753 , 

TABMST=»4. .5. *6. ,7. ,8. .9, 1 . . 1 . 1 , 1 ,2. 


TABPT2=. 1648. .2235. .3409. ,5038. .6685. ,84 . 1 , . 1 , 1 05. 1 ,1695. 

TABWl ( 1 )=. 1959. .276 . ,4027. ,5975. ,7751 . ,9428,1 ,0359. 1,108.1,165. 

, 1 948. ,2743. .4005, ,5954. .7751 , .9428. 1 ,0359,1 ,108.1,165. .1938. .2724, 

.3994 , ,5943. .7751 , ,9428. 1 ,0359. 1.108.1 ,165., 1927. ,2704 . ,3979. . 591 . ,7748. ,9428. 
1,0316.1,108.1.165, ,1905. ,2685. ,394. .5878. .771 8. ,9396. 1 ,031 6. 1 . 1 08 . 1 , 1 65 . ,1884. 
,2663. ,3908. ,5845. .7675, .9353, 1 ,0305. 1 , 1 08 . 1 , 1 65 . , 1 862. ,2639. , 3875 ♦, 577 .. 762 1 ♦ 
,9299. I .0305.1 .1 08 . 1 , 1 65 , , 1 84 ♦ .2607 3843 .. 571 6. , 7532 , .9221 . 1 ,0305.1 ,106.1,165. 
, 1819. .2565. ,3778, .5629, ,7415. ,9060, 1 ,0284 .1,103,1 , 1 63 . . 1 797 . , 251 1 ,, 3735 ...552 1 . 
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TABZ=0 




TABH( 1 )= 0««12i.40» 1 92 .8 ♦ 204 • 7 . 2 1 6o 65 * 228 a62 » 240 0 62 « 252 o64 » 264 , 7 . 276 , 79 , 288. 88 ♦ 
301 «O5*313.27»325.52«337086*35O027»362«75»375.29»387.92»4 OOo63941 3.42«426.28, 

4 39.24 9 452.26 9 465.38 « 478. 56 9 491 .81 9 505. 1 4 9 51 8 . 56 » 532. 02 9545.56 9 S59 . 1 7 9 572 . 84 9 
566.57.600.35 9 614.20 9 628. 1 0 . 642 . 04 9 656. 04 9 670 . 09 . 684. 20 . 698. 35 9 71 2 .53 9 726 0 76 9 
741 .039755.359769.709784.089798.509812.969827.459841 .98 9856.53 9871 .129385.739 

900. 38 99 15. 05 9 929.75 9944. 48 9959 .23 9974. 00 9 986.80. 1 003.63 . 1 0 1 8. 48 9 1033.36. 

1 048.25.1 063.16.1 078.09. 1093.05. 1108.03.1123.02.1138.04.1153.08.1168.13. 

1183.20.1198.29. 1350.09. 1503.36. 

TA8PHI <1 )=0. . 1 .30757.1 .49849. 1 *53031 . 1,55845.1 .58368.1.60656. 1 .62748.1 .64679. 

1 .66471 .1.68146.1 .69717. 1 ,71200.1 .72598. 1 .73932*1 .75206. 1 .76423. 1.77591. 

1 .78713. I .79796. 1 .80840. 1 ,81846. I *82823. 1 .83771 .1 .84691 . I .85584. 1 .86454. 

1 .87300. 1,88128. 1 .88932. 1 ,89716. 1 ,90486.1 .91235.1 .91967. 1 .92684.1 .93384. 

1 .94071.1 .94743. 1 .95402. 1 .96048. 1 .96682. 1 .97304 . 1 .97914 . 1.98513.1 .99102. 

1 .99680. 2. 00249. 2.00807. 2. 01 357. 2.0 1897. 2. 02494. 2.02950* 2 .03467. 2. 03974. 

2 .04474 .2*04 967.2 . 05452 .2 ,05930 .2.06401 « 2 • 06866 . 2 • 07324 « 2 .07776 *2 . 08222 ♦ 
2.08662.2.09096.2.09524.2,09948.2, 10365.2. 10778-. 2.1 1185.2,1 1588.2.11986. 

2 .12379.2.1 2767.2.13151 .2,1 3531 .2.1 7106.2.20336. 

TABTEM (1 1=0. .50. .350 • . 400 . . 450. . 500. .550 • ♦ 600. .650. .700 • .750. . 300. .850 • .900, . 

950. . 1000. . 1 050. ♦ 1 1 00, . 1 150. . 1200. . 1250. . 1 300 . ♦ 1 350 . , 1 4OO . , 1 450 . * 1500. . 1550. * 

1 600. . 1650. . 1 700. . 1 750 . . 1800. . 1850. . 1900* .1950. .2000. . 2050 . . 2 1 00 . . 21 50* .2200, . 

2250. . 2300. . 2350. .2400 . .2450. .2500. .2550. .2600. .2650. .2700. . 2750 .. 2800. .2850. . 

2900. . 2950. . 3000. .3050. .3100. .3150. .3200. .3250. .3300. .3350. .3400. .3450* .3500. . 

3550. . 3600. . 3650 . .3700 . . 3750. . 3800. . 3850. .3900. .3950. .4000. .4500. . 5000 « . 

TABA7=1 0 01. !• 05. 1.1. 1,25. 1.1 75. 1.2. 1.25. 1.3. 1.35 .1.4. 1.45. I, 5. 

r 

. rABN=l . . , 975 . . 95 , o 9 , 0 . .85. ,6. .75. . 7 , ,65,. 6. . 55 . .5. 


TABDLDB ( 1 ) =6 . 1 5 » 6 . 25 .6 , 55 . 6 . 7 . 6 • 95 . 7 , 3 . 7 . 85 * 8 . 5 . 9, 3 » 1 1035 . 13 . 5 . 14 . 5 . 

1 5 . 35 . 15. 95 . 16.35. 16.95. l7o2. 

TABVJ ( I )=1000, , 1200. . 1400, . 1500o . 1600, ♦ 1700, ,1800. , 1900. .2000, .2200. . 

2400. .2500. .2600. .2700 . .2800, .3000. ,3200. . 

IDES=0» IDDBS=0. IPC=2.ZTHRUST-0. . IFUEL=1 . A4A4DN=1 . oA7A7DN=l . , IPS=5 . NF INAL= 1 2 . 
NSTART=1 .NUM3ER=1 ,P80P8=1 *. IPUNCH=0. 



£ 



i?-a 


It should be noted that to establish an engine design that it must be 
done at sea level static standard day conditions and all parameters 
associated with installation effects must be zeroed out, that is the 
following parameters must be as specified below; for a non-afterburning 
engine. 


ALT = 0. 
DTEMR = 0. 
EMACH = 0. 
HPEXT = 0. 
PT65AD = 1 . 
PT8PT6 = 1 . 


WB = 0. 

TABPTl = 19*1. 
IDES = 0. 
IDDS = 0. 

IPC = 2 
ZTHRUST = o! 


For the afterburning design case values which are typical of after- 
burner pressure losses must be specified. That is PT65AD and PT8PT6 
should have a value less than one and the control IPC = 1 . 


Example 1 is for a non-afterburning design point case with the long 
form printed output as follows: 


-■ • INPUT CASE 

1 

PC 2 

P3 



I 




AIT 0. 

. .. — 

HACH 0. 


HPEXT 

0 0 

■ ■ IDOBS 0 

- 


■ WCOWl .1000CC.00E 

“03 

WB 0. 


CV 

.9S50C000E+00 ZTHRUST 0. 



— DESIGN CASE 

PTZTID 

.150000COE+02' 

WIDCOR .10193376S+0] 

Z . 

89247312E+00 ' ONSRT 

.lOCOOOOOE+Ol 

ETAZ .67500000E 

+00 PTIPTO .lOOOOOOOE+Ol 

TTO 

.51867000E+03 

TT2 ,lie98770£+04 

TT4 

,306C0q00E' 

DH2-1. .16A95593E+03 

ETAS ;90000000E+00 

, 

_ 




1 . - 

' - HTl • .23310080E-(-03 

- r 

HT2I 

.377A372AE+03 




■* ’ 

* ‘ 

FN 73^66.09100 


.97996 

Fn/W1 

10^,95156 

PTl/PTO 

l.COOOO 

PT2/PT1 15,00000 ' 

ETA2 

,87500 

WIK I 

N/SRT 

■ i.OOCOO 

TTO 

51S. 67000 

TT2 

1189,87698 

TT4 ■ 3060.00000 

PT4/PT2 

.95000 

TTS 2'^GQ. 03123 

PT5/PT4 

',35665 

TT5A 

2339,70308 

TT6 

2384.70808 

PT6/PT5A- 1,00000 

F6 

.02857 

T8 1623,0<*063 

PT87P8 

5.08232 ■ 

VO 

O.CQOOO 

V8 

3332.91509 

M8 1.73554 

j ^ 

1,00000 

AA7A^iD ■ '1.00000 

A7/A7D ' 

■ l.COOOO 

AS/ A7 

1,41268 

A3 ■' 

1920,38657' 

■" "CAPAS 3.S4<^-45 

Wi 

700,00000 

ET.AP' ■ "0,00000 

ETAC ■ 

0,00000 ' 

ETAQ 

' ■ 0.0000'J 

FM/WIO 

104.95156 

HF 71994. 10768 ' 

■■ ■ WG" 

719,99836 

■ PT7P7 ■ ■ 1.82262 

PT1PT7 

.19676 PT5T5A 

1,00000 

■ ■ PT6T5A 

1.00000" 

PT8PT6 1.00000 

■ ‘ ■ TTl 

■ 513,67000 

WB O.OOOwO 

WBDWl 

■ ■ 0.00000 

DLDB • 

■ 0.00000 

DFGS' 

I 

FGl 73466,09100 ' 

FNl 

73466.09100 " 

■ SFCl' '■ .97996' 

■ DTEHP 

■ ' Oo 00000 ■ ■ ■ 

WC 

■■ 63.00000 

F4 

■ '■ ,03139 



, ' 



Example 2 (Case 2) - The second step in establishing an engine design 
is to lock in or fix the flow matching parameters established in the 
design point case. This is accomplished by setting IDES = 1. 

SENPUT 

ICASE = 2t IDES = 1 , 

s 

At the same time that IDES is set equal to 1, the off-design data 
calculations may also be initiated, however, for simplicity off-design 
performance has been delayed until Example 3. 


Case 2, which follows, is 

also in a 

long form printed output. 





INPUT CASE 

2 

PC 

2 P3 



I 


- 


ALT C. 


HACH 0. 


HPEXT 0 


2D0BS 

0 



— WCOWl .lOOOOOOOE 

-03 

WS 0. 


CV 

.93500000E+00 ZThRUST 

0. 




FN 73A63«69230 

SFC 

.97995 ' FK/s'l 

K-A.951PA 

PTl/PTO 

1.00000 PT2/P71 

] 5.003 17 

£TA2 

,67500 ■ 

" UiK 700.C0032 

N/SRT 

1.00000 

TTQ 

516.6'/C00 

TT2 

1189.90631 TTA 

3060,00000 

PT4/PT2 

.95001 ■” 

■' TT5 a-iSSoOOASB 

PT5/PT4 

.35669 

TT5A 

233Aofc8525 

TT6 

2369.68525 PT6/PT5A 

1.09000 

F6 

.02857 “ 

T8 1623.010<56 

PT8/P3 

5.062^9 

VC 

O.OOOCO 

V3 

3332.92536 H8 

1.73556 

N 

1.00000 ■ 

AA/AAD 1.00000 

A7/A7D 

.99995 

A8/A7 

1.A1270 

A8 

1920.36919 CAPA8 

6.84418 

“ yi 

699.99432"“ 

~ ETAP 0.00000 

ETAC 

■■ 0.00000 

STAO 

O.COCOO ■ 

FN/WIO 

104,95099 WF 

71992,49154 ' 

■ i;g 

719,99224“ 

■“PT7P7 ■ '1.62262 

PT1PT7 

.19675 

PT5T5A 

■■ 1.00000 

PT6T5A 

i.OOClOO PT3PT6 

1.00000 

TTl 

516.67C00 ' 

“ ■ WB O.COOCO 

kBOL'l 

0,00000 

DLDD 

0.00000 

DFGS 

1 'FGl 

73465.69230 ' 

FNl 

73465.69230 

~ SFCl ■ ■■ .97995 

DTENP 

0.00000 

WC 

62.999A9 

FA 

.03139 






Example 3 (Case 3) - The next step in the normal operating sequence is 
to introduce the installation effects of service airbleed, power 
extraction inlet recovery, non-standard day temperatures and exercise 
the variable turbine geometry mode of operation. This is accomplished 
by changing the parameters shown below: 

■ SENPUT 

ICASE=3« IPC=4f 
WB=l • ♦HPEXT=200, < 

TABPTl = .95« .968» »977 * *982 « • 982 » .978. .973. .967 . .96 . . 952 . .943. .934 . .924. .914 ♦ 
.902. .891 . .877. .859. .831 . 

DTEMP=8. ♦ 

A4A4DN=o8. 

S 

Engine performance generated by introducing the above parameters is for 
full and part power without afterburning at sea level static standard 
+8°C atmospheric conditions. It should be noted that this step could 
have been combined with Step 2 (Example 2) above. Output is in the 
long format as follows: 

CO 

I 


■ - INPUT CASE • 

3 

■ PC 

4 P3 



I 


" 


' * ■ 

ALT 

■ 0 , 


hACH' 0. 


HPEXT 

■ .20C0G000E-f-03 ICDBS ' 

■ 0 

* .... , . 

^ - - -- “ * 

“ WCOWl 

" .lOCOOOOOE 

-03 

1<B .100QC000E-^01 

cv 

.9S500000E-^00 ZTHRUST ' 

0 ,' 


■• •• ' 

■" • FN 

65792. 7860A 

SFC 

■ ■ .997-45 

' FN/Wl 

101.96592 

PTl/PTO 

♦95000 

PT2/PT1 

■ 14.54179 

ETA2 

» ii 7 6 8 4 

WIK" 

-682. 57289 

N/SRT 

.96640 

TTO 

533.07000 

TT2 

1209.09110 

TT4 

3060.OCOOO 

PT4/PT2 

.949 JO ■“ 

■ " TT5 

'2A8C.7J097 ‘ 

PT5/PT4 

.35170 

TT5A 

2379.13398 

TT 6 

2379.13398 

PT6/PT5A 

1,00000 

F 6 

.02830 

■ T 8 

i658.3A116 

PTO/PS 

■ 4.61135 

VD 

C.OCOOO 

V 8 

3243,98935 

H 8 

1.67162 

•• • - ^ 

'■ " 1, 00000 

•A4/AAD 

1.00000 ■ 

A7/A7D 

1.01289 

A8/A7 

1,34675 

A0 

1940,63186 

CAPA 8 

8,54264 

■■ ■ ■•VhI 

-645, 24292'“ 

“ etap 

■ , 0.00000 

ETAC 

O.COOOO 

ETAO 

O.COOOO 

FN/ulD 

93.93969 

WF 

65625,19385 

■ ■ ■ WG 

"663.50043 “ 

“PT7P7 

1.82500 

PT1PT7 

.20601 

PT5T5A 

l.OCOOO 

PT6T5A 

“■ 1,00000 

PT 8 PT 6 

l.OOCCO 

“ TTl 

533,07000 '■ 

WB 

l.COOCO 

WBCWl 

. 00155 

DLD 8 

0 . 00000 

DFGS 

I 

F61 

65792.73604 

" ■ FNi 

65792.76604 " 

““ SFCl • 

.99745 ■ 

■ DTtKP 

• 8 . cocao 

IvC 

58.07186 

F4 

,03110 





“ FN 

'60548,63215 

SFC 

.92864 

FN/l/1 

93,63851 

PTl/PTO 

.95000 

’PT2/PT1 

14.54179 

ETA2 

“ ,87684" 

■ WIK 

' 683.57289 

N/SRT 

■ .98640 

TTD 

533.07000 

■■ TT2 

1209.09110 

TT4 

276^., 60000 

PT4/PT2 

■ ;94910 " 

--- TTrs ■ 

2204,15092 

PT5/PT4 

,32019 

•TT5A 

2147.92018 ‘ 

ri 6 

2147,92018 

PT6/PT5A 

i.oooco ■ 

F6 

" ■ .02424 '' 

• T 8 

‘1519.05364 

PT 8 /P 8 

■ ' 4.19820 

',u 

0.00000 

V3 

2997,23055 

MS 

1.607S1 

■ ■ N 

1,00000 ■■■ 

"AA/AAD 

,97707'^ 

A7/A7D 

1,05296 

A6/a7 

1.28304 

A3 

186^,17185 

' ■■ CAPAS 

3 0 43438 ” 

■" ‘ W 1 

645,24292 " ' 

etap 

•" O.COOOO 

ETAC 

■ o.oocoo ■ 

■ ETAO 

0.00000 

FN/WID 

66t>^9S05 

■" ■ WF 

56220.11165 

'■ ■ WG 

66 O 0 S 86 OO 

“PT7P7 " 

1.82500 

• PT1PT7' 

• ■ ,22629 

PT5T5A " 

i.COOOO 

' PT 6 T 8 A 

1.00000 

PTSPT 6 

1.00000 


■■"533‘.07000 ”■ 

v;b 

•” 1.00000 ' 

■ USDWl 

■ • ,00155 

OLDS ■ 

O.OOCOO 

DFGS 

I 

FGl 

60548,63215 

FNI 

60548,63215 ■“ 

■“SFCl' 

.92364 

DTE.MP 

'C.OOCCO 

■ WC 

40.30137 

... . 

■ ■■■ ,02566 

■ “ 

“ 

— . . .. 

V" 



CO 


S'f'i'SA. 19207 


WIK 

nb 

T3 

4ir/A4D 
■• 5TAP 
■“ PT7P7 
~ "WB 
-- SFCl 


603,57209 
■ 1927,963'<.0 
1371.05672 
.95039 
0.00000 
1.84A21 
■" ‘ 1.00000 
• • ■ .85132 


SFC" 
K/SRT 
PT5/PT4 ■ 
PT8/P3 ' 
A7/A7D ■ ■ 
ETAC 
PT1PT7 
WBCWl 
■ DTtMP 


.65132 
.936‘tO 
,26750 
3.7o955 
’1.09993 
0 « 0 0 C 0 0 
.23202 
.00155 
8.00000 ■ 


FN/Ij1 


TfO 533.07000 
TT5A ■ 1904.96503 


84,e5’i25 PTl/PTO 


VO 

A6/A7 

ETAO 

PT5T5A 

OLDB 

WC 


O.COOOO 

1.21471 

0.00003 

l.OCGOO 

0.00000 

22.531S3 


TT2 1209.09110 


.95000 PT2/PT1 


14.54179 


ETA2' 


TT4 2509.20000 ’ PT4/PT2 


TT6 1904.96503 PT6/PT5A 


A3 

FN/IJID 

PT6T5A 

DF6S 

F4 


2721.33205 

1776.32303 

73.22027 

1.00000 

I 

.02084 


:*18 

CAPA3 


l.OCOOO 
1, 53030' 
8.34980 


KF 46640.76335 


PT3PT6 


l.OCfOOO 


FGl 54754.19207 


; F6; 

W1 • 

■ WG 
TTi 


.’876 3 4" 

■ “.94910Z? 

.02011 -" 
1,00000 ■" 
645.24292 “ ’ 
658,21380 r' 
533.0700C , 


FNl 54754.19207 


FN 

48050.80353 

SFC ■“ 

.76787 

FN/Wl 

74,46932 

PTl/PTO 

- .95000 

PT2/PT1 

14.54179 

ETA2 

,876847 

WIK ■■ 

■ 688,57269 

N/SRT 

“.98640 

TTO 

533.070'’0 

“ TT2 

1209.09110 

TT4 

2233,60000 

PT4/PT2 

,94910 -’ 

TT5 

1647.72963 

PT5/PT4 ■ ■ 

,24753 

TT5A 

1644.69101 ■ 

■ " TT6' 

1644.69101 

PT6/PT5A 

’ 1,00000 

■ F6 

,01591 

„ _ . ^ g - 

'1217.47082 ■ 

PT8/P8 

3.24556 

VO • 

0.00090 

V8“ 

2393.08822 

M3 

1.42496 

- ■- N 

1,00000 ■■ 

A4/A4D • 

-■ ■ .91343 

’ A7/A7D ’ 

1,13214 

A6/A7 ’ 

1.1 4248 

A8 

1682,91550 

CAPA3’ 

8.37820 

Wl 

645.24292 ’” 

— ETAP 

O.OOOOO 

ETAC ’” 

O.OOOOO 

ETAQ 

0.CCC09 

FN/WIO 

63,64401 

WF 

36896.60453 

WG 

" 655.50768 ‘ 

•- PT7P7 ■ 

1.64543 

■ PT1PT7 ’ ‘ 

,29271 

PT5T5A 

l.CCOOO 

PT6T5A 

l.COOOC 

PT8PT6 

1,00000 ‘ 

7T1 

5-33.07000 • 

yB“ 

1,00000 

W80K1 ■ ■ 

.00155 ■ 

' DLD3 

o.occoo 

“ DFGS ' 

I ' 

FGl 

48050,80358 

■■ FNl 

43050.80358' ' 

SFCl" 

’.76787 

DTEKP • ’ 

0.COOOO ’ 

WC ’ 

4.76189 

F4 

.01603 




- 1,. 

FN 

39150.25499 

SFC 

.67692 

‘ FN/Wl 

60,67522 

■ PTI/I'Td' 

.95000 

PT2/PT1 

'■ 14.54179 

ETA2 

.87684 

- WIK 

- 683,57289 

N/SRT 

.93640 

TTO 

533. 07000 

TT2 

1209.09110 

■ TT4 

1958,40000 

PT4/PT2 

’■’ . .94910 ’ 

' ’T75 

1352.13672 

PT5/PT4 

.19392 

TfSA 

1352.188’!’2 

’ ' TTO 

1352,16372 

PT6/PT5A 

1.00000 

F6 

.01143 ■ 

T3 

1055.37005 

PT8/P8 

2.34255 

VG 

’ ’O.OOOOO 

V8 

1962.64370 

M8 

1.24635 

~ ' N 

1.00000 

• A4/A40 

.36243 

' A7/A’?D 

1.36221 

AS/A7 

“ 1,04312 ” 

A3 

1574.62824 

■ CAPAS 

• 3.33325 

Wl' 

’’ 645.24292 

"• ETAP 

0.00000 

'■ ETAC 

0.00000 

■ ETAO 

■ 0.00000 

■’’FN/WIO ■ 

’■ 55.92894 ’ 

' ■ WF 

26501,75746 ' 

WG 

’ 652.61594" 

■ PT7P7 

- ” 1.86511 

PT1PT7 ■ ■ 

■ .37364 

PT5T5A 

1.00009 

’ PT6T3A ’ 

1.00000 

'PT6PT6 

1.00000 

TTI 

533.07000 ’■ 

WB 

1,00000 

■ WBOWl 

,00155 

DLDE 

O.OOOOO 

■" OFGS 

1 

• FGl 

39150.25499 

’■ ' FNl 

3915C. 25499 ’’ 

SFCl 

.67o92 

DTcKP 

8.00CCO 

WC 

4.7C189 ’ 

" ■ F4 

.01143 



.... 

. .. , 1 


FN 

23105.72120 

SFC 

.60227 

FN/Wl 

■ 43.55836 

PTl/PTQ 

,95000 

PT2/PT1 

14.54179 

ETA2 

.87684 

— ' ■ WIK 

688,57289 

N/SRT 

' .93640 ' 

‘ TTO 

533.07000 

‘ TT2 

1209.C9110 

TT4 

1699.64340 

PT4/P72 

”” -.94910 

" TT5 

1066.41963 

PT5/PT4 

.13720 

TT5A 

1066,41963 

■ ■ TT6 

1006.41963 

PT6/PT5A 

1.00000 

F6 

■ .00720 ■ 

ta 

907,.33307 

PT8/PB 

1,79594 

VO 

O.OCCOO 

■ V0 

1414,67195 

M8 

.96490 

N 

" i.cooeo' 

A4/A40 

.80021 

A7/A7D 

’ 1.70461 

A8/A7 

1, OOOOO 

AO 

1466.13749 

CAPA8 

10.49749 

Wl 

645.24292 1 

• ETA? 

0 .oocco 

ETAC 

0.03000 

ETAQ 

0.00000 

FN/Win 

40.15103 

K? 

16927,14210 

WG 

649.95220 

PT7P7 • 

* 1,86511 

PT1PT7 

.52309 

PT5T5A 

’ 1,00000 

PT6T5A 

1.00000 

" PT8PT6 

1.00000 

TTI 

533.07000 

WB ‘ 

1,00000 

■ waowi 

" .00155 

DLDB 

O.C'CCOO 

OFGS 

I 

FGl 

28105.72120 

FNl 

28105.72120 

SFCl 

.60227 

DTEHP 

8.00000 

WC 

4.76139 

F4 

.00730 




.. 

FM 

27354.63604 

SFC 

.58923 

FN/Wl 

42.36349 

PTl/PTO 

.95000 

PT2/PT1 

14.45092 

ETA2 

.5;710 

WIK 

689.07525 

N/SRT 

.93640 

TTO 

533,07000 

TT2 

1206.74075 

7T4 

1674,64340 

PT4/PT2 

■ ■ .94349 

TT5 

1045.50730 

PT5/PT4 

.13533 

TT5A 

1045.50730 

TT6 

1045.50730 

PT6/PT5A 

1.00000 

F6 

• ‘ ,00694 

’’ • T6 

895.14919 

PT3/P6 

1,76219 

VD 

c.ccooo 

V8 

1376.34752 

' !I8 

,94399 

■ N 

l.OOOCO 

A4/A4D- 

,30021 

A7/A7D 

1.726S4 

A8/A7 

l.oocno 

A3 

1453,00560 

CAPAS 

10.64295 

Wl’ 

645,71367 

ETAP 

’ 0.00000 

ETAC 

0,00000 

ETAO 

' 0.00000 

FN/Wl 0 

39.07032 

WF 

16113.21302 

' "■ WG 

650.19739 

■ PT7P7 

1.86511 

PT1PT7 

’ .53910 

PT5T5A 

1,00000 

PT6T5A 

I. OOOOO 

PT3PT6 

l.OCOOO ' 

■ TTI’ 

‘ 533.07000 

WB 

1,00000 

WBDWl 

.00155 

DLDB 

G.OCOCO 

OFGS 

1 

FGl 

27354,68604 

FNl 

27354.60604 

SFCl 

.53923 

OTEMP 

3.00000 

WC 

4.76149 

r*'f 

.00694 






ORIGINAL PAGE Ib 
OF POOR QUATjITY 



B-3 


FN 

24122.20640 

SFC • 

.58308 

FN/Wl 

38.53141 

PTl/PTQ 

.95000 

“PT2/PT1 

• 13.75074 

ETA2 

.83032 

"" ■ WIK 

668.08044 

• N/SRT • 

.96174 

.. 

533.07000 " 

' TT2 

1187.27637 

TT4 

1611.57693 

PT4/PT2 

.94738 

TT5 

995.72072 

PT5/PT4 

.13162 

TT5A 

995.72072 

TT6 

995.72072 

PT6/PT5A 

liOOOOO 

F6 

.00625 

" ■ T8 

870.71162 

PT8/P6 

1.62894 

VO 

0.00000 ■■ 

V3 

1252,77130 

M8 

.67108 

N 

■■■■•• ■',9 7500 

■ A^t/A^D 

■ ■ • .60021 

A7/A70 

1.78776 

A8/A7 

i.oocoo ■ 

A8 

1433.13897 

CAPA3 

11. 01866 

“■ Wl' 

“"626.04000 . 

r~ ETAP 

O.GOOOO 

fcTAC 

0,00000 

ETAO 

o.occoo 

FN/WID 

34.46029 

WF 

14065.27010 

WG ■ 

629.95327 

— PT7P7 ■■ 

■■ • 1.86511 

PT1PT7 

.53320 

PT5T5A 

1.00C09 

PT6T5A 

1,00000 

PT8PT6 

l.OOOCO 

■ ■■ TTl 

533.07000 

ViB 

1.00000 

WBQWl 

.00160 

DLOB 

o.ooboo 

DPGS 

I 

FGl 

24122.20640 

FNl 

24122.20640 

SFCl " 

.58308 

OTfcKP 

8.00000 

WC 

4.76189 

F4 

.00625 




"* * 

■ " FN 

21052.53062 

SFC 

.57973 

FN/Wl 

34.71029 

PTl/PTO 

.95000 

PT2/PT1 

13.07423 

■ ■■ ETA2 

.88364 

WIK 

647.25107 

N/SRT 

.93708 

TTQ 

533.07000 

, TT2 

1167.78453 

TT4 

1550.59924 

PT4/PT2 

.9462'/ ' 

- •• -TTS 

946,49931 

PT5/PT4 ■ 

■ .12814 

'* TT5A 

948,49931 

•TT6 

943.49931 

PT6/PT5A 

l.OOCOO 

F6 

.00560 

T8 

846.71615 

PT8/P8 

■ 1.50607 

VO 

0.00000 ■ 

■ V3 

1129.32693 

MS 

.79573 

N 

.95000 ' 

■ AA/AAD ■ 

.60021 

A7/A70 ■ 

1.87223 

A8/A7 

l.COOOO “ 

A8 

l419,229-/3 

CAPAS 

11.50735 

Wl 

606.52136 

ETAP 

0.00000 

ETAC 

• 0.00000 

ETAO 

0.00000 

FN/WID 

30.07504 

WF 

12204.73660 ' 

WG 

609.91716 

“ PT7P7 

1.86511 

PT1PT7 

,63078 

PT5T5A 

1.00000 

PT6T5A 

l.COOOO 

PT6PT6 

l.OOCOO 

TTl 

533.07000 

WB ■ 

1.00000 

WBOWi 

.00165 

DLD3 

O.OOOCO 

DF6S 

1 

FGl 

21052.53082 

FNl 

21052,53082 

“■•SFCl" 

■ ■ .57973 

DTEMP 

6 .00000 

WC ■ 

• 4.76189 

F4 

■_ .00560 



. . 

... 

" • FN 

16921.46838 

SFC 

“■ .53768 

■ FN/Wl 

29.99917 ' 

PTl/PTO 

.95000 

PT2/PT1 

11,84306 

ETA2 

,88793 

WIX 

601.94308 

N/SRT 

.83776 

TTO 

533.C7000 

TT2 

1132.60961 

TT4 

1471.34021 

PT4/PT2 

" .94512 

’■ TT5 

397.71741 

PT5/PT4 

.12916 

TT5A 

697.71741 

TT6 

897,71741 

PT6/PT5A 

l.OOOCO 

F6 

■■“ ■ .00491 

T8 

■ 821.66869 

PT6/P8 

1.37398 

VO 

O.OOOCO 

V3 

976.83984 

H3 

.69859 

N ■ 

,900'00 

■ 'AA/A40 ■ 

.80021 

■ A7/A7D 

■■• 1.95438 

A0/A7 

1.00000 

A3 

1398,30199 

CAPAS' 

11,99664 ' 

Wl 

“' 564,06448 

~ • ETAP 

0*00000 

ETAC 

0.00000 

■ ETAO 

0.00000 

FN/WIO 

24.17353 ■ 

WF 

9944.40820 

WG 

■ 566.63172 

- PT7P7 

1.66511 

PT1PT7 

.69142 

PT5T5A 

l.OOCOO “ 

PT6T5A 

■ 1. 00000 

PT8PT6 

1.00000 

TTl 

533.07000 

“ •■ WB 

1.00000 

• WBCWl 

"■■ .00177 

DLDB 

“ 0,00000 " 

’ ■ OFGS ■ 

I 

■■ FGl 

16921.46838 

FNl 

16921,46388 ‘ 

" ■ SFCl • 

.58768 

DTEMP 

■ e.oocoo 

WC 

4.76189 ■■ 

■ F4 

■ ■ .00491“ 


■ 



FN 

13373.72161 

SFC 

■ ■■ .61379 

■ FN/Wl 

26.04554 ■ 

PTl/PTQ 

.95000 

PT2/PT1 

10.54602 

■ ■' ETA2 

.886 51’ 

WIK 

547.95588 

N/SRT 

” .83844 

TTO 

533.07000 “ 

• TT2 

1096.11391 

TT4 

1406.06951 

PT4/PT2 

.94444 ■ 

'• • TT5 

865.35002 

•PT5/PT4 

.13550 

TT5A 

86S.35C02 

■ TT6 • 

665.35002' 

PT6/PT5A' 

1.00000 ■■ 

■ F6 

" ■ “ .00445 ■ 

T8 

• 806.83532 ' 

PT8/P8 

■ 1.28211 

VO ■ 

0.00000 

■ ■ Vfl 

848,63466. 

H8 

.61239 ■ 

N 

■.85000' ' 

A4/A4D 

.60021 

A7/A70 

2.00525 

• A6/A7 

i.cccno 

A3 

1385.78292 ' 

CAPA8 

12.33607' 

' Wl 

513.47454' 

ETAP ’ 

O.OOOCO 

ETAC 

0.00000 

ETAO 

O.COOOO 

FN/WIO 

* 19.10532 ■ 

WF 

8208, 65120 

■" “' WG 

'515,.75917 

■■ PT7P7 

■ ■■ 1.86511 

PT1PT7 

.74097 

PT5T5A 

1.00000 ■ 

PT6T5A 

' . 1.00000 

PT8PT6 

1.00000 

■ ““'TTl 

■ ■ 533,07000 O O 

7 ■ W8 

1.00000 

W3CW1 

.00195 

DLDB 

■ o.ocooo 

DFGS 

1 

■ FGl 

13373.72161 ' 

FNl 

13373.72161 tO 

—“SFCl ■ 

■•■ .61379 

■ DTEMP 

■ 8.00000 

WC 

“ 4,76189 

F4 

.00445' 




. _.^g 

FN 

10496.02714 

SFC- 

• ■ .65817 

• FN/Wl 

'• 22.74348 ’ 

PTl/PTO 

.95000 

PT2/P71 

9.29191 

■ ETA2 

'.58212 S 

^ WIK 

492.48698 

N/SRT 

.78912 

TTO 

533.07000 

TT2 

1058.74038 

TT4 

1351.32592 

PT4/PT2 ■ 

■ g 94416 " ^ 

. TT5 

“ 844.33375 

PT5/PT4 

■ .14566 

TT5A • 

• 844.38375 

■ TT6 

344,36375 

PT6/PT5A 

• l.OOOCO ■ 

F6 ■ 

'.00417jO 

1 Ts - 

““799.71876 

PT8/P8 ' 

1.21417 

VO 

0.00000 "■ 

V8 

741.41572 

MS ■ 

.53731 “ 

■ N ■ 

" ■' ■ '.eoooocj ^ 

rA4/A40 ■ 

.00021 

A7/A7D 

2.04027 

■ A8/A7 ■ 

1.00000 '• 

■“ • A8 

1379,85954 ‘ 

■ ■ CAPAS 

' 12.57155 “■ 

• Wl 

“ 461.49615P> 

-■■ “ETAP ■ 

•— 0.00000 

ETAC 

“ ■ 0.00000 

‘ ETAO ■ 

" 0.00000 ■■ 

'FN/WID ■ 

14.99432 

WF 

6908,12951 “■ 

WG'" 

" ' 463i41924£2 fil 

- PT7PV - 

1.86511 

PT1PT7 

• ■ .73242 

PT5T5A 

1,00000 ‘ 

PT6T5A, 

' 1.00000 

PT8PT6 

1.00000 ■ 

TTl 

“*533.07000*=3 • 

WB 

• 1.00000 

WBQWl 

.00217 

DLDB 

0,00000 

• DFGS 

I 

FGl 

X0496. 02714 " 

FNl 

10496.02714^ K 

■“" SFCl" 

■— .65817 

OTEMP 

' 6.00000 

WC 

• 4.76189 

F4 ■ 

■ .00417 

— • ' 

• 

“ * •“* * 

’**“**■ * ‘ 



FN 

■ 7603.69'J8^ 

src 

.75159 

FN/Wl 

• 19.21616 

“ ’ ■ WiK' 

■ A33^66057 

N/SRT 

■ .73930 

•iTO 

S33.C70-:'0 

- TX5' 

■ ■■ 333.39027 

PT5/PTA 

.15925 

TT5A 

' 833,39027 

T8 

■' 50i,AS116 

■ PT8/P8 

■ 1»1A608 

VO 

0.00000 

"AA/A^D ' 

■■ .80021 

• A7/A7D • 

2.13192 

A6/A7 

■ l.COOOO 

• - ETAP' 

0.00000 

ETAC 

■ 0.00000 

ETAO 

0.00000 

- PX7P7 

1.06511 

PT1PT7 

■ .32797 

PT3T5A 

■■ 1.000^0 

WB 

loOOOOO 

■ WBDWl 

.00296 

DLDB 

”0.00000 

SFCl 

- .7315A 

DTEKP 

■ 8.00000 

WC 

•• 9.76189 


" “NCrZLE PRESSUf^c R^TID LESS THAN 1*005 


PTl/PTO 

.95000 

PT2/PT1" 

8.05708 

■ ■ ETA2 

.66799 

TT2 

1022.50712 

TT9 

i206.7S936 

PT9/PT2' 

• 9 A H 1 2 

TT6 

833.39027 

PT6/PT5A 

1.00000 


.00902 

V8 

626.70398 

H6 

,95372 

”■ ■ H' 

“ ■ .76000 

A8 

1381.25936 

CAPA6 

■ 13.11916 

~ W1 

■■ 906,37151 

Fk}/W10 

11,15557 

WF 

5868,68838 

‘ WG 

■ '908.00573 

PT6T5A 

1.00000 

■ PT6PT6 

1.00000 

-TTl 

"■“5 33.07000 

DFGS 

I 

■ FGl 

■■7803,89989 

■ FNl 

“7806,39939 

F9 

.00902 



“ “ • •"* 

■ ' . X.... 


CO 



Example 4 (Case 4} - To generate engine performance with afterburning 
requires changing the afterburner pressure losses and the power control 
parameter (I PC), as shown below. 


SENPUT 

I CASE=4* IPC-1 ♦PT8PT6=»98»PT 65A0=*95» 
$ 


DO 

I 

o 


The values introduced for the pressure losses, above, are typical of 
current afterburners. It was not necessary in this case to introduce 
afterburner combustion efficiency or the maximum afterburner temperature 
since these were input in Example 1. It was possible to do this in Case 
since these parameters are used in the cycle match calculation only if 
the power control calls for an afterburning case. Pressures, on the 
other hand, could not be input in Case 1 since they are used in the cal- 
culation reguardless of whether or not the power control calls for an 
afterburning case. 

Shown below is -afterburning engine performance at sea level static 
standard +8°C atmospheric conditions from maximum afterburning to maxi- 
mum non-afterburning, with the output in the long form format. 


INPUT - CASE 

• • 4 

- PC 

1 P3 

HPE>!T 


1 





ALT - 0. - • •• 


KACH ■ 0. 

OCOOGE+01 

,20000000Ev03 IDD8S 

0 ■ 



WCOWl • .iOOOODOOE 

-03 ■ • 

WB .IOC 

CV 

■ ,9850C-OO0E'*C0 ' 2THRUST 

0, 

■ . . 

■■■ FNS7^48.A5146' 

' ■■ SFC 

1.75800 

' ■ FN/Wl ■ 

135.52795 

PTl/PTD 

.95000 

PT2/PT1 

14.54179 

ETA2 

. S768A 

— - WIK 68e.57a89 

■ N/SRT ■ 

■ .98640 

■■ TTO 

533.07000 ■ 

TT2 

1209.09110 

TT4' 

3060.000C0 

PT4/PT2 

,94910 

TT5” 2480.70097 

PT5/PT4 

.35170 

■ ' TT.5A 

2379.13398 ' 

.... . 

4060,00000' 

PT6/PT5A ■ 

“ ■ .91738 ■ 

’■ F6"' 

~ ■ ■■.06629— 

T8 "’SOaS. 52731 ' 

PT3/P3 

— • 4.14576 

VO • 

O.COCOO 

• V8 

4158.12654 ' 

,K8 

■1,61971-’ ■ 

... 

— i.coooo'- 

■'•A4/A4D ■■■ a.OUCOO ' 

•A7/A7D 

" - 1.52614 

■ A8/A7 

1.22880 

A3 ■ 

2567.20614 

■ CAPAS 

12.66504 

wi ■ 

"■ 64 5.2429.2"' ■ 

- ETAP" 0.00000 

ETAC 

■■ '0. 00000' 

‘ • ETAO 

-• O.COOOO 

FN/WIO 

124.92636 

WF153734. 79550 

" - ■ WG ■ 

■”688.01331 " 

- PX7P7 -1,30000 ■ 

- PT1PT7 ■ 

- — .22915 ■ 

PTDT5A 

• l.CCOOO 

■'■PT6Tf A 

" "i 91738’: 

PT8PT6 

.98000 

’TTl' 

■”533i07000 "■ 

WB ■ 1.00000 ■ 

" hBOWl ■ 

.00155 

— - DLDB' 

0.00000 

■ DFGS ■■ 

.. J 

•PGI 

87448.45146 

FNl 

87446. 45.146""" 

— ~SFC1 - “ 1.758C0 

‘ DTEMP- 

8.00000 

wc ■ 

58.07186 

-- - - P4- 

'.03110 

' 


. 

. ... 

• 80355.60331 ■ 

SFC 

1.52572 

■ FN/Wl 

" 124.53543 

■■ PTl/PTO 

.95000 

PT2/PT1 

■ 14.54179 

ETA2 

.87684 ' 

WIK ■ 683.57289 

N/SHT 

.98640 

TTO 

533.07000 

TT2 

1209.09110 ■ 

TT4 

3060,00000 ' 

PT4/PT2' 

■■ .94910 

TT5~" 2460.70097 

PT5/PT4 

' ■ .35170 

■■ TT5A 

2379.13398 

■■ ■ TTo 

3516.37799 

PT6/PT5A 

.92826 " 

F6- 

.05286 

~ ■ T8 2591.53855 

“Te/pe 

■ 4,19489 

\/o ■ 

O.OGCOO 

V8 

3669.58325 

rt8 

1.62154 

N- 

l.OOOOO ' 

•A4/A4D- ■ 1.00000 

• A7/A70 

1,33599 

A8/A7 

1.32746 

A8’ 

2386,36657 

’■ CAPA3 

11,46768 ■ 

Wl 

■■ 645.24292 

“ ’ETAP ~ 0.00000 

ETAC 

• 0.00000 

ETAQ 

C. 00000 

FN/WIO 

114.79372 

WF122599. 76960 

WG 

“ 679.35127 ■ 

~ PT7P7 1.81019 

PT1PT7 

.22647 

PT5TDA 

■ 1,00000 

PT6T5A 

,92826 

PT8PT6 

,98000 ‘ 

• ■ TTl 

533.07000 

■ ixB ■ 1.00000 

VjBOtgi 

'■ .00156 

DL03 

0.00000 

DFGS 

I 

FGi 

80355.60331 ’ 

FUi 

80355.60331 

■■ SFCl " ■ 1.52572 

DTch? 

8.00000 

WC 

5 So 071 86 

F4 

.03110 






o o 
Si 

gl 

KK- ' 



“ ■■ FN 

729«1^5iV69 

SFC 

1.28873 

FN/Wl 

113.10751 

PT1/P70 

.95000 

PT2/PT1- 14.54179 

■ ETA2 

' .876S4 

WIK 

63D»$72S9 

ft/'SRT 

.93640 

TTO 

533.07000 

TT2 

1269,09110 

.TT4 

3060.00000 

PT4/PT2 

.94910 ' 

TT5 

2480.70097 

PT5/P74 

.35170 

TT5A 

2379.13398 

TT6 

2972.75599 

P7o/PT5A 

.93913 

‘ F6 

■ .04055 " 

T3' 

2155.09933 

PT3/P8 

4.24403 

va 

0,80009 

V3 

3556.06396 

(13 

1.62375 

N 

l.OOOOO 

“A4/A4D ' 

l.OOOOO 

A7/A70 

‘ 1.24488 

■ A8/A7 

1,33100 

A8 

2 190,6330 3 

CAPA8 

10.25210 

■ ■ W1 ' 

645.24292 

” ■ ETAP " 

OiOOOOO 

ETAC 

0.00000 

ETAG 

0. 00000 

FN/WIO 

104.25974 

WF 

94053.98272 

WG 

671.40953 ' 

~ PT7P7 

■ 1.8i424 

PTIP77 

• 223 84 

PT5T5A ■ 

l.COOOC 

PT6T5A 

■ .93913 

PT8PT6 

.96000 " 

TTl 

533,07000 ■ 

WB 

■ 1.00000 

WBQVil 

.00155 

OLOB 

o.oocoo 

OFOS 

1 

FGl 

72981.31769 

■ FMl 

72981.81769 

■ SFCl • 

1.20873 

■ DTENP 

■ 8.00000 

wc 

5S.C7136 

F4 

a 0 3 J. X 0 





■ • FN 

65194.59372 

SFC 

1.04169 

■ FN/Wi 

101.C3884 

PTl/PTO 

.95050 

PT2/PT1 

14.54179 

ETA2 

.87684 ' 

■ WIK 

668.57269 

N/.SRT 

,9’3640 

■ ■ TTO 

533,07000 

TT2 

1209.09110 

•fT4 

3060.00000 

PT47PT2 

,94910 

■ TT5 

2480.70097 

PT5/PT4 

.65170 

"TT5A 

23T9.13398' 

TT6 

2429,13.393 

PT6/PT5A 

,95000 

F6 

.02928 

T8 ■ 

1726.20456 

PT6/PS 

<». 29316 

VO 

0.00000 

ve 

3211.41503 

M8 

1.62401 

N 

1.00000 

• AA/AAD 

l.COCOO 

A7/A7D 

1.10039 

A3/A7 

1. 62316 

AS 

1977.46563 

■ CAPA8 

8.99133 

■ . W1 ' 

645,24292 

ETAP 

'• 0.00800 

ETAC 

0.00000 

ETAO 

o.coooc 

■ FN/WID 

93.13513 

WF 

67912,63519 

• WG 

664.13663"' 

■ PT7P7 

1.62500 

PT1PT7 

.22128 

PT5T5A 

l.COOOO 

PT6T8A 

.95000 

PT8PT6 

.98000 

' TTl 

533.07000 

W8 ■ 

l.OOOOO 

WBOMl 

.00155 

OLOe 

o.coooo 

OFCS 

I 

FGl 

65194.59372 

" FNl 

65194.59372 ”” 

SFCl 

1.04169 

OTtNP 

e . 00000 

WC 

53,07186 

F4 

,03110 



. 


FN 

64493.53907 

SFC 

1.01755 

FM/U'l 

99. 9.5? 34 

FTl/PTO 

.95000 

PT2/PT1 

14.54179 

' ETA2 

.87634“ 

WIK ■ 

653 .572o9 

H/SRT 

.98540 

TIG 

533.C7000 

TT2 

1269,09110 

TT4 

3060.00000 

PT4/PT2 

.94910 

TT5 

24S0. 70097 

PT5/PT4 

.35170 

• TT5« 

2379.13398 

TT6 

2379.13393 

PT6/P73A 

.95000 

F6 ' 

.02830' ■ 

78 

1687.50930 

PT3/FS ■ 

4.29316 ' 

VO 

0,00000 

V3 

3179,92897 

1-18 

1.62437 

N 

"■ l.OOOOO”"' 

■ 'AA/A4D 

1.00000 

A7/A70 

1.03796 

A8/A7 

1.32616 

AS 

1957.64151 

CAPA3 

8.86317 

■■ Wi 

645.24292 ■■ 

■" ETAP 

0.00000 

ETAC 

0.00000 

ETAG 

'O.COOOO 

FN/WIO 

92,133o3 

WF 

65625.19385 

' WG" 

663.50043 

-- PT7P7 

■ 1.325U0 

PT1FT7 

.22123 

PT5T5A 

1.C60C0 

PT6T5A 

.95000 

■ P78PT& 

.98000 

■■ TTl 

533.07000 ■ 

W3 

l.OOOOO 

U30v;i 

.00155 

DLOB 

0. CC 000 

DFGS 

1 

FGl 

64493.53907" 

■*" FNl 

64493.53907" 

— SFCl 

1.01755 

DTEMP 

8.00000 

WC 

58.07186 

F4 

■ .03110 







Example 5 (.Case 5) - This case is similar to Case-'4 except that this 
case is for a Mach number of 2.62 at an altitude of 65000 ft. with a 
fixed geometry turbine as determined by the input’ parameters shown 
bel ow: 

SENPUT 

ICASE=5« IPC=3» IPUNCH=2tA4A4DN=l ,0 tALT=65000, »EMACH=2*62 » 

s 

Engine performance for this case is generated from maximum afterburning 
through maximum non-afterburning to minimum non-afterburning power as 
determined by the program. Output is in the short form format. 

INPUT CASE 5 PC ■ 3 P3 ,500000006+01 

ALf .65C00O00E+05 HACH .262COOCOE+01 ' HPEXT .20COOOOOE+C3 10D3S "O'’ 

WCOWl .lOOOOOOOE-03 UB . lOOOCOCOE^Ol ‘ CV ,03500000£+C0 ZTHRUST ” — 0, ' 


FN 

37673,7338 

FG 

65271,4233 

' FO 

27547.6896 

SFC ' 

1.96014 

fef?35o5,G690 

Wl 343,7901' 

WG' " 

364.362V 

FNl 

37673,7338 

FGl 

65271,4233 

N 

' 1 coooo 

■ SFCl ■ 

1.96014 

WIK 439.3817 

V6 5863,4462 

"CAPA3 ■ 

54.0292 

7TA 

3060,00000 

A4/A40 

1.00000 

■ • T71 

954.22444 

" ’ . 








rH 

31762.38^1 

FG 

59390.0737 

FD 

27597.6696 

'' SFC 

1.79135 

■ ■" " 'riF56951.9319 

■ Ui 343-7901 

" ' ' WG ■ 

359.6561' 

FIJI 

31792,3341 

FGl 

5939C.0737 

N 

l.OCOOO 

■’SFCl ■■ 

1.79135 

"■ WIK 439.3617' 

' ■ Ve 5409.5405 

' "CAPAS ■ 

47.9313 

TT.% 

3060.00000 

A4/A4D 

l.OOOCO 

■ TTl 

954.22444 

■ 






FN 

25737,6530 

FG 

53335.3476 

FD 

27597.6896 

SFC ■ 

1.61290 

WF41509.3103 

Wl 343.7901 

WG 

'355.3542 

FNi 

25737.0580 

FGl 

53335,3476 

N 

1.00900 

SFCl 

1.61280 

'WIK 439.3817 

■ V8 '4916.8568 

CAPAS 

41.762S 

TTA( 

3060.00COO 

A4/A4D 

1.00000 

TTl 

954.22444 





' 1. . 


FN 

19367,3264 

FG 

46965.0160 

• " FO 

27597,6896 

■ ■' SFC 

1,41545 

■ ■" WF27413.5730 

■ ""'Ml 343.7901 

WG' 

351.4272 

FNl 

19367,3264 

FGl 

46965.0160 

N 

1.00000 

' SFCl' 

1,41545 

WIK 439o3317 

“V8 4377,9726 

" ■ CAPAS " ' 

35,5354 

TTs 

3060,00000 

A4/A40 

1.00000 

TTl 

954,22444 







FN 

f 

18810.2355 

FG 

46407.9250 

' ■ FD 

27597.6896 

SFC 

1.39437 

■ ■ ■ WF26223.S612 

'“ ■" Wl 343.7901 

WG'" 

351,0970 

FNl 

10810,2355 

FGl 

46407.9250 

"■ N 

1,00000 

SFCl 

1.39437 ■ 

■ WIK 439.381? 

■ ve 4330,1103 ■ 

"' CAPA3"" ■ 

34,9767 

TTA 

3060.00000 

A4/A4D 

1.00000 

TTl 

954.22444 


_ 

I'll. 






FN 

18254.0938 

FG. 

45660,4162 

• ■ FO 

27605,6224 

- -SFC" 

■' 1,38795 

WF25336.'95e0 

' Wl 343,8877 

WG""' 

350.9462 

FNl 

16254.6936 

FCl 

45860,4162 

N 

1.00000 

■ ' SFCl 

1.38795 ■ 

-"W1K“ 439.5064 

ve 4260.8594 

"■"CAPA8 " 

34,7560 

TTA 

3010.00000 

A4/A4D 

1.00000 

TTl 

954,22494 

. 






FN 

15351,5265 

FG 

43492,5910 

■ FD 

27641.0646 

■ ■ ‘ SFC" 

■ 1,36532 

- WF21642.3091 

Wl '344'. 3304 

WG"" 

'350.3597, 

FNl 

15851,5265 

- FGl 

43492,5910 

N 

1.00000 

■■ ■ SFCl ■ ' 

l;36532 

WIK 440,0722 

V8 4066.6150 ' 

"•■"CAPAS 

■33.9167" 

TT-j 

2604,43746 

A4/A40 

1,00000 

TTi 

954,22444 




rv 



FN 


rG 

41023.0049 

FO' 

'276'!53,6812 

• ' src 

3 , 34914 

WF17996,6543' 

Wl 344.8613 

■' ■ "WG" 

349,3749 

•FM 

13339,3236 

FGl 

41023,0049 

U 

l.OCOOO 

■' SFCl ■ 

1.34014 

' Wl!< 440,750? 

"V8 3841.0030 

CAPA6__ 

33,1212^ _ 

T7<> 

2603,44037 

A4/A4D 

1.90000 

■ TTl 

959,22444 

• 



. ^ . • * 





FN 

11339.1437 

FC> 

39057.3279 

FO 

27718.1842 

SFC 

1.34500 

WF1S251.1399 

W1 345,2911 

■ ■ WG - 

FNl 

11339,1437 

F61 

39057.3279 

N 

1.00000 

OFCl 

•• 1.34500 

WIK 

441.3001 

' ■ VS 3660.4481 

CAFA3 ■ 

TT4 

2453.44037 

A4/44D 

1.00000 

■ • TTl 

554,22444 

' 






- FN 

1C649.8952 

FG 

383SO.3S0O 

• FO 

2773C.4847 

S!=C 

1. 34582 

yF14332.3C41 

W1 345.4443 

WG' 

FNl 

10649.6952 

FGl 

36380.3300 

■ - N 

1.0 0000 

■■■ 'SFCl 

1.34582 

■■ WIK 

441.4959 

■ V8 3598.0562 ' 

CAPA8 " 

;tt4 

2403,44037 

A4/A4D 

l.OOCOO 

■ TTl 

954.22444 







' r-N 

9589,0507 

FG 

37333.6314' 

F3 

27749.5727 

"• SFC 

■ 1,35793 

■■ • ' WF12963.&959 

rWl 345.6821 

■■ WG" 

FNl 

9569.0567 

FGl 

37333.6314 

N 

l.OOOOO 

SFCl 

1.35193 ■ 

■ 'UIK 

441.7998 

V8 5501. 6275 ■ 

' “■ CAPA3 ■ 

TT4 

2328.99281 

A4/A4D 

l.COOOO 

- "TTl 

954.22444 







FN 

3116,6698 

FG 

35895. 599C 

■ FD 

27776.7292 

SFC 

1. 36960 

WF11119.6021 

■ ■■ W1 346.0204 

■ ■ ■ WG 

FNl 

8118.6698 

FGl 

35895.5990 

■ ■" N 

1.00000 

SFCl 

1.36960 

■■■ WIK 

442,2321 

" ■" ■ Ve 3368.1741 

' ' CAPA8 ■■ 

TT4 

2228.99231 

A4/A4D 

1.00000 

TTl 

954,22444 







FN 

6755.3663 

FG 

34561,7014 

FO 

276C6.2352 

SFC 

1.40251 

WF 

9474,^ 'f47 

■■ W1 346.3892 

'WG' 

FNl 

6755.3653 

' FGl 

34561.701'* 

N 

l.OOOOC 

■ SFCl 

1.40251 ■' 

WIK 

442.7025 

■ Vo 3243.6328 ' 

■ CAPA8 ■ 

TT4 

2140,12151 

A4/A40 

1.00000 

TTl 

954.22444 







FN 

5762.7769 

FG 

33598.7224 

rD 

27635.9455 

■■ SFC 

1.43293 

Wf 

8372.9401 

W1 346.7581 

WG " 

FNl 

5762.7769 

FGl 

33598.7224 

■ ■ (1 

l.COOOO 

SFCl 

■ 1.45293 

WIK 

443.1749 • 

V8 3152.8815 

CAPAS ■■ 

TT4 

2C-3C. 62404 

A4/A4D 

1.00000 

TTl 

954.22444 







• FN 

4922.5401 

FG 

32763,6636 

FO 

27861.1235 

SFC 

1.51473 

WF 

7456.3203 

••• ' V/1 347,0717 

WG' 

FNl 

4922.5401 

FGl 

52783,6636 

N 

i.COOOO 

■ SFCl 

1.51473 

WIK 

443.5758 

VS 3075.8752 

CAPAS 

TT4 

2030.62404 

A4/A4D 

l.OOC-00 

TTl 

954.22444 






" . . . 

FN 

4075.0262 

FG 

31952.6335 

FD 

2'7886.6074 

SFC 

1.60494 

WF 

6541.7854 

W1 347.3392 

WG 

FNl 

4076.0262 

FGl 

31962.6335 

■ ' N 

1 » 00000 

SFCl 

1.60494 ■ 

■• ■ WIK 

443.9815 

VS 2998. 2v72 ' 

■ CAPAS ■ 

■TT4 

1980.62404 

A4/A40 

1.00000 

TTl 

954,22444 







FN 

3188.7874 

FG 

31101.1387 

FO 

27512.3512 

SFC 

1.76727 

WF 

5635.4473 

W1 347.7099 

WG 

FNl 

3133.7374 

FGl 

31101.1387 

N 

1.00300 

SFCl 

1.76727 

■ WIK 

444.3914 

V8 2916.9061 

■ CAPAS 

TT4 

1930.62404 

A47A4D 

l.OOOOC 

TU 

954.22444 








II 

Is 


■•t3^;9,5398 
'?32,.55 37 


34,9;. 4372 
r32r.36i7 ■ 


'349.2935 

32.112C 


349.1162 

31.7750 


349.0285 

■31,3046 


349.09C6 

■31,3132 


349,1409 

31.5144 


349.2116 

31.4918 


349*2798 

31.5783 



i7i-a 


Example 6 (Case 6) - This case is an exainple;-of*the use of the ZTHRUST 
control mode. When ZTHRUST is specified at an Altitude and Mach 
number as shown below, and with a part power- control operational, IPC = 3 
from Case 5, the program will interpolate the -.generated data at the 
specified ZTHRUST, altitude and Mach number 'ttx determine specific data 
used in engine acoustic analysis. *- 

SENPUT ' 

ICASE=6«ZTHRUST=70000, ,ALT=lOOO» 9EMACH=o3» 

$ 

Output is for standard +8°C atmospheric conditions in the format shown 
below: 


■ INPUT CASE 6 

PC 

3 ?5 

• 

50000000E+01 



' ■■ ALT .10000C00E + J4 

HACH 

« 30 jOCOOC'E+UO 

HPcXT 

.20000COOE+03 

■ ID033 

0 

• '■ WCQWl .lOOOOOOOE-03 

HB 

.lOOOOOOCiE + Cl 

CV • 

.95500COOE<-00 

ZTHRUST 

.70E+05 


■' NOZZLE 'PRESSURE RATIO LESS THAN 1.005 


OUTPUT ENGINE PARAMETERS CORRESPONOIMG TO THRUST VALUE OF 70000.0LSS, 


■“GASFLDW » ''697.4 L3S/SSC JET VELOCITY « 3614.37FT/SEC, JET AREA '» 10 .609SQ. FT. j- RAM DRAG « ' 7Q53.U3S. " 

'“NET THRUST » 70000. OLBS. TURBINE IN TEMP. •= 3060. OC.p. KDISF REOUCTIDN » O.OOODB. ATMOSPHERE « STD + 6.0 C. 

“SPECIFIC FUEL CONSUMPTION •= l,378CLSK7HR/L8F. TURBINE EXIT TEMP ■> 2104.5 R ' " ' 



page is 



